Sutures form an integral part of the functioning skull, but their role has long been debated among vertebrate morphologists and palaeontologists. Furthermore, the relationship between typical skull sutures, and those involved in cranial kinesis, is poorly understood. In a series of computational modelling studies, complex loading conditions obtained through multibody dynamics analysis were imposed on a finite element model of the skull ofUromastyx hardwickii, an akinetic herbivorous lizard. A finite element analysis (FEA) of a skull with no sutures revealed higher patterns of strain in regions where cranial sutures are located in the skull. From these findings, FEAs were performed on skulls with sutures (individual and groups of sutures) to investigate their role and function more thoroughly. Our results showed that individual sutures relieved strain locally, but only at the expense of elevated strain in other regions of the skull. These findings provide an insight into the behaviour of sutures and show how they are adapted to work together to distribute strain around the skull. Premature fusion of one suture could therefore lead to increased abnormal loading on other regions of the skull causing irregular bone growth and deformities. This detailed investigation also revealed that the frontal–parietal suture of theUromastyxskull played a substantial role in relieving strain compared with the other sutures. This raises questions about the original role of mesokinesis in squamate evolution.

Additive manufacturing has been used to develop a variety of scaffold designs for clinical and industrial applications. Mechanical properties (i.e., compression, tension, bending, and torsion response) of these scaffolds are significantly important for load-bearing orthopaedic implants. In this study, we designed and additively manufactured porous metallic biomaterials based on two different types of triply periodic minimal surface structures (i.e., gyroid and diamond) that mimic the mechanical properties of bone, such as porosity, stiffness, and strength. Physical and mechanical properties, including compressive, tensile, bending, and torsional stiffness and strength of the developed scaffolds, were then characterised experimentally and numerically using finite element method. Sheet thickness was constant at 300 μm, and the unit cell size was varied to generate different pore sizes and porosities. Gyroid scaffolds had a pore size in the range of 600–1200 μm and a porosity in the range of 54–72%, respectively. Corresponding values for the diamond were 900–1500 μm and 56–70%. Both structure types were validated experimentally, and a wide range of mechanical properties (including stiffness and yield strength) were predicted using the finite element method. The stiffness and strength of both structures are comparable to that of cortical bone, hence reducing the risks of scaffold failure. The results demonstrate that the developed scaffolds mimic the physical and mechanical properties of cortical bone and can be suitable for bone replacement and orthopaedic implants. However, an optimal design should be chosen based on specific performance requirements.

This study aimed to clarify the influence of various speeding strategies (i.e. adjustments of cadence and stride length) on external joint moments. This study investigated the gait of 52 healthy subjects who performed self-selected normal and fast speed walking trials in a motion analysis laboratory. Subjects were classified into three separate groups based on how they increased their speed from normal to fast walking: (i) subjects who increased their cadence, (ii) subjects who increased their stride length and (iii) subjects who simultaneously increased both stride length and cadence. Joint moments were calculated using inverse dynamics and then compared between normal and fast speed trials within and between three groups using spatial parameter mapping. Individuals who increased cadence, but not stride length, to walk faster did not experience a significant increase in the lower limb joint moments. Conversely, subjects who increased their stride length or both stride length and cadence, experienced a significant increase in all joint moments. Additionally, our findings revealed that increasing the stride length had a higher impact on joint moments in the sagittal plane than those in the frontal plane. However, both sagittal and frontal plane moments were still more responsive to the gait speed change than transverse plane moments. This study suggests that the role of speed in altering the joint moment patterns depends on the individual's speed-regulating strategy, i.e. an increase in cadence or stride length. Since the confounding effect of walking speed is a major consideration in human gait research, future studies may investigate whether stride length is the confounding variable of interest.

Evaluating stress and strain fields in anatomical structures is a way to test hypotheses that relate specific features of facial and skeletal morphology to mechanical loading. Engineering techniques such as finite element analysis are now commonly used to calculate stress and strain fields, but if we are to fully accept these methods we must be confident that the applied loading regimens are reasonable. Multibody dynamics analysis (MDA) is a relatively new three dimensional computer modeling technique that can be used to apply varying muscle forces to predict joint and bite forces during static and dynamic motions. The method ensures that equilibrium of the structure is maintained at all times, even for complex statically indeterminate problems, eliminating nonphysiological constraint conditions often seen with other approaches. This study describes the novel use of MDA to investigate the influence of different muscle representations on a macaque skull model (Macaca fascicularis), where muscle groups were represented by either a single, multiple, or wrapped muscle fibers. The impact of varying muscle representation on stress fields was assessed through additional finite element simulations. The MDA models highlighted that muscle forces varied with gape and that forces within individual muscle groups also varied; for example, the anterior strands of the superficial masseter were loaded to a greater extent than the posterior strands. The direction of the muscle force was altered when temporalis muscle wrapping was modeled, and was coupled with compressive contact forces applied to the frontal, parietal and temporal bones of the cranium during biting. Anat Rec, 291:491–501, 2008.

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.