Validity and sensitivity of a human cranial finite element model: implications for comparative studies of biting performance

Biomechanics in anthropology

Biomechanics is the set of tools that explain organismal movement and mechanical behavior and links the organism to the physicality of the world. As such, biomechanics can relate behaviors and culture to the physicality of the organism. Scale is critical to biomechanical analyses, as the constitutive equations that matter differ depending on the scale of the question. Within anthropology, biomechanics has had a wide range of applications, from understanding how we and other primates evolved to understanding the effects of technologies, such as the atlatl, and the relationship between identity, society, culture, and medical interventions, such as prosthetics. Like any other model, there is great utility in biomechanical models, but models should be used primarily for hypothesis testing and not data generation except in the rare case where models can be robustly validated. The application of biomechanics within anthropology has been extensive, and holds great potential for the future.

Food texture and vitamin D influence mouse mandible form and molar roots

Industrialization influenced several facets of lifestyle, including softer nutrient-poor diets that contributed to vitamin D deficiency in post-industrzialized populations, with concomitantly increased dental problems. Here we simulated a post-industrialized diet in a mouse model to test the effects of diet texture and vitamin D level on mandible and third molar (M3) forms. Mice were raised on a soft diet with vitamin D (VitD) or without it (NoD), or on a hard diet with vitamin D. We hypothesized that a VitD/hard diet is optimal for normal mandible and tooth root form, as well as for timely M3 initiation. Subsets of adult NoD/soft and VitD/soft groups were bred to produce embryos that were micro-computed tomography (μCT) scanned to stage M3 development. M3 stage did not differ between embryos from mothers fed VitD and NoD diets, indicating that vitamin D does not affect timing of M3 onset. Sacrificed adult mice were μCT-scanned, their mandibles 3D-landmarked and M3 roots were measured. Principal component (PC) analysis described the largest proportion of mandible shape variance (PC1, 30.1%) related to diet texture, and nominal shape variance (PC2, 13.8%) related to vitamin D. Mice fed a soft diet had shorter, relatively narrower, and somewhat differently shaped mandibles that recapitulated findings in human populations. ANOVA and other multivariate tests found significantly wider M3 roots and larger root canals in mice fed a soft diet, with vitamin D having little effect. Altogether our experiments using a mouse model contribute new insights about how a post-industrial diet may influence human craniodental variation.

Effect of short implant crown-to-implant ratio on stress distribution in anisotropic bone with different osseointegration rates

OBJECTIVE

This study aimed to provide evidence for the clinical application of single short implants by establishing an anisotropic, three-dimensional (3D) finite element mandible model and simulating the effect of crown-to-implant ratio (CIR) on biomechanics around short implants with different osseointegration rates.

METHODS

Assuming that the bone is transversely isotropic by finite element method, we created four distinct models of implants for the mandibular first molar. Subsequently, axial and oblique forces were applied to the occlusal surface of these models. Ultimately, the Abaqus 2020 software was employed to compute various mechanical parameters, including the maximum von Mises stress, tensile stress, compressive stress, shear stress, displacement, and strains in the peri-implant bone tissue.

RESULTS

Upon establishing consistent osseointegration rates, the distribution of stress exhibited similarities across models with varying CIRs when subjected to vertical loads. However, when exposed to inclined loads, the maximum von Mises stress within the cortical bone escalated as the CIR heightened. Among both loading scenarios, notable escalation in the maximum von Mises stress occurred in the model featuring a CIR of 2.5 and an osseointegration rate of 25%. Conversely, other models displayed comparable strength. Notably, stress and strain values uniformly increased with augmented osseointegration across all models. Furthermore, an increase in osseointegration rate correlated with reduced maximum displacement for both cortical bone and implants.

CONCLUSIONS

After fixing osseointegration rates, the stress around shorter implants increased as the CIR increased under inclined loads. Thus, the effect of lateral forces should be considered when selecting shorter implants. Moreover, an implant failure risk was present in cases with a CIR ≥ 2.5 and low osseointegration rates. Additionally, the higher the osseointegration rate, the more readily the implant can achieve robust stability.

The biomechanical significance of the elongated rodent incisor root in the mandible during incision

Rodents are characterised by a distinctive masticatory apparatus which includes a single pair of enlarged and continually growing incisors. This morphology, termed diprotodonty, has also independently evolved in a number of other mammals, including the aye-aye. This study examined the functional significance of the internal "root" of the elongated rodent-like incisor. The mandibles of four rodents and an aye-aye were modelled to exhibit incrementally shorter incisor roots. Finite element analysis was used to predict stress and strain patterns across the jaw to determine whether the length of the incisor root contributes to the resistance of mechanical forces encountered in the mandible during incision. It was found that von Mises stresses increase in the region of the mandible local to where the incisor is removed, but that the stress distribution across the wider mandible is only minimally affected. Thus, the long internal incisor appears to play a small role in resisting bending forces close to the incisor alveolus, and may act with the arch-like mandibular shape to strengthen the mandible in this region. However, the impact across the whole mandible is relatively limited, suggesting the highly elongate incisor in diprotodont mammals may be principally driven by other factors such as rapid incisor wear.

Anthroengineering: an independent interdisciplinary field

In recent decades, funding agencies, institutes and professional bodies have recognized the profound benefits of transdisciplinarity in tackling targeted research questions. However, once questions are answered, the previously abundant support often dissolves. As such, the long-term benefits of these transdisciplinary approaches are never fully achieved. Over the last several decades, the integration of anthropology and engineering through inter- and multidisciplinary work has led to advances in fields such as design, human evolution and medical technologies. The lack of formal recognition, however, of this transdisciplinary approach as a unique entity rather than a useful tool or a subfield makes it difficult for researchers to establish laboratories, secure permanent jobs, fund long-term research programmes and train students in this approach. To facilitate the growth and development and witness the long-term benefits of this approach, we propose the integration of anthropology and engineering be recognized as a new, independent field known as anthroengineering. We present a working definition for anthroengineering and examples of how anthroengineering has been used. We discuss the necessity of recognizing anthroengineering as a unique field and explore potential novel applications. Finally, we discuss the future of anthroengineering, highlighting avenues for moving the field forward.

Effect of crown-to-implant ratio and crown height space on marginal bone stress: a finite element analysis

BACKGROUND

Crown-to-implant ratio and crown height space, associated with the use of short implants, have been related with marginal bone loss. However, it is unclear which of the two entities would play the most important role on the bone remodelling process. Using a finite element analysis, the present work aims to help clarifying how those two factors contribute for the stress generation at the marginal bone level. A numerical model (reference model), with a crown-to-implant ratio of 4, was double validated and submitted to a numerical calculation. Then, it was modified in two different ways: (a) by decreasing the prosthetic height obtaining crown-to-implant ratios of 3, 2.5 and 2 and (b) by increasing the implants length obtaining a crown-to-implant ratio of 2.08. The new models were also submitted to numerical calculations.

RESULTS

The reference model showed a marginal bone stress of 96.9 MPa. The increase in the implants' length did not show statistically significant differences in the marginal bone stress (p-value = 0.2364). The decrease in the prosthetic height was accompanied with a statistically significant decrease in the marginal bone stresses (p-value = 2.2e- 16).

CONCLUSIONS

The results represent a paradigm change as the crown height space appears to be more responsible for marginal bone stress than the high crown-to-implant ratios or the implants' length. New prosthetic designs should be attempted to decrease the stress generated at the marginal bone level.

Efficacy of shear strain gradients as an osteogenic stimulus

The question of which mechanical variables are responsible for inducing osteogenic activity is unresolved despite extensive experimental and theoretical investigation. Candidate variables include strain magnitude, loading frequency, the interaction of magnitude and frequency (strain rate), and strain gradients. An additional challenge is discerning the coordination of periosteal and endosteal expansion during growth, and whether this coordination (or lack thereof) is fully dependent or partially independent of the local mechanical environment. In this study, under the assumption that calculated stresses correspond to relative strain magnitudes, we specify alternative growth algorithms of bone cross-sectional size and geometry to explore skeletal growth under alternative scenarios of osteogenic activity that are tracking 1) an attractor stress, 2) local stress magnitude or 3) steepness of stress gradients. These developmental simulations are initiated from two initial geometries (symmetrical and asymmetrical ellipses) under a time-varying torsional load whose magnitude is proportional to body size growth in a model primate. In addition, we model endosteal expansion under three conditions hypothesized in the literature, in which endosteal expansion is 1) independent of the mechanical milieu, 2) completely dependent on the mechanical milieu, and 3) a "hybrid" model in which intrinsic biological (independent) growth is operative early but gives way to mechanically-sensitive (dependent) growth at later ages. Three variables were recorded over each growth simulation: the safety factor (ratio of yield stress to actual stress), an efficiency ratio (invested bone area per unit of stress), and proximity to an isostress condition (an optimal design criterion in which stress is invariant throughout the structure). The attractor stress algorithm produces the most "adapted" bones in terms of mechanical competence and economy of material. Localized osteogenic activity that is guided in direct proportion to stress magnitude produces competent bones but with variable adult geometries depending on conditions of endosteal expansion. Stress gradients also produce functional but relatively inefficient bones, with widely variable safety factors during growth and heterogeneous stress fields. If, in fact, the osteocyte network monitors strain gradients to generate osteogenic signals, the resulting morphology is competent but falls well short of an optimal mechanical solution.

Clarifying relationships between cranial form and function in tapirs, with implications for the dietary ecology of early hominins

Paleontologists and paleoanthropologists have long debated relationships between cranial morphology and diet in a broad diversity of organisms. While the presence of larger temporalis muscle attachment area (via the presence of sagittal crests) in carnivorans is correlated with durophagy (i.e. hard-object feeding), many primates with similar morphologies consume an array of tough and hard foods—complicating dietary inferences of early hominins. We posit that tapirs, large herbivorous mammals showing variable sagittal crest development across species, are ideal models for examining correlations between textural properties of food and sagittal crest morphology. Here, we integrate dietary data, dental microwear texture analysis, and finite element analysis to clarify the functional significance of the sagittal crest in tapirs. Most notably, pronounced sagittal crests are negatively correlated with hard-object feeding in extant, and several extinct, tapirs and can actually increase stress and strain energy. Collectively, these data suggest that musculature associated with pronounced sagittal crests—and accompanied increases in muscle volume—assists with the processing of tough food items in tapirs and may yield similar benefits in other mammals including early hominins.

Muscle-Bone Crosstalk in the Masticatory System: From Biomechanical to Molecular Interactions

The masticatory system is a complex and highly organized group of structures, including craniofacial bones (maxillae and mandible), muscles, teeth, joints, and neurovascular elements. While the musculoskeletal structures of the head and neck are known to have a different embryonic origin, morphology, biomechanical demands, and biochemical characteristics than the trunk and limbs, their particular molecular basis and cell biology have been much less explored. In the last decade, the concept of muscle-bone crosstalk has emerged, comprising both the loads generated during muscle contraction and a biochemical component through soluble molecules. Bone cells embedded in the mineralized tissue respond to the biomechanical input by releasing molecular factors that impact the homeostasis of the attaching skeletal muscle. In the same way, muscle-derived factors act as soluble signals that modulate the remodeling process of the underlying bones. This concept of muscle-bone crosstalk at a molecular level is particularly interesting in the mandible, due to its tight anatomical relationship with one of the biggest and strongest masticatory muscles, the masseter. However, despite the close physical and physiological interaction of both tissues for proper functioning, this topic has been poorly addressed. Here we present one of the most detailed reviews of the literature to date regarding the biomechanical and biochemical interaction between muscles and bones of the masticatory system, both during development and in physiological or pathological remodeling processes. Evidence related to how masticatory function shapes the craniofacial bones is discussed, and a proposal presented that the masticatory muscles and craniofacial bones serve as secretory tissues. We furthermore discuss our current findings of myokines-release from masseter muscle in physiological conditions, during functional adaptation or pathology, and their putative role as bone-modulators in the craniofacial system. Finally, we address the physiological implications of the crosstalk between muscles and bones in the masticatory system, analyzing pathologies or clinical procedures in which the alteration of one of them affects the homeostasis of the other. Unveiling the mechanisms of muscle-bone crosstalk in the masticatory system opens broad possibilities for understanding and treating temporomandibular disorders, which severely impair the quality of life, with a high cost for diagnosis and management.

The Mechanical Effect of the Periodontal Ligament on Bone Strain Regimes in a Validated Finite Element Model of a Macaque Mandible

The primary anatomical function of the periodontal ligament (PDL) is to attach teeth to their sockets. However, theoretical and constitutive mechanical models have proposed that during mastication the PDL redistributes local occlusal loads and reduces the jaw's resistance to torsional deformations. These hypotheses imply that accurately modeling the PDL's material properties and geometry in finite element analysis (FEA) is a prerequisite to obtaining precise strain and deformation data. Yet, many finite element studies of the human and non-human primate masticatory apparatus exclude the PDL or model it with simplicity, in part due to limitations in μCT/CT scan resolution and material property assignment. Previous studies testing the sensitivity of finite element models (FEMs) to the PDL have yielded contradictory results, however a major limitation of these studies is that FEMs were not validated against in vivo bone strain data. Hence, this study uses a validated and subject specific FEM to assess the effect of the PDL on strain and deformation regimes in the lower jaw of a rhesus macaque (Macaca mulatta) during simulated unilateral post-canine chewing. Our findings demonstrate that the presence of the PDL does influence local and global surface strain magnitudes (principal and shear) in the jaw. However, the PDL's effect is limited (diff. ~200-300 με) in areas away from the alveoli. Our results also show that varying the PDL's Young's Modulus within the range of published values (0.07-1750 MPa) has very little effect on global surface strains. These findings suggest that the mechanical importance of the PDL in FEMs of the mandible during chewing is dependent on the scope of the hypotheses being tested. If researchers are comparing strain gradients across species/taxa, the PDL may be excluded with minimal effect on results, but, if researchers are concerned with absolute strain values, sensitivity analysis is required.

An Overview of Modelling Craniosynostosis Using the Finite Element Method

Craniosynostosis is a medical condition caused by the early fusion of the cranial joint. The finite element method (FEM) is a computational technique that can answer a variety of "what if" questions in relation to the biomechanics of this condition. The aim of this study was to review the current literature that has used FEM to investigate the biomechanics of any aspect of craniosynostosis, being its development or its reconstruction. This review highlights that a relatively small number of studies (n = 10) has used FEM to investigate the biomechanics of craniosynostosis. Current studies set a good foundation for the future to take advantage of this method and optimize reconstruction of various forms of craniosynostosis.

Modelling human skull growth: a validated computational model

During the first year of life, the brain grows rapidly and the neurocranium increases to about 65% of its adult size. Our understanding of the relationship between the biomechanical forces, especially from the growing brain, the craniofacial soft tissue structures and the individual bone plates of the skull vault is still limited. This basic knowledge could help in the future planning of craniofacial surgical operations. The aim of this study was to develop a validated computational model of skull growth, based on the finite-element (FE) method, to help understand the biomechanics of skull growth. To do this, a two-step validation study was carried out. First, an in vitro physical three-dimensional printed model and an in silico FE model were created from the same micro-CT scan of an infant skull and loaded with forces from the growing brain from zero to two months of age. The results from the in vitro model validated the FE model before it was further developed to expand from 0 to 12 months of age. This second FE model was compared directly with in vivo clinical CT scans of infants without craniofacial conditions (n = 56). The various models were compared in terms of predicted skull width, length and circumference, while the overall shape was quantified using three-dimensional distance plots. Statistical analysis yielded no significant differences between the male skull models. All size measurements from the FE model versus the in vitro physical model were within 5%, with one exception showing a 7.6% difference. The FE model and in vivo data also correlated well, with the largest percentage difference in size being 8.3%. Overall, the FE model results matched well with both the in vitro and in vivo data. With further development and model refinement, this modelling method could be used to assist in preoperative planning of craniofacial surgery procedures and could help to reduce reoperation rates.

Predicting calvarial growth in normal and craniosynostotic mice using a computational approach

During postnatal calvarial growth the brain grows gradually and the overlying bones and sutures accommodate that growth until the later juvenile stages. The whole process is coordinated through a complex series of biological, chemical and perhaps mechanical signals between various elements of the craniofacial system. The aim of this study was to investigate to what extent a computational model can accurately predict the calvarial growth in wild-type (WT) and mutant type (MT) Fgfr2C342Y/+ mice displaying bicoronal suture fusion. A series of morphological studies were carried out to quantify the calvarial growth at P3, P10 and P20 in both mouse types. MicroCT images of a P3 specimen were used to develop a finite element model of skull growth to predict the calvarial shape of WT and MT mice at P10. Sensitivity tests were performed and the results compared with ex vivo P10 data. Although the models were sensitive to the choice of input parameters, they predicted the overall skull growth in the WT and MT mice. The models also captured the difference between the ex vivoWT and MT mice. This modelling approach has the potential to be translated to human skull growth and to enhance our understanding of the different reconstruction methods used to manage clinically the different forms of craniosynostosis, and in the long term possibly reduce the number of re-operations in children displaying this condition and thereby enhance their quality of life.

The jaw is a second-class lever in Pedetes capensis (Rodentia: Pedetidae)

The mammalian jaw is often modelled as a third-class lever for the purposes of biomechanical analyses, owing to the position of the resultant muscle force between the jaw joint and the teeth. However, it has been proposed that in some rodents the jaws operate as a second-class lever during distal molar bites, owing to the rostral position of the masticatory musculature. In particular, the infraorbital portion of the zygomatico-mandibularis (IOZM) has been suggested to be of major importance in converting the masticatory system from a third-class to a second-class lever. The presence of the IOZM is diagnostic of the hystricomorph rodents, and is particularly well-developed in Pedetes capensis, the South African springhare. In this study, finite element analysis (FEA) was used to assess the lever mechanics of the springhare masticatory system, and to determine the function of the IOZM. An FE model of the skull of P. capensis was constructed and loaded with all masticatory muscles, and then solved for biting at each tooth in turn. Further load cases were created in which each masticatory muscle was removed in turn. The analyses showed that the mechanical advantage of the springhare jaws was above one at all molar bites and very close to one during the premolar bite. Removing the IOZM or masseter caused a drop in mechanical advantage at all bites, but affected strain patterns and cranial deformation very little. Removing the ZM had only a small effect on mechanical advantage, but produced a substantial reduction in strain and deformation across the skull. It was concluded that the masticatory system of P. capensis acts as a second class lever during bites along almost the entire cheek tooth row. The IOZM is clearly a major contributor to this effect, but the masseter also has a part to play. The benefit of the IOZM is that it adds force without substantially contributing to strain or deformation of the skull. This may help explain why the hystricomorphous morphology has evolved multiple times independently within Rodentia.

Form-function relationships in dragonfly mandibles under an evolutionary perspective

Functional requirements may constrain phenotypic diversification or foster it. For insect mouthparts, the quantification of the relationship between shape and function in an evolutionary framework remained largely unexplored. Here, the question of a functional influence on phenotypic diversification for dragonfly mandibles is assessed with a large-scale biomechanical analysis covering nearly all anisopteran families, using finite element analysis in combination with geometric morphometrics. A constraining effect of phylogeny could be found for shape, the mandibular mechanical advantage (MA), and certain mechanical joint parameters, while stresses and strains, the majority of joint parameters and size are influenced by shared ancestry. Furthermore, joint mechanics are correlated with neither strain nor mandibular MA and size effects have virtually play no role for shape or mechanical variation. The presence of mandibular strengthening ridges shows no phylogenetic signal except for one ridge peculiar to Libelluloidea, and ridge presence is also not correlated with each other. The results suggest that functional traits are more variable at this taxonomic level and that they are not influenced by shared ancestry. At the same time, the results contradict the widespread idea that mandibular morphology mainly reflects functional demands at least at this taxonomic level. The varying functional factors rather lead to the same mandibular performance as expressed by the MA, which suggests a many-to-one mapping of the investigated parameters onto the same narrow mandibular performance space.

Human feeding biomechanics: performance, variation, and functional constraints

The evolution of the modern human (Homo sapiens) cranium is characterized by a reduction in the size of the feeding system, including reductions in the size of the facial skeleton, postcanine teeth, and the muscles involved in biting and chewing. The conventional view hypothesizes that gracilization of the human feeding system is related to a shift toward eating foods that were less mechanically challenging to consume and/or foods that were processed using tools before being ingested. This hypothesis predicts that human feeding systems should not be well-configured to produce forceful bites and that the cranium should be structurally weak. An alternate hypothesis, based on the observation that humans have mechanically efficient jaw adductors, states that the modern human face is adapted to generate and withstand high biting forces. We used finite element analysis (FEA) to test two opposing mechanical hypotheses: that compared to our closest living relative, chimpanzees (Pan troglodytes), the modern human craniofacial skeleton is (1) less well configured, or (2) better configured to generate and withstand high magnitude bite forces. We considered intraspecific variation in our examination of human feeding biomechanics by examining a sample of geographically diverse crania that differed notably in shape. We found that our biomechanical models of human crania had broadly similar mechanical behavior despite their shape variation and were, on average, less structurally stiff than the crania of chimpanzees during unilateral biting when loaded with physiologically-scaled muscle loads. Our results also show that modern humans are efficient producers of bite force, consistent with previous analyses. However, highly tensile reaction forces were generated at the working (biting) side jaw joint during unilateral molar bites in which the chewing muscles were recruited with bilateral symmetry. In life, such a configuration would have increased the risk of joint dislocation and constrained the maximum recruitment levels of the masticatory muscles on the balancing (non-biting) side of the head. Our results do not necessarily conflict with the hypothesis that anterior tooth (incisors, canines, premolars) biting could have been selectively important in humans, although the reduced size of the premolars in humans has been shown to increase the risk of tooth crown fracture. We interpret our results to suggest that human craniofacial evolution was probably not driven by selection for high magnitude unilateral biting, and that increased masticatory muscle efficiency in humans is likely to be a secondary byproduct of selection for some function unrelated to forceful biting behaviors. These results are consistent with the hypothesis that a shift to softer foods and/or the innovation of pre-oral food processing techniques relaxed selective pressures maintaining craniofacial features that favor forceful biting and chewing behaviors, leading to the characteristically small and gracile faces of modern humans.

The Effect of Varying Jaw-elevator Muscle Forces on a Finite Element Model of a Human Cranium

Finite element analyses simulating masticatory system loading are increasingly undertaken in primates, hominin fossils and modern humans. Simplifications of models and loadcases are often required given the limits of data and technology. One such area of uncertainty concerns the forces applied to cranial models and their sensitivity to variations in these forces. We assessed the effect of varying force magnitudes among jaw-elevator muscles applied to a finite element model of a human cranium. The model was loaded to simulate incisor and molar bites using different combinations of muscle forces. Symmetric, asymmetric, homogeneous, and heterogeneous muscle activations were simulated by scaling maximal forces. The effects were compared with respect to strain distribution (i.e., modes of deformation) and magnitudes; bite forces and temporomandibular joint (TMJ) reaction forces. Predicted modes of deformation, strain magnitudes and bite forces were directly proportional to total applied muscle force and relatively insensitive to the degree of heterogeneity of muscle activation. However, TMJ reaction forces and mandibular fossa strains decrease and increase on the balancing and working sides according to the degree of asymmetry of loading. These results indicate that when modes, rather than magnitudes, of facial deformation are of interest, errors in applied muscle forces have limited effects. However the degree of asymmetric loading does impact on TMJ reaction forces and mandibular fossa strains. These findings are of particular interest in relation to studies of skeletal and fossil material, where muscle data are not available and estimation of muscle forces from skeletal proxies is prone to error. Anat Rec, 299:828-839, 2016. © 2016 Wiley Periodicals, Inc.