The functional role of the rabbit digastric muscle during mastication.

Muscle spindle abundance is highly variable in vertebrates, but the functional determinants of this variation are unclear. Recent work has shown that human leg muscles with the lowest abundance of muscle spindles primarily function to lengthen and absorb energy, while muscles with a greater spindle abundance perform active-stretch-shorten cycles with no net work, suggesting that muscle spindle abundance may be underpinned by muscle function. Compared with other mammalian muscles, the digastric muscle contains the lowest abundance of muscle spindles and, therefore, might be expected to generate substantial negative work. However, it is widely hypothesised that as a jaw-opener (anatomically) the digastric muscle would primarily function to depress the jaw, and consequently do positive work. Through a combination of X-ray reconstruction of moving morphology (XROMM), electromyography and fluoromicrometry, we characterised the 3D kinematics of the jaw and digastric muscle during feeding in rabbits. Subsequently, the work loop technique was used to simulate in vivo muscle behaviour in situ, enabling muscle force to be quantified in relation to muscle strain and hence determine the muscle's function during mastication. When functioning on either the working or balancing side, the digastric muscle generates a large amount of positive work during jaw opening, and a large amount of negative work during jaw closing, on average producing a relatively small amount of net negative work. Our data therefore further support the hypothesis that muscle spindle abundance is linked to muscle function; specifically, muscles that absorb a relatively large amount of negative work have a low spindle abundance.

The role of miniaturization in the evolution of the mammalian jaw and middle ear.

The evolution of the mammalian jaw is one of the most important innovations in vertebrate history, and underpins the exceptional radiation and diversification of mammals over the last 220 million years1,2. In particular, the transformation of the mandible into a single tooth-bearing bone and the emergence of a novel jaw joint-while incorporating some of the ancestral jaw bones into the mammalian middle ear-is often cited as a classic example of the repurposing of morphological structures3,4. Although it is remarkably well-documented in the fossil record, the evolution of the mammalian jaw still poses the paradox of how the bones of the ancestral jaw joint could function both as a joint hinge for powerful load-bearing mastication and as a mandibular middle ear that was delicate enough for hearing. Here we use digital reconstructions, computational modelling and biomechanical analyses to demonstrate that the miniaturization of the early mammalian jaw was the primary driver for the transformation of the jaw joint. We show that there is no evidence for a concurrent reduction in jaw-joint stress and increase in bite force in key non-mammaliaform taxa in the cynodont-mammaliaform transition, as previously thought5-8. Although a shift in the recruitment of the jaw musculature occurred during the evolution of modern mammals, the optimization of mandibular function to increase bite force while reducing joint loads did not occur until after the emergence of the neomorphic mammalian jaw joint. This suggests that miniaturization provided a selective regime for the evolution of the mammalian jaw joint, followed by the integration of the postdentary bones into the mammalian middle ear.

Bite force and cranial bone strain in four species of lizards.

In vivo bone strain data provide direct evidence of strain patterns in the cranium during biting. Compared with those in mammals, in vivo bone strains in lizard skulls are poorly documented. This paper presents strain data from the skulls of Anolis equestris, Gekko gecko, Iguana iguana and Salvator merianae during transducer biting. Analysis of variance was used to investigate effects of bite force, bite point, diet, cranial morphology and cranial kinesis on strain magnitude. Within individuals, the most consistent determinants of variance in bone strain magnitude were gauge location and bite point, with the importance of bite force varying between individuals. Inter-site variance in strain magnitude - strain gradient - was present in all individuals and varied with bite point. Between individuals within species, variance in strain magnitude was driven primarily by variation in bite force, not gauge location or bite point, suggesting that inter-individual variation in patterns of strain magnitude is minimal. Between species, variation in strain magnitude was significantly impacted by bite force and species membership, as well as by interactions between gauge location, species and bite point. Independent of bite force, species differences in cranial strain magnitude may reflect selection for different cranial morphology in relation to feeding function, but what these performance criteria are is not clear. The relatively low strain magnitudes in Iguana and Uromastyx compared with those in other lizards may be related to their herbivorous diet. Cranial kinesis and the presence or absence of postorbital and supratemporal bars are not important determinants of inter-specific variation in strain magnitude.

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

Finite element analysis (FEA) is a modelling technique increasingly used in anatomical studies investigating skeletal form and function. In the case of the cranium this approach has been applied to both living and fossil taxa to (for example) investigate how form relates to function or infer diet or behaviour. However, FE models of complex musculoskeletal structures always rely on simplified representations because it is impossible completely to image and represent every detail of skeletal morphology, variations in material properties and the complexities of loading at all spatial and temporal scales. The effects of necessary simplifications merit investigation. To this end, this study focuses on one aspect, model geometry, which is particularly pertinent to fossil material where taphonomic processes often destroy the finer details of anatomy or in models built from clinical CTs where the resolution is limited and anatomical details are lost. We manipulated the details of a finite element (FE) model of an adult human male cranium and examined the impact on model performance. First, using digital speckle interferometry, we directly measured strains from the infraorbital region and frontal process of the maxilla of the physical cranium under simplified loading conditions, simulating incisor biting. These measured strains were then compared with predicted values from FE models with simplified geometries that included modifications to model resolution, and how cancellous bone and the thin bones of the circum-nasal and maxillary regions were represented. Distributions of regions of relatively high and low principal strains and principal strain vector magnitudes and directions, predicted by the most detailed FE model, are generally similar to those achieved in vitro. Representing cancellous bone as solid cortical bone lowers strain magnitudes substantially but the mode of deformation of the FE model is relatively constant. In contrast, omitting thin plates of bone in the circum-nasal region affects both mode and magnitude of deformation. Our findings provide a useful frame of reference with regard to the effects of simplifications on the performance of FE models of the cranium and call for caution in the interpretation and comparison of FEA results.

In vivo cranial bone strain and bite force in the agamid lizard Uromastyx geyri.

In vivo bone strain data are the most direct evidence of deformation and strain regimes in the vertebrate cranium during feeding and can provide important insights into skull morphology. Strain data have been collected during feeding across a wide range of mammals; in contrast, in vivo cranial bone strain data have been collected from few sauropsid taxa. Here we present bone strain data recorded from the jugal of the herbivorous agamid lizard Uromastyx geyri along with simultaneously recorded bite force. Principal and shear strain magnitudes in Uromastyx geyri were lower than cranial bone strains recorded in Alligator mississippiensis, but higher than those reported from herbivorous mammals. Our results suggest that variations in principal strain orientations in the facial skeleton are largely due to differences in feeding behavior and bite location, whereas food type has little impact on strain orientations. Furthermore, mean principal strain orientations differ between male and female Uromastyx during feeding, potentially because of sexual dimorphism in skull morphology.

Functional reorganisation of the cranial skeleton during the cynodont-mammaliaform transition.

Skeletal simplification occurred in multiple vertebrate clades over the last 500 million years, including the evolution from premammalian cynodonts to mammals. This transition is characterised by the loss and reduction of cranial bones, the emergence of a novel jaw joint, and the rearrangement of the jaw musculature. These modifications have long been hypothesised to increase skull strength and efficiency during feeding. Here, we combine digital reconstruction and biomechanical modelling to show that there is no evidence for an increase in cranial strength and biomechanical performance. Our analyses demonstrate the selective functional reorganisation of the cranial skeleton, leading to reduced stresses in the braincase and the skull roof but increased stresses in the zygomatic region through this transition. This cranial functional reorganisation, reduction in mechanical advantage, and overall miniaturisation in body size are linked with a dietary specialisation to insectivory, permitting the subsequent morphological and ecological diversification of the mammalian lineage.

Developing a musculoskeletal model of the primate skull: predicting muscle activations, bite force, and joint reaction forces using multibody dynamics analysis and advanced optimisation methods.

An accurate, dynamic, functional model of the skull that can be used to predict muscle forces, bite forces, and joint reaction forces would have many uses across a broad range of disciplines. One major issue however with musculoskeletal analyses is that of muscle activation pattern indeterminacy. A very large number of possible muscle force combinations will satisfy a particular functional task. This makes predicting physiological muscle recruitment patterns difficult. Here we describe in detail the process of development of a complex multibody computer model of a primate skull (Macaca fascicularis), that aims to predict muscle recruitment patterns during biting. Using optimisation criteria based on minimisation of muscle stress we predict working to balancing side muscle force ratios, peak bite forces, and joint reaction forces during unilateral biting. Validation of such models is problematic; however we have shown comparable working to balancing muscle activity and TMJ reaction ratios during biting to those observed in vivo and that peak predicted bite forces compare well to published experimental data. To our knowledge the complexity of the musculoskeletal model is greater than any previously reported for a primate. This complexity, when compared to more simple representations provides more nuanced insights into the functioning of masticatory muscles. Thus, we have shown muscle activity to vary throughout individual muscle groups, which enables them to function optimally during specific masticatory tasks. This model will be utilised in future studies into the functioning of the masticatory apparatus.

Biomechanical evaluation of the unilateral crossbite on the asymmetrical development of the craniofacial complex. A mechano-morphological approach.

BACKGROUND AND OBJECTIVE: The occlusion effect on the craniofacial development is a controversial topic that has attracted the interest of many researchers but that remains unclear, mainly due to the difficulties on measure its mechanical response experimentally. This mechano-morphological relationship of the craniofacial growth is often explained by the periosteal and capsular matrices of the functional matrix hypothesis (FMH); however, its outcomes have not been analytically demonstrated yet. This computational study aims, therefore, to analytically demonstrate the mechano-morphological relationship in the craniofacial development of children with unilateral crossbite (UXB) using the finite element (FE) method. METHODS: The craniofacial complex asymmetry of ten children, five of whom exhibit UXB, was 3D-analysed and compared with the biomechanical response computed from a FE analysis of each patient's occlusion. Due to the complexity of the geometry and the multitude of contacts involved, the inherent limitations of the model were evaluated by comparing computed occlusal patterns with those recorded by an occlusal analysis on 3D printed copies. RESULTS: Comparison's outcomes proved the reliability of our models with just a deviation error below 6% between both approaches. Out of validation process, computational results showed that the significant elongation of mandibular branch in the contralateral side could be related to the mandibular shift and increase of thickness on the crossed side, and particularly of the posterior region. These morphological changes could be associated with periodontal overpressure (>4.7 kPa) and mandibular over deformation (0.002 ε) in that side, in agreement with the periosteal matrix's principles. Furthermore, the maxilla's transversal narrowing and the elevation of the maxillary and zygomatic regions on the crossed side were statistically demonstrated and seem to be related with their respective micro displacements at occlusion, as accounted by their specific capsule matrices. Our results were consistent with those reported clinically and demonstrated analytically the mechano-morphological relationship of children's craniofacial development based on the FMH's functional matrices. CONCLUSIONS: This study is a first step in the understanding of the occlusion's effect on the craniofacial development by computational methods. Our approach could help future engineers, researchers and clinicians to understand better the aetiology of some dental malocclusions and functional disorders improve the diagnosis or even predict the craniofacial development.

Why do humans have chins? Testing the mechanical significance of modern human symphyseal morphology with finite element analysis.

The modern human mandibular symphysis differs from those of all other primates in being vertically orientated and possessing a chin, but the functional significance of this unique morphology is not well understood. Some hypotheses propose that it is an adaptation to specific loads occurring during masticatory function. This study uses finite element analysis to examine these symphyseal loads in a model of a modern human mandible. By modifying the symphyseal cross-sectional form, the mechanical significance of the presence of the chin and symphyseal orientation is tested, and modern human and Neanderthal symphyseal cross-sections are compared with regard to their ability to withstand different loads. The results show that changes in symphyseal form have profound effects on the strains. The presence of a chin leads to lower symphyseal strains overall, whereas a vertical orientation of the symphysis results in higher strains under wishboning, but not under vertical bending in the coronal plane and dorsoventral shear. Compared to Neanderthals, the modern human symphysis shows higher strains during dorsoventral shear and wishboning, but is as effective as the Neanderthal symphysis in resisting vertical bending in the coronal plane and the loads resulting from simulated incision and unilateral molar biting. In general, the results of this study corroborate prior hypotheses about the mechanical effects of the human chin and vertical symphyseal orientation and support the idea that the relative importance of wishboning and vertical bending in the coronal plane might have played a role in the evolution of modern human symphyseal morphology.

Back to the bones: do muscle area assessment techniques predict functional evolution across a macroevolutionary radiation?

Measures of attachment or accommodation area on the skeleton are a popular means of rapidly generating estimates of muscle proportions and functional performance for use in large-scale macroevolutionary studies. Herein, we provide the first evaluation of the accuracy of these muscle area assessment (MAA) techniques for estimating muscle proportions, force outputs and bone loading in a comparative macroevolutionary context using the rodent masticatory system as a case study. We find that MAA approaches perform poorly, yielding large absolute errors in muscle properties, bite force and particularly bone stress. Perhaps more fundamentally, these methods regularly fail to correctly capture many qualitative differences between rodent morphotypes, particularly in stress patterns in finite-element models. Our findings cast doubts on the validity of these approaches as means to provide input data for biomechanical models applied to understand functional transitions in the fossil record, and perhaps even in taxon-rich statistical models that examine broad-scale macroevolutionary patterns. We suggest that future work should go back to the bones to test if correlations between attachment area and muscle size within homologous muscles across a large number of species yield strong predictive relationships that could be used to deliver more accurate predictions for macroevolutionary and functional studies.

Computational biomechanical modelling of the rabbit cranium during mastication.

Although a functional relationship between bone structure and mastication has been shown in some regions of the rabbit skull, the biomechanics of the whole cranium during mastication have yet to be fully explored. In terms of cranial biomechanics, the rabbit is a particularly interesting species due to its uniquely fenestrated rostrum, the mechanical function of which is debated. In addition, the rabbit processes food through incisor and molar biting within a single bite cycle, and the potential influence of these bite modes on skull biomechanics remains unknown. This study combined the in silico methods of multi-body dynamics and finite element analysis to compute musculoskeletal forces associated with a range of incisor and molar biting, and to predict the associated strains. The results show that the majority of the cranium, including the fenestrated rostrum, transmits masticatory strains. The peak strains generated over all bites were found to be attributed to both incisor and molar biting. This could be a consequence of a skull shape adapted to promote an even strain distribution for a combination of infrequent incisor bites and cyclic molar bites. However, some regions, such as the supraorbital process, experienced low peak strain for all masticatory loads considered, suggesting such regions are not designed to resist masticatory forces.

Towards an early 3D-diagnosis of craniofacial asymmetry by computing the accurate midplane: A PCA-based method.

BACKGROUND AND OBJECTIVE: Craniofacial asymmetry is a common growth disorder often caused by unilateral chewing. Although an early orthodontic treatment would avoid surgical procedures later in life, the uncertainty of defining the accurate sagittal midplane potentially leads to misdiagnosis and therefore inaccurate orthodontic treatment plans. This novel study aims to 3D-diagnose craniofacial complex malformations in children with unilateral crossbite (UXB) considering a midplane which compensates the asymmetric morphology. METHODS: The sagittal midplane of 20 children, fifteen of whom exhibited UXB, was computed by a PCA-based method which compensates the asymmetry mirroring the 3D models obtained from cone-beam computed tomography data. Once determined, one side of the data was mirrored using the computed midplane to visualize the malformations on the hard and soft tissues by 3D-computing the distances between both halves. Additionally, 31 skull's landmarks were manually placed in each model to study the principal variation modes and the significant differences in the group of subjects with and without UXB through PCA and Mann-Whitney U test analyses respectively. RESULTS: Morphological 3D-analysis showed pronounced deformities and aesthetic implications for patients with severe asymmetry (jaw deviation > 0.8 mm) in whole craniofacial system, while initial signs of asymmetry were found indistinctly in the mandible or maxilla. We detected significant (p < 0.05) malformations for example in mandibular ramus length (0.0086), maxillary palate width (0.0481) and condylar head width (0.0408). Craniofacial malformations increased the landmarks' variability in the group of patients with UXB over the control group requiring 8 variation modes more to define 99% of the sample' variability. CONCLUSIONS: Our findings demonstrated the viability of early diagnosis of craniofacial asymmetry through computing the accurate sagittal midplane which compensates the individual's asymmetrical morphology. Furthermore, this study provides important computational insights into the determination of craniofacial deformities which are caused by UXB, following some empirical findings of previous clinical studies. Hence, this computational approach can be useful for the development of new software in craniofacial surgery or for its use in biomedical research and clinical practice.

Mechanical adaptation of trabecular bone morphology in the mammalian mandible.

Alveolar bone, together with the underlying trabecular bone, fulfils an important role in providing structural support against masticatory forces. Diseases such as osteoporosis or periodontitis cause alveolar bone resorption which weakens this structural support and is a major cause of tooth loss. However, the functional relationship between alveolar bone remodelling within the molar region and masticatory forces is not well understood. This study investigated this relationship by comparing mammalian species with different diets and functional loading (Felis catus, Cercocebus atys, Homo sapiens, Sus scrofa, Oryctolagus cuniculus, Ovis aries). We performed histomorphometric analyses of trabecular bone morphology (bone volume fraction, trabecular thickness and trabecular spacing) and quantified the variation of bone and tooth root volumes along the tooth row. A principal component analysis and non-parametric MANOVA showed statistically significant differences in trabecular bone morphology between species with contrasting functional loading, but these differences were not seen in sub-adult specimens. Our results support a strong, but complex link between masticatory function and trabecular bone morphology. Further understanding of a potential functional relationship could aid the diagnosis and treatment of mandibular diseases causing alveolar bone resorption, and guide the design and evaluation of dental implants.

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.

The biomechanical role of the chondrocranium and sutures in a lizard cranium.

The role of soft tissues in skull biomechanics remains poorly understood. Not least, the chondrocranium, the portion of the braincase which persists as cartilage with varying degrees of mineralization. It also remains commonplace to overlook the biomechanical role of sutures despite evidence that they alter strain distribution. Here, we examine the role of both the sutures and the chondrocranium in the South American tegu lizard Salvator merianae We use multi-body dynamics analysis (MDA) to provide realistic loading conditions for anterior and posterior unilateral biting and a detailed finite element model to examine strain magnitude and distribution. We find that strains within the chondrocranium are greatest during anterior biting and are primarily tensile; also that strain within the cranium is not greatly reduced by the presence of the chondrocranium unless it is given the same material properties as bone. This result contradicts previous suggestions that the anterior portion (the nasal septum) acts as a supporting structure. Inclusion of sutures to the cranium model not only increases overall strain magnitudes but also leads to a more complex distribution of tension and compression rather than that of a beam under sagittal bending.

Inclusion of periodontal ligament fibres in mandibular finite element models leads to an increase in alveolar bone strains.

Alveolar bone remodelling is vital for the success of dental implants and orthodontic treatments. However, the underlying biomechanical mechanisms, in particular the function of the periodontal ligament (PDL) in bone loading and remodelling, are not well understood. The PDL is a soft fibrous connective tissue that joins the tooth root to the alveolar bone and plays a critical role in the transmission of loads from the tooth to the surrounding bone. However, due to its complex structure, small size and location within the tooth socket it is difficult to study in vivo. Finite element analysis (FEA) is an ideal tool with which to investigate the role of the PDL, however inclusion of the PDL in FE models is complex and time consuming, therefore consideration must be given to how it is included. The aim of this study was to investigate the effects of including the PDL and its fibrous structure in mandibular finite element models. A high-resolution model of a human molar region was created from micro-computed tomography scans. This is the first time that the fibrous structure of the PDL has been included in a model with realistic tooth and bone geometry. The results show that omission of the PDL creates a more rigid model, reducing the strains observed in the mandibular corpus which are of interest when considering mandibular functional morphology. How the PDL is modelled also affects the strains. The inclusion of PDL fibres alters the strains in the mandibular bone, increasing the strains in the tooth socket compared to PDL modelled without fibres. As strains in the alveolar bone are thought to play a key role in bone remodelling during orthodontic tooth movement, future FE analyses aimed at improving our understanding and management of orthodontic treatment should include the fibrous structure of the PDL.

Morphological evolution of the mammalian jaw adductor complex.

The evolution of the mammalian jaw during the transition from non-mammalian synapsids to crown mammals is a key event in vertebrate history and characterised by the gradual reduction of its individual bones into a single element and the concomitant transformation of the jaw joint and its incorporation into the middle ear complex. This osteological transformation is accompanied by a rearrangement and modification of the jaw adductor musculature, which is thought to have allowed the evolution of a more-efficient masticatory system in comparison to the plesiomorphic synapsid condition. While osteological characters relating to this transition are well documented in the fossil record, the exact arrangement and modifications of the individual adductor muscles during the cynodont-mammaliaform transition have been debated for nearly a century. We review the existing knowledge about the musculoskeletal evolution of the mammalian jaw adductor complex and evaluate previous hypotheses in the light of recently documented fossils that represent new specimens of existing species, which are of central importance to the mammalian origins debate. By employing computed tomography (CT) and digital reconstruction techniques to create three-dimensional models of the jaw adductor musculature in a number of representative non-mammalian cynodonts and mammaliaforms, we provide an updated perspective on mammalian jaw muscle evolution. As an emerging consensus, current evidence suggests that the mammal-like division of the jaw adductor musculature (into deep and superficial components of the m. masseter, the m. temporalis and the m. pterygoideus) was completed in Eucynodontia. The arrangement of the jaw adductor musculature in a mammalian fashion, with the m. pterygoideus group inserting on the dentary was completed in basal Mammaliaformes as suggested by the muscle reconstruction of Morganucodon oehleri. Consequently, transformation of the jaw adductor musculature from the ancestral ('reptilian') to the mammalian condition must have preceded the emergence of Mammalia and the full formation of the mammalian jaw joint. This suggests that the modification of the jaw adductor system played a pivotal role in the functional morphology and biomechanical stability of the jaw joint.

The biomechanical function of periodontal ligament fibres in orthodontic tooth movement.

Orthodontic tooth movement occurs as a result of resorption and formation of the alveolar bone due to an applied load, but the stimulus responsible for triggering orthodontic tooth movement remains the subject of debate. It has been suggested that the periodontal ligament (PDL) plays a key role. However, the mechanical function of the PDL in orthodontic tooth movement is not well understood as most mechanical models of the PDL to date have ignored the fibrous structure of the PDL. In this study we use finite element (FE) analysis to investigate the strains in the alveolar bone due to occlusal and orthodontic loads when PDL is modelled as a fibrous structure as compared to modelling PDL as a layer of solid material. The results show that the tension-only nature of the fibres essentially suspends the tooth in the tooth socket and their inclusion in FE models makes a significant difference to both the magnitude and distribution of strains produced in the surrounding bone. The results indicate that the PDL fibres have a very important role in load transfer between the teeth and alveolar bone and should be considered in FE studies investigating the biomechanics of orthodontic tooth movement.

Masticatory biomechanics in the rabbit: a multi-body dynamics analysis.

Multi-body dynamics is a powerful engineering tool which is becoming increasingly popular for the simulation and analysis of skull biomechanics. This paper presents the first application of multi-body dynamics to analyse the biomechanics of the rabbit skull. A model has been constructed through the combination of manual dissection and three-dimensional imaging techniques (magnetic resonance imaging and micro-computed tomography). Individual muscles are represented with multiple layers, thus more accurately modelling muscle fibres with complex lines of action. Model validity was sought through comparing experimentally measured maximum incisor bite forces with those predicted by the model. Simulations of molar biting highlighted the ability of the masticatory system to alter recruitment of two muscle groups, in order to generate shearing or crushing movements. Molar shearing is capable of processing a food bolus in all three orthogonal directions, whereas molar crushing and incisor biting are predominately directed vertically. Simulations also show that the masticatory system is adapted to process foods through several cycles with low muscle activations, presumably in order to prevent rapidly fatiguing fast fibres during repeated chewing cycles. Our study demonstrates the usefulness of a validated multi-body dynamics model for investigating feeding biomechanics in the rabbit, and shows the potential for complementing and eventually reducing in vivo experiments.

The importance of accurate muscle modelling for biomechanical analyses: a case study with a lizard skull.

Computer-based simulation techniques such as multi-body dynamics analysis are becoming increasingly popular in the field of skull mechanics. Multi-body models can be used for studying the relationships between skull architecture, muscle morphology and feeding performance. However, to be confident in the modelling results, models need to be validated against experimental data, and the effects of uncertainties or inaccuracies in the chosen model attributes need to be assessed with sensitivity analyses. Here, we compare the bite forces predicted by a multi-body model of a lizard (Tupinambis merianae) with in vivo measurements, using anatomical data collected from the same specimen. This subject-specific model predicts bite forces that are very close to the in vivo measurements and also shows a consistent increase in bite force as the bite position is moved posteriorly on the jaw. However, the model is very sensitive to changes in muscle attributes such as fibre length, intrinsic muscle strength and force orientation, with bite force predictions varying considerably when these three variables are altered. We conclude that accurate muscle measurements are crucial to building realistic multi-body models and that subject-specific data should be used whenever possible.