The importance of craniofacial sutures in biomechanical finite element models of the domestic pig

Regional variation of the cortical and trabecular bone material properties in the rabbit skull

The material properties of some bones are known to vary with anatomical location, orientation and position within the bone (e.g., cortical and trabecular bone). Details of the heterogeneity and anisotropy of bone is an important consideration for biomechanical studies that apply techniques such as finite element analysis, as the outcomes will be influenced by the choice of material properties used. Datasets detailing the regional variation of material properties in the bones of the skull are sparse, leaving many finite element analyses of skulls no choice but to employ homogeneous, isotropic material properties, often using data from a different species to the one under investigation. Due to the growing significance of investigating the cranial biomechanics of the rabbit in basic science and clinical research, this study used nanoindentation to measure the elastic modulus of cortical and trabecular bone throughout the skull. The elastic moduli of cortical bone measured in the mediolateral and ventrodorsal direction were found to decrease posteriorly through the skull, while it was evenly distributed when measured in the anteroposterior direction. Furthermore, statistical tests showed that the variation of elastic moduli between separate regions (anterior, middle and posterior) of the skull were significantly different in cortical bone, but was not in trabecular bone. Elastic moduli measured in different orthotropic planes were also significantly different, with the moduli measured in the mediolateral direction consistently lower than that measured in either the anteroposterior or ventrodorsal direction. These findings demonstrate the significance of regional and directional variation in cortical bone elastic modulus, and therefore material properties in finite element models of the skull, particularly those of the rabbit, should consider the heterogeneous and orthotropic properties of skull bone when possible.

Assessment of the mechanical role of cranial sutures in the mammalian skull: Computational biomechanical modelling of the rat skull

Cranial sutures are fibrocellular joints between the skull bones that are progressively replaced with bone throughout ontogeny, facilitating growth and cranial shape change. This transition from soft tissue to bone is reflected in the biomechanical properties of the craniofacial complex. However, the mechanical significance of cranial sutures has only been explored at a few localised areas within the mammalian skull, and as such our understanding of suture function in overall skull biomechanics is still limited. Here, we sought to determine how the overall strain environment is affected by the complex network of cranial sutures in the mammal skull. We combined two computational biomechanical methods, multibody dynamics analysis and finite element analysis, to simulate biting in a rat skull and compared models with and without cranial sutures. Our results show that including complex sutures in the rat model does not substantially change overall strain gradients across the cranium, particularly strain magnitudes in the bones overlying the brain. However, local variations in strain magnitudes and patterns can be observed in areas close to the sutures. These results show that, during feeding, sutures may be more important in some regions than others. Sutures should therefore be included in models that require accurate local strain magnitudes and patterns of cranial strain, particularly if models are developed for analysis of specific regions, such as the temporomandibular joint or zygomatic arch. Our results suggest that, for mammalian skulls, cranial sutures might be more important for allowing brain expansion during growth than redistributing biting loads across the cranium in adults.

Multiple pathways to herbivory underpinned deep divergences in ornithischian evolution

The extent to which evolution is deterministic is a key question in biology,^{1,2,3,4,5,6,7,8,9} with intensive debate on how adaptation^{6,10,11,12,13} and constraints^14,15,16 might canalize solutions to ecological challenges.^4,5,6 Alternatively, unique adaptations^1,9,17 and phylogenetic contingency^1,3,18 may render evolution fundamentally unpredictable.³ Information from the fossil record is critical to this debate,^1,2,11 but performance data for extinct taxa are limited.⁷ This knowledge gap is significant, as general morphology may be a poor predictor of biomechanical performance.^17,19,20 High-fiber herbivory originated multiple times within ornithischian dinosaurs,²¹ making them an ideal clade for investigating evolutionary responses to similar ecological pressures.²² However, previous biomechanical modeling studies on ornithischian crania^17,23,24,25 have not compared early-diverging taxa spanning independent acquisitions of herbivory. Here, we perform finite-element analysis on the skull of five early-diverging members of the major ornithischian clades to characterize morphofunctional pathways to herbivory. Results reveal limited functional convergence among ornithischian clades, with each instead achieving comparable performance, in terms of reconstructed patterns and magnitudes of functionally induced stress, through different adaptations of the feeding apparatus. Thyreophorans compensated for plesiomorphic low performance through increased absolute size, heterodontosaurids expanded jaw adductor muscle volume, ornithopods increased jaw system efficiency, and ceratopsians combined these approaches. These distinct solutions to the challenges of herbivory within Ornithischia underpinned the success of this diverse clade. Furthermore, the resolution of multiple solutions to equivalent problems within a single clade through macroevolutionary time demonstrates that phenotypic evolution is not necessarily predictable, instead arising from the interplay of adaptation, innovation, contingency, and constraints.^{1,2,3,7,8,9,18}.

The efficacy of computed tomography scanning versus surface scanning in 3D finite element analysis

Finite element analysis (FEA) is a commonly used application in biomechanical studies of both extant and fossil taxa to assess stress and strain in solid structures such as bone. FEA can be performed on 3D structures that are generated using various methods, including computed tomography (CT) scans and surface scans. While previous palaeobiological studies have used both CT scanned models and surface scanned models, little research has evaluated to what degree FE results may vary when CT scans and surface scans of the same object are compared. Surface scans do not preserve the internal geometries of 3D structures, which are typically preserved in CT scans. Here, we created 3D models from CT scans and surface scans of the same specimens (crania and mandibles of a Nile crocodile, a green sea turtle, and a monitor lizard) and performed FEA under identical loading parameters. It was found that once surface scanned models are solidified, they output stress and strain distributions and model deformations comparable to their CT scanned counterparts, though differing by notable stress and strain magnitudes in some cases, depending on morphology of the specimen and the degree of reconstruction applied. Despite similarities in overall mechanical behaviour, surface scanned models can differ in exterior shape compared to CT scanned models due to inaccuracies that can occur during scanning and reconstruction, resulting in local differences in stress distribution. Solid-fill surface scanned models generally output lower stresses compared to CT scanned models due to their compact interiors, which must be accounted for in studies that use both types of scans.

Enhancing biomedical data validity with standardized segmentation finite element analysis

Finite element analysis is a powerful computational technique for augmenting biomedical research, prosthetics design, and preoperative surgical assessment. However, the validity of biomechanical data obtained from finite element analysis is dependent on the quality of the preceding data processing. Until now, little information was available about the effect of the segmentation process on finite element models and biomechanical data. The current investigation applied 4 segmentation approaches to 129 femur specimens, yielding a total of 516 finite element models. Biomechanical data including average displacement, pressure, stress, and strain were collected from experimental groups based on the different segmentation approaches. The results indicate that only a 5.0% variation in the segmentation process leads to statistically significant differences in all 4 biomechanical measurements. These results suggest that it is crucial for consistent segmentation procedures to be applied to all specimens within a study. This methodological advancement will help to ensure that finite element data will be more accurate and that research conclusions will have greater validity.

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.

Comparative cranial biomechanics in two lizard species: impact of variation in cranial design

Cranial morphology in lepidosaurs is highly disparate and characterised by the frequent loss or reduction of bony elements. In varanids and geckos, the loss of the postorbital bar is associated with changes in skull shape, but the mechanical principles underlying this variation remain poorly understood. Here, we sought to determine how the overall cranial architecture and the presence of the postorbital bar relate to the loading and deformation of the cranial bones during biting in lepidosaurs. Using computer-based simulation techniques, we compared cranial biomechanics in the varanid Varanus niloticus and the teiid Salvator merianae, two large, active foragers. The overall strain magnitude and distribution across the cranium were similar in the two species, despite lower strain gradients in V. niloticus. In S. merianae, the postorbital bar is important for resistance of the cranium to feeding loads. The postorbital ligament, which in varanids partially replaces the postorbital bar, does not affect bone strain. Our results suggest that the reduction of the postorbital bar impaired neither biting performance nor the structural resistance of the cranium to feeding loads in V. niloticus. Differences in bone strain between the two species might reflect demands imposed by feeding and non-feeding functions on cranial shape. Beyond variation in cranial bone strain related to species-specific morphological differences, our results reveal that similar mechanical behaviour is shared by lizards with distinct cranial shapes. Contrary to the situation in mammals, the morphology of the circumorbital region, calvaria and palate appears to be important for withstanding high feeding loads in these lizards.

Summary:In vivo measurements and computer-based simulations of the cranial mechanics of two large lizards indicate that similar mechanical behaviour is shared by lizards with distinct cranial architecture, and show the importance of the postorbital bar in resisting the feeding loads.

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.

Variability, morphometrics, and co-variation of the os lacrimale in Cervidae

In Ruminantia, the lacrimal bone forms a considerable part of the facial skeleton, and the morphology of its facial facet is highly variable when compared to other mammals. In this study, we quantify the species-specific variability in size and shape of the lacrimal facial facet in species of Cervidae (deer) and relate it to systematics and various aspects of their ecology and behavior. We sampled 143 skull specimens from 10 genera; 12 Moschus and 3 Tragulus specimens were used as outgroups. We find that size and shape of the lacrimal facial facet allow differentiating most species analyzed here, except for Mazama gouazoubira and Capreolus capreolus. Size and shape of the lacrimal facial facet vary widely across Cervidae regardless of their systematic relationships, ecology or behavior. Thus, we could not detect a unique signature of adaptational criteria in lacrimal morphology. Our data indicate that the lacrimal facial facet scales allometrically with skull size, in particular, the lacrimojugal length scales positively and the lacrimomaxillar length scales negatively. However, correlation analyses did not reveal any differences in the integration of the lacrimal bone with any specific skull module in any of the species compared. Lastly, we could not ascertain any correlation between the size and position of the preorbital depression with the size and shape of the lacrimal facial facet. We conclude that the lacrimal facial facet is highly flexible and may rapidly adjust to its surrounding bones. Its allometric growth appears to be an example of exaptation: changes in size and shape in the context of the increase of the skull length provide lacrimal contacts, in particular, a lacrimojugal one, which may serve to reduce mechanical loads resulting from increasingly larger antlers in large cervids.

ALLOMETRIC GROWTH IN BRAZILIAN NATURALIZED PIGS

Abstract The allometric growth of male and female naturalized pig breeds (Nilo, Piau, Moura, and Monteiro) was evaluated regarding hip height, snout length, and head length using the allometric equation. Body measurements generally show later growth in relation to the snout length and earlier growth in relation to the head length and hip height. When comparing breeds of pigs, female Moura animals showed later growth of the characteristics linked to muscle (thoracic circumference, body length, and longitudinal diameter) than of head length. Nilo females showed characteristics such as shoulder length, rump length and interschiadic distance earlier than snout length. There were differences between the breed and gender in terms of development times for the different parts of the body, showing that breeding systems and use of these animals have to be specific for each breed.

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 years^1,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 structures^3,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 thought^5-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.

Merging cranial histology and 3D-computational biomechanics: a review of the feeding ecology of a Late Triassic temnospondyl amphibian

Finite Element Analysis (FEA) is a useful method for understanding form and function. However, modelling of fossil taxa invariably involves assumptions as a result of preservation-induced loss of information in the fossil record. To test the validity of predictions from FEA, given such assumptions, these results could be compared to independent lines of evidence for cranial mechanics. In the present study a new concept of using bone microstructure to predict stress distribution in the skull during feeding is put forward and a correlation between bone microstructure and results of computational biomechanics (FEA) is carried out. The bony framework is a product of biological optimisation; bone structure is created to meet local mechanical conditions. To test how well results from FEA correlate to cranial mechanics predicted from bone structure, the well-known temnospondyl Metoposaurus krasiejowensis was used as a model. A crucial issue to Temnospondyli is their feeding mode: did they suction feed or employ direct biting, or both? Metoposaurids have previously been characterised either as active hunters or passive bottom dwellers. In order to test the correlation between results from FEA and bone microstructure, two skulls of Metoposaurus were used, one modelled under FE analyses, while for the second one 17 dermal bone microstructure were analysed. Thus, for the first time, results predicting cranial mechanical behaviour using both methods are merged to understand the feeding strategy of Metoposaurus. Metoposaurus appears to have been an aquatic animal that exhibited a generalist feeding behaviour. This taxon may have used two foraging techniques in hunting; mainly bilateral biting and, to a lesser extent, lateral strikes. However, bone microstructure suggests that lateral biting was more frequent than suggested by Finite Element Analysis (FEA). One of the potential factors that determined its mode of life may have been water levels. During optimum water conditions, metoposaurids may have been more active ambush predators that were capable of lateral strikes of the head. The dry season required a less active mode of life when bilateral biting is particularly efficient. This, combined with their characteristically anteriorly positioned orbits, was optimal for ambush strategy. This ability to use alternative modes of food acquisition, independent of environmental conditions, might hold the key in explaining the very common occurrence of metoposaurids during the Late Triassic.

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.

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.

Validation experiments on finite element models of an ostrich (Struthio camelus) cranium

The first finite element (FE) validation of a complete avian cranium was performed on an extant palaeognath, the ostrich (Struthio camelus). Ex-vivo strains were collected from the cranial bone and rhamphotheca. These experimental strains were then compared to convergence tested, specimen-specific finite element (FE) models. The FE models contained segmented cortical and trabecular bone, sutures and the keratinous rhamphotheca as identified from micro-CT scan data. Each of these individual materials was assigned isotropic material properties either from the literature or from nanoindentation, and the FE models compared to the ex-vivo results. The FE models generally replicate the location of peak strains and reflect the correct mode of deformation in the rostral region. The models are too stiff in regions of experimentally recorded high strain and too elastic in regions of low experimentally recorded low strain. The mode of deformation in the low strain neurocranial region is not replicated by the FE models, and although the models replicate strain orientations to within 10° in some regions, in most regions the correlation is not strong. Cranial sutures, as has previously been found in other taxa, are important for modifying both strain magnitude and strain patterns across the entire skull, but especially between opposing the sutural junctions. Experimentally, we find that the strains on the surface of the rhamphotheca are much lower than those found on nearby bone. The FE models produce much higher principal strains despite similar strain ratios across the entirety of the rhamphotheca. This study emphasises the importance of attempting to validate FE models, modelling sutures and rhamphothecae in birds, and shows that whilst location of peak strain and patterns of deformation can be modelled, replicating experimental data in digital models of avian crania remains problematic.

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.

The impact of simplifications on the performance of a finite element model of a Macaca fascicularis cranium

In recent years finite element analysis (FEA) has emerged as a useful tool for the analysis of skeletal form-function relationships. While this approach has obvious appeal for the study of fossil specimens, such material is often fragmentary with disrupted internal architecture and can contain matrix that leads to errors in accurate segmentation. Here we examine the effects of varying the detail of segmentation and material properties of teeth on the performance of a finite element model of a Macaca fascicularis cranium within a comparative functional framework. Cranial deformations were compared using strain maps to assess differences in strain contours and Procrustes size and shape analyses, from geometric morphometrics, were employed to compare large scale deformations. We show that a macaque model subjected to biting can be made solid, and teeth altered in material properties, with minimal impact on large scale modes of deformation. The models clustered tightly by bite point rather than by modeling simplification approach, and fell out as being distinct from another species. However localized fluctuations in predicted strain magnitudes were recorded with different modeling approaches, particularly over the alveolar region. This study indicates that, while any model simplification should be undertaken with care and attention to its effects, future applications of FEA to fossils with unknown internal architecture may produce reliable results with regard to general modes of deformation, even when detail of internal bone architecture cannot be reliably modeled.

The relationship between cranial structure, biomechanical performance and ecological diversity in varanoid lizards

Skull structure is intimately associated with feeding ability in vertebrates, both in terms of specific performance measures and general ecological characteristics. This study quantitatively assessed variation in the shape of the cranium and mandible in varanoid lizards, and its relationship to structural performance (von Mises strain) and interspecific differences in feeding ecology. Geometric morphometric and linear morphometric analyses were used to evaluate morphological differences, and finite element analysis was used to quantify variation in structural performance (strain during simulated biting, shaking and pulling). This data was then integrated with ecological classes compiled from relevant scientific literature on each species in order to establish structure-function relationships. Finite element modelling results showed that variation in cranial morphology resulted in large differences in the magnitudes and locations of strain in biting, shaking and pulling load cases. Gracile species such as Varanus salvadorii displayed high strain levels during shaking, especially in the areas between the orbits. All models exhibit less strain during pull back loading compared to shake loading, even though a larger force was applied (pull =30N, shake = 20N). Relationships were identified between the morphology, performance, and ecology. Species that did not feed on hard prey clustered in the gracile region of cranial morphospace and exhibited significantly higher levels of strain during biting (P = 0.0106). Species that fed on large prey clustered in the elongate area of mandible morphospace. This relationship differs from those that have been identified in other taxonomic groups such as crocodiles and mammals. This difference may be due to a combination of the open 'space-frame' structure of the varanoid lizard skull, and the 'pull back' behaviour that some species use for processing large prey.

The sensitivity of biological finite element models to the resolution of surface geometry: a case study of crocodilian crania

The reliability of finite element analysis (FEA) in biomechanical investigations depends upon understanding the influence of model assumptions. In producing finite element models, surface mesh resolution is influenced by the resolution of input geometry, and influences the resolution of the ensuing solid mesh used for numerical analysis. Despite a large number of studies incorporating sensitivity studies of the effects of solid mesh resolution there has not yet been any investigation into the effect of surface mesh resolution upon results in a comparative context. Here we use a dataset of crocodile crania to examine the effects of surface resolution on FEA results in a comparative context. Seven high-resolution surface meshes were each down-sampled to varying degrees while keeping the resulting number of solid elements constant. These models were then subjected to bite and shake load cases using finite element analysis. The results show that incremental decreases in surface resolution can result in fluctuations in strain magnitudes, but that it is possible to obtain stable results using lower resolution surface in a comparative FEA study. As surface mesh resolution links input geometry with the resulting solid mesh, the implication of these results is that low resolution input geometry and solid meshes may provide valid results in a comparative context.

Descriptive anatomy and three-dimensional reconstruction of the skull of the early tetrapod Acanthostega gunnari Jarvik, 1952

The early tetrapod Acanthostega gunnari is an iconic fossil taxon exhibiting skeletal morphology reflecting the transition of vertebrates from water onto land. Computed tomography data of two Acanthostega skulls was segmented using visualization software to digitally separate bone from matrix and individual bones of the skull from each other. A revised description of cranial and lower jaw anatomy in this taxon based on CT data includes new details of sutural morphology, the previously undescribed quadrate and articular bones, and the mandibular symphysis. Sutural morphology is used to infer loading regime in the skull during feeding, and suggests Acanthostega used its anterior jaws to initially seize prey while smaller posterior teeth were used to restrain struggling prey during ingestion. Novel methods were used to repair and retrodeform the skull, resulting in a three-dimensional digital reconstruction that features a longer postorbital region and more strongly hooked anterior lower jaw than previous attempts while supporting the presence of a midline gap between the nasals and median rostrals.

Beware the black box: investigating the sensitivity of FEA simulations to modelling factors in comparative biomechanics

Finite element analysis (FEA) is a computational technique of growing popularity in the field of comparative biomechanics, and is an easily accessible platform for form-function analyses of biological structures. However, its rapid evolution in recent years from a novel approach to common practice demands some scrutiny in regards to the validity of results and the appropriateness of assumptions inherent in setting up simulations. Both validation and sensitivity analyses remain unexplored in many comparative analyses, and assumptions considered to be ‘reasonable’ are often assumed to have little influence on the results and their interpretation.

Here we report an extensive sensitivity analysis where high resolution finite element (FE) models of mandibles from seven species of crocodile were analysed under loads typical for comparative analysis: biting, shaking, and twisting. Simulations explored the effect on both the absolute response and the interspecies pattern of results to variations in commonly used input parameters. Our sensitivity analysis focuses on assumptions relating to the selection of material properties (heterogeneous or homogeneous), scaling (standardising volume, surface area, or length), tooth position (front, mid, or back tooth engagement), and linear load case (type of loading for each feeding type).

Our findings show that in a comparative context, FE models are far less sensitive to the selection of material property values and scaling to either volume or surface area than they are to those assumptions relating to the functional aspects of the simulation, such as tooth position and linear load case. Results show a complex interaction between simulation assumptions, depending on the combination of assumptions and the overall shape of each specimen. Keeping assumptions consistent between models in an analysis does not ensure that results can be generalised beyond the specific set of assumptions used. Logically, different comparative datasets would also be sensitive to identical simulation assumptions; hence, modelling assumptions should undergo rigorous selection. The accuracy of input data is paramount, and simulations should focus on taking biological context into account. Ideally, validation of simulations should be addressed; however, where validation is impossible or unfeasible, sensitivity analyses should be performed to identify which assumptions have the greatest influence upon the results.

Cranial sutures work collectively to distribute strain throughout the reptile skull

The skull is composed of many bones that come together at sutures. These sutures are important sites of growth, and as growth ceases some become fused while others remain patent. Their mechanical behaviour and how they interact with changing form and loadings to ensure balanced craniofacial development is still poorly understood. Early suture fusion often leads to disfiguring syndromes, thus is it imperative that we understand the function of sutures more clearly. By applying advanced engineering modelling techniques, we reveal for the first time that patent sutures generate a more widely distributed, high level of strain throughout the reptile skull. Without patent sutures, large regions of the skull are only subjected to infrequent low-level strains that could weaken the bone and result in abnormal development. Sutures are therefore not only sites of bone growth, but could also be essential for the modulation of strains necessary for normal growth and development in reptiles.

Blunt liver injury with intact ribs under impacts on the abdomen: a biomechanical investigation

Abdominal trauma accounts for nearly 20% of all severe traffic injuries and can often result from intentional physical violence, from which blunt liver injury is regarded as the most common result and is associated with a high mortality rate. Liver injury may be caused by a direct impact with a certain velocity and energy on the abdomen, which may result in a lacerated liver by penetration of fractured ribs. However, liver ruptures without rib cage fractures were found in autopsies in a series of cases. All the victims sustained punches on the abdomen by fist. Many studies have been dedicated to determining the mechanism underlying hepatic injury following abdominal trauma, but most have been empirical. The actual process and biomechanism of liver injury induced by blunt impact on the abdomen, especially with intact ribs remained, are still inexhaustive. In order to investigate this, finite element methods and numerical simulation technology were used. A finite element human torso model was developed from high resolution CT data. The model consists of geometrically-detailed liver and rib cage models and simplified models of soft tissues, thoracic and abdominal organs. Then, the torso model was used in simulations in which the right hypochondrium was punched by a fist from the frontal, lateral, and rear directions, and in each direction with several impact velocities. Overall, the results showed that liver rupture was primarily caused by a direct strike of the ribs induced by blunt impact to the abdomen. Among three impact directions, a lateral impact was most likely to cause liver injury with a minimum punch speed of 5 m/s (the momentum was about 2.447 kg.m/s). Liver injuries could occur in isolation and were not accompanied by rib fractures due to different material characteristics and injury tolerance.

A sensitivity analysis to the role of the fronto-parietal suture in Lacerta bilineata: a preliminary finite element study

Cranial sutures are sites of bone growth and development but micromovements at these sites may distribute the load across the skull more evenly. Computational studies have incorporated sutures into finite element (FE) models to assess various hypotheses related to their function. However, less attention has been paid to the sensitivity of the FE results to the shape, size, and stiffness of the modeled sutures. Here, we assessed the sensitivity of the strain predictions to the aforementioned parameters in several models of fronto-parietal (FP) suture in Lacerta bilineata. For the purpose of this study, simplifications were made in relation to modeling the bone properties and the skull loading. Results highlighted that modeling the FP as either an interdigitated suture or a simplified butt suture, did not reduce the strain distribution in the FP region. Sensitivity tests showed that similar patterns of strain distribution can be obtained regardless of the size of the suture, or assigned stiffness, yet the exact magnitudes of strains are highly sensitive to these parameters. This study raises the question whether the morphogenesis of epidermic scales in the FP region in the Lacertidae is related to high strain fields in this region, because of micromovement in the FP suture.