Crucial to the interpretation of the results of any finite element analysis of a skeletal system is a test of the validity of the results and an assessment of the sensitivity of the model parameters. We have therefore developed finite element models of two crania of Macaca fascicularis and investigated their sensitivity to variations in bone material properties, the zygomatico-temporal suture and the loading regimen applied to the zygomatic arch. Maximum principal strains were validated against data derived from ex vivo strain gauge experiments using non-physiological loads applied to the macaque zygomatic arch. Elastic properties of the zygomatic arch bone and the zygomaticotemporal suture obtained by nanoindentation resulted in a high degree of congruence between experimental and simulated strains. The findings also indicated that the presence of a zygomatico-temporal suture in the model produced strains more similar to experimental values than a completely separated or fused arch. Strains were distinctly higher when the load was applied through the modelled superficial masseter compared with loading an array of nodes on the arch. This study demonstrates the importance of the accurate selection of the material properties involved in predicting strains in a finite element model. Furthermore, our findings strongly highlight the influence of the presence of craniofacial sutures on strains experienced in the face. This has implications when investigating craniofacial growth and masticatory function but should generally be taken into account in functional analyses of the craniofacial system of both extant and extinct species.

Research on the evolution and adaptive significance of primate craniofacial morphologies has focused on adult, fully developed individuals. Here, we investigate the possible relationship between the local stress environment arising from masticatory loadings and the emergence of the supraorbital torus in the developing face of the crab-eating macaque Macaca fascicularis. By using finite element analysis (FEA), we are able to evaluate the hypothesis that strain energy density (SED) magnitudes are high in subadult individuals with resulting bone growth in the supraorbital torus. We developed three micro-CT-based FEA models of M. fascicularis skulls ranging in dental age from deciduous to permanent dentitions and validated them against published experimental data. Applied masticatory muscle forces were estimated from physiological cross-sectional areas of macaque cadaveric specimens. The models were sequentially constrained at each working side tooth to simulate the variation of the bite point applied during masticatory function. Custom FEA software was used to solve the voxel-based models and SED and principal strains were computed. A physiological superposition SED map throughout the face was created by allocating to each element the maximum SED value from each of the load cases. SED values were found to be low in the supraorbital torus region throughout ontogeny, while they were consistently high in the zygomatic arch and infraorbital region. Thus, if the supraorbital torus arises to resist masticatory loads, it is either already adapted in each of our subadult models so that we do not observe high SED or a lower site-specific bone deposition threshold must apply.

The use of finite element analysis (FEA) to investigate the biomechanics of anatomical systems critically relies on the specification of physiologically representative boundary conditions. The biomechanics of the pelvis has been the specific focus of a number of FEA studies previously, but it is also a key aspect in other investigations of, for example, the hip joint or new design of hip prostheses. In those studies, the pelvis has been modelled in a number of ways with a variety of boundary conditions, ranging from a model of the whole pelvic girdle including soft tissue attachments to a model of an isolated hemi-pelvis. The current study constructed a series of FEA models of the same human pelvis to investigate the sensitivity of the predicted stress distributions to the type of boundary conditions applied, in particular to represent the sacro-iliac joint and pubic symphysis. Varying the method of modelling the sacro-iliac joint did not produce significant variations in the stress distribution, however changes to the modelling of the pubic symphysis were observed to have a greater effect on the results. Over-constraint of the symphysis prevented the bending of the pelvis about the greater sciatic notch, and underestimated high stresses within the ilium. However, permitting medio-lateral translation to mimic widening of the pelvis addressed this problem. These findings underline the importance of applying the appropriate boundary conditions to FEA models, and provide guidance on suitable methods of constraining the pelvis when, for example, scan data has not captured the full pelvic girdle. The results also suggest a valid method for performing hemi-pelvic modelling of cadaveric or archaeological remains which are either damaged or incomplete.