Strain rate dependent properties of human craniovertebral ligaments

Simulative investigation of the required level of geometrical individualization of the lumbar spines to predict fractures

Injury mechanisms of the lumbar spine under dynamic loading are dependent on spine curvature and anatomical variation. Impact simulation with finite element (FE) models can assist the reconstruction and prediction of injuries. The objective of this study was to determine which level of individualization of a baseline FE lumbar spine model is necessary to replicate experimental responses and fracture locations in a dynamic experiment.Experimental X-rays from 26 dynamic drop tower tests were used to create three configurations of a lumbar spine model (T12 to L5): baseline, with aligned vertebrae (positioned), and with aligned and morphed vertebrae (morphed). Each model was simulated with the corresponding loading and boundary conditions from dynamic lumbar spine experiments. Force, moment, and kinematic responses were compared to the experimental data. Cosine similarity was computed to assess how well simulation responses match the experimental data. The pressure distribution within the vertebrae was used to compare fracture risk and fracture location between the different models.The positioned models replicated the injured spinal level and the fracture patterns quite well, though the morphed models provided slightly more accuracy. However, for impact reconstruction or injury prediction, the authors recommend pure positioning for whole-body models, as the gain in accuracy was relatively small, while the morphing modifications of the model require considerably higher efforts. These results improve the understanding of the application of human body models to investigate lumbar injury mechanisms with FE models.

Subject-Specific Geometry of FE Lumbar Spine Models for the Replication of Fracture Locations Using Dynamic Drop Tests

For traumatic lumbar spine injuries, the mechanisms and influence of anthropometrical variation are not yet fully understood under dynamic loading. Our objective was to evaluate whether geometrically subject-specific explicit finite element (FE) lumbar spine models based on state-of-the-art clinical CT data combined with general material properties from the literature could replicate the experimental responses and the fracture locations via a dynamic drop tower-test setup. The experimental CT datasets from a dynamic drop tower-test setup were used to create anatomical details of four lumbar spine models (T12 to L5). The soft tissues from THUMS v4.1 were integrated by morphing. Each model was simulated with the corresponding loading and boundary conditions from the dynamic lumbar spine tests that produced differing injuries and injury locations. The simulations resulted in force, moment, and kinematic responses that effectively matched the experimental data. The pressure distribution within the models was used to compare the fracture occurrence and location. The spinal levels that sustained vertebral body fracture in the experiment showed higher simulation pressure values in the anterior elements than those in the levels that did not fracture in the reference experiments. Similarly, the spinal levels that sustained posterior element fracture in the experiments showed higher simulation pressure values in the vertebral posterior structures compared to those in the levels that did not sustain fracture. Our study showed that the incorporation of the spinal geometry and orientation could be used to replicate the fracture type and location under dynamic loading. Our results provided an understanding of the lumbar injury mechanisms and knowledge on the load thresholds that could be used for injury prediction with explicit FE lumbar spine models.

Applied use of biomechanical measurements from human tissues for the development of medical skills trainers: a scoping review

OBJECTIVE

The objective of this review was to identify quantitative biomechanical measurements of human tissues, the methods for obtaining these measurements, and the primary motivations for conducting biomechanical research.

INTRODUCTION

Medical skills trainers are a safe and useful tool for clinicians to use when learning or practicing medical procedures. The haptic fidelity of these devices is often poor, which may be because the synthetic materials chosen for these devices do not have the same mechanical properties as human tissues. This review investigates a heterogeneous body of literature to identify which biomechanical properties are available for human tissues, the methods for obtaining these values, and the primary motivations behind conducting biomechanical tests.

INCLUSION CRITERIA

Studies containing quantitative measurements of the biomechanical properties of human tissues were included. Studies that primarily focused on dynamic and fluid mechanical properties were excluded. Additionally, studies only containing animal, in silico , or synthetic materials were excluded from this review.

METHODS

This scoping review followed the JBI methodology for scoping reviews and the Preferred Reporting Items for Systematic Reviews and Meta-Analyses extension for Scoping Reviews (PRISMA-ScR). Sources of evidence were extracted from CINAHL (EBSCO), IEEE Xplore, MEDLINE (PubMed), Scopus, and engineering conference proceedings. The search was limited to the English language. Two independent reviewers screened titles and abstracts as well as full-text reviews. Any conflicts that arose during screening and full-text review were mediated by a third reviewer. Data extraction was conducted by 2 independent reviewers and discrepancies were mediated through discussion. The results are presented in tabular, figure, and narrative formats.

RESULTS

Data were extracted from a total of 186 full-text publications. All of the studies, except for 1, were experimental. Included studies came from 33 countries, with the majority coming from the United States. Ex vivo methods were the predominant approach for extracting human tissue samples, and the most commonly studied tissue type was musculoskeletal. In this study, nearly 200 unique biomechanical values were reported, and the most commonly reported value was Young's (elastic) modulus. The most common type of mechanical test performed was tensile testing, and the most common reason for testing human tissues was to characterize biomechanical properties. Although the number of published studies on biomechanical properties of human tissues has increased over the past 20 years, there are many gaps in the literature. Of the 186 included studies, only 7 used human tissues for the design or validation of medical skills training devices. Furthermore, in studies where biomechanical values for human tissues have been obtained, a lack of standardization in engineering assumptions, methodologies, and tissue preparation may implicate the usefulness of these values.

CONCLUSIONS

This review is the first of its kind to give a broad overview of the biomechanics of human tissues in the published literature. With respect to high-fidelity haptics, there is a large gap in the published literature. Even in instances where biomechanical values are available, comparing or using these values is difficult. This is likely due to the lack of standardization in engineering assumptions, testing methodology, and reporting of the results. It is recommended that journals and experts in engineering fields conduct further research to investigate the feasibility of implementing reporting standards.

REVIEW REGISTRATION

Open Science Framework https://osf.io/fgb34.

Quantitatively biomechanical response analysis of posterior musculature reconstruction in cervical single-door laminoplasty

BACKGROUND AND OBJECTIVE

The current trend of laminoplasty is developing toward the goal of muscle preservation and minimum tissue damage. Given this, muscle-preserving techniques in cervical single-door laminoplasty have been modified with protecting the spinous processes at the sites of C2 and/or C7 muscle attachment and reconstruct the posterior musculature in recent years. To date, no study has reported the effect of preserving the posterior musculature during the reconstruction. The purpose of this study is to quantitatively evaluate the biomechanical effect of multiple modified single-door laminoplasty procedures for restoring stability and reducing response level on the cervical spine.

METHODS

Different cervical laminoplasty models were established for evaluating kinematics and response simulations based on a detailed finite element (FE) head-neck active model (HNAM), including ① C3 - C7 laminoplasty (LP_C37), ② C3 - C6 laminoplasty with C7 spinous process preservation (LP_C36), ③ C3 laminectomy hybrid decompression with C4 - C6 laminoplasty (LT_C3 + LP_C46) and ④ C3 - C7 laminoplasty with unilateral musculature preservation (LP_C37 + UMP). The laminoplasty model was validated by the global range of motion (ROM) and percentage changes relative to the intact state. The C2 - T1 ROM, axial muscle tensile force, and stress/strain levels of functional spinal units were compared among the different laminoplasty groups. The obtained effects were further analysed by comparison with a review of clinical data on cervical laminoplasty scenarios.

RESULTS

Analysis of the locations of concentration of muscle load showed that the C2 muscle attachment sustained more tensile loading than the C7 muscle attachment, primarily in flexion-extension (FE) and in lateral bending (LB) and axial rotation (AR), respectively. Simulated results further quantified that LP_C36 primarily produced 10% decreases in LB and AR modes relative to LP_C37. Compared with LP_C36, LT_C3 + LP_C46 resulted in approximately 30% decreases in FE motion; LP C37 + UMP also showed a similar trend. Additionally, when compared to LP_C37, LT_C3 + LP_C46 and LP C37 + UMP reduced the peak stress level at the intervertebral disc by at most 2-fold as well as the peak strain level of the facet joint capsule by 2-3-fold. All these findings were well correlated with the result of clinical studies comparing modified laminoplasty and classic laminoplasty.

CONCLUSIONS

Modified muscle-preserving laminoplasty is superior to classic laminoplasty due to the biomechanical effect of the posterior musculature reconstruction, with a retained postoperative ROM and loading response levels of the functional spinal units. More motion-sparing is beneficial for increasing cervical stability, which probably accelerates the recovery of postoperative neck movement and reduces the risk of the complication for eventual kyphosis and axial pain. Surgeons are encouraged to make every effort to preserve the attachment of the C2 whenever feasible in laminoplasty.

Head supported mass, moment of inertia, neck loads and stability: A simulation study

Occupations or activities where donning head-supported mass (HSM) is commonplace put operators at an elevated risk of chronic neck pain. Yet, there is no consensus about what features of HSM influence the relative contributions to neck loads. Therefore, we tested four hypotheses that could increase neck loads: (i) HSM increases gravitational moments; (ii) more muscle activation is required to stabilize the head with HSM; (iii) the position of the HSM centre of mass (COM) induces gravitational moments; and (iv) the added moment of inertia (MOI) from HSM increases neck loads during head repositioning tasks. We performed a sensitivity analysis on the C5-C6 compression evaluated from a 24-degree freedom cervical spine model in OpenSim for static and dynamic movement trials. For static trials, we varied the magnitude of HSM, the position of its COM, and developed a novel stability constraint for static optimization. In dynamic trials, we varied HSM and the three principle MOIs. HSM magnitude and compression were linearly related to one another for both static and dynamic trials, with amplification factors varying between 1.9 and 3.9. Similar relationships were found for the COM position, although the relationship between C5-C6 peak compression and MOI in dynamic trials was generally nonlinear. This sensitivity analysis uncovered evidence in favour of hypotheses (i), (ii) and (iii). However, the model's prediction of C5-C6 compression was not overly sensitive to the magnitude of MOI. Therefore, the HSM mass properties may be more influential on neck compression than MOI properties, even during dynamic tasks.

Quantitative cervical spine injury responses in whiplash loading with a numerical method of natural neural reflex consideration

BACKGROUND AND OBJECTIVE

Neural reflex is hypothesized as a regulating step in spine stabilizing system. However, neural reflex control is still in its infancy to consider in the previous finite element analysis of head-neck system for various applications. The purpose of this study is to investigate the influences of neural reflex control on neck biomechanical responses, then provide a new way to achieve an accurate biomechanical analysis for head-neck system with a finite element model.

METHODS

A new FE head-neck model with detailed active muscles and spinal cord modeling was established and globally validated at multi-levels. Then, it was coupled with our previously developed neuromuscular head-neck model to analyze the effects of vestibular and proprioceptive reflexes on biomechanical responses of head-neck system in a typical spinal injury loading condition (whiplash). The obtained effects were further analyzed by comparing a review of epidemiologic data on cervical spine injury situations.

RESULT

The results showed that the active model (AM) with neural reflex control obviously presented both rational head-neck kinematics and tissue injury risk referring to the previous experimental and epidemiologic studies, when compared with the passive model (PM) without it. Tissue load concentration locations as well as stress/strain levels were both changed due to the muscle activation forces caused by neural reflex control during the whole loading process. For the bony structures, the AM showed a peak stress level accounting for only about 25% of the PM. For the discs, the stress concentrated location was transferred from C2-C6 in the PM to C4-C6 in the AM. For the spinal cord, the strain concentrated locations were transferred from C1 segment to around C4 segment when the effects of neural reflex control were implemented, while the gray matter and white matter peak strains were reduced to 1/3 and 1/2 of the PM, respectively. All these were well correlated with epidemiological studies on clinical cervical spine injuries.

CONCLUSION

In summary, the present work demonstrated necessity of considering neural reflex in FE analysis of a head-neck system as well as our model biofidelity. Overall results also verified the previous hypothesis and further quantitatively indicated that the muscle activation caused by neural reflex is providing a protection for the neck in impact loading by decreasing the strain level and changing the possible injury to lower spinal cord level to reduce injury severity.

An Electromyographically Driven Cervical Spine Model in OpenSim

Relatively few biomechanical models exist aimed at quantifying the mechanical risk factors associated with neck pain. In addition, there is a need to validate spinal-rhythm techniques for inverse dynamics spine models. Therefore, the present investigation was 3-fold: (1) the development of a cervical spine model in OpenSim, (2) a test of a novel spinal-rhythm technique based on minimizing the potential energy in the passive tissues, and (3) comparison of an electromyographically driven approach to estimating compression and shear to other cervical spine models. The authors developed ligament force-deflection and intervertebral joint moment-angle curves from published data. The 218 Hill-type muscle elements, representing 58 muscles, were included and their passive forces validated against in vivo data. Our novel spinal-rhythm technique, based on minimizing the potential energy in the passive tissues, disproportionately assigned motion to the upper cervical spine that was not physiological. Finally, using kinematics and electromyography collected from 8 healthy male volunteers, the authors calculated the compression at C7-T1 as a function of the head-trunk Euler angles. Differences from other models varied from 25.5 to 368.1 N. These differences in forces may result in differences in model geometry, passive components, number of degrees of freedom, or objective functions.

A Hyper-Viscoelastic Continuum-Level Finite Element Model of the Spinal Cord Assessed for Transverse Indentation and Impact Loading

Finite Element (FE) modelling of spinal cord response to impact can provide unique insights into the neural tissue response and injury risk potential. Yet, contemporary human body models (HBMs) used to examine injury risk and prevention across a wide range of impact scenarios often lack detailed integration of the spinal cord and surrounding tissues. The integration of a spinal cord in contemporary HBMs has been limited by the need for a continuum-level model owing to the relatively large element size required to be compatible with HBM, and the requirement for model development based on published material properties and validation using relevant non-linear material data. The goals of this study were to develop and assess non-linear material model parameters for the spinal cord parenchyma and pia mater, and incorporate these models into a continuum-level model of the spinal cord with a mesh size conducive to integration in HBM. First, hyper-viscoelastic material properties based on tissue-level mechanical test data for the spinal cord and hyperelastic material properties for the pia mater were determined. Secondly, the constitutive models were integrated in a spinal cord segment FE model validated against independent experimental data representing transverse compression of the spinal cord-pia mater complex (SCP) under quasi-static indentation and dynamic impact loading. The constitutive model parameters were fit to a quasi-linear viscoelastic model with an Ogden hyperelastic function, and then verified using single element test cases corresponding to the experimental strain rates for the spinal cord (0.32–77.22 s⁻¹) and pia mater (0.05 s⁻¹). Validation of the spinal cord model was then performed by re-creating, in an explicit FE code, two independent ex-vivo experimental setups: 1) transverse indentation of a porcine spinal cord-pia mater complex and 2) dynamic transverse impact of a bovine SCP. The indentation model accurately matched the experimental results up to 60% compression of the SCP, while the impact model predicted the loading phase and the maximum deformation (within 7%) of the SCP experimental data. This study quantified the important biomechanical contribution of the pia mater tissue during spinal cord deformation. The validated material models established in this study can be implemented in computational HBM.

Predicting Cervical Spine Compression and Shear in Helicopter Helmeted Conditions Using Artificial Neural Networks

OCCUPATIONAL APPLICATIONSMilitary helicopter pilots around the globe are at high risk of neck pain related to their use of helmet-mounted night vision goggles. Unfortunately, it is difficult to design alternative helmet configurations that reduce the biomechanical exposures on the cervical spine during flight because the time and resource costs associated with assessing these exposures in vivo are prohibitive. Instead, we developed artificial neural networks (ANNs) to predict cervical spine compression and shear given head-trunk kinematics and joint moments in the lower neck, data readily available from digital human models. The ANNs detected differences in cervical spine compression and anteroposterior shear between helmet configuration conditions during flight-relevant head movement, consistent with results from a detailed model based on in vivo electromyographic data. These ANNs may be useful in helping to prevent neck pain related to military helicopter flight by facilitating virtual biomechanical assessment of helmet configurations upstream in the design process.

The ratio of the posterior atlanto-occipital interval (PAOI): a novel radiographic ratio method evaluating the risk of cervical spondylotic myelopathy-a case-control study

Background

Our study aims to introduce a dynamic interval ratio method calculated using cervical hyperextension-flexion X-ray films. Secondarily, we aim to evaluate the relationship between the posterior atlanto-occipital interval ratio and cervical spondylotic myelopathy and explain the rationale.

Methods

We reviewed 83 cases with visible cervical dynamic X-ray films in our hospital from February 2015 to December 2018. Cases were divided into 2 groups according to their diagnosis (with or without spondylotic myelopathy). Radiographic measurements included the shortest distance between the posterior arch of the atlas and the occipital bone and cervical range of motion, and demographic data such as gender, age, and body mass index were also extracted. The posterior atlanto-occipital interval ratio (distance at hyperextension position/distance at hyperflexion position) was determined using logistic regression analysis models between the 2 groups.

Results

We included 40 cases in the disease group and 43 cases in the control group. The mean posterior atlanto-occipital interval ratio was 0.65±0.30 (mean ± standard deviation) in the disease group and 0.30±0.28 in the control group, with a significant difference (P<0.01). There was no correlation between the posterior atlanto-occipital interval ratio and gender or body mass index. However, the interval ratio had strong correlations with age, cervical spondylotic myelopathy, and Japanese Orthopaedic Association scores. Age, posterior atlanto-occipital interval ratio, and interval distance at hyperextension in the disease group were higher than those of the control group. Contrastingly, range of motion, Japanese Orthopaedic Association scores, and interval distance at the disease group's hyperflexion position were lower than in the control group. In all cases, the risk of cervical spondylotic myelopathy in the T2 group (cases with middle posterior atlanto-occipital interval ratio, according to the tertiles) was 6 times more than the T1 group (cases with lower ratio), and the T3 group (cases with higher ratio) had a 26.4 times greater risk than the T1 group.

Conclusions

Our results suggest that the posterior atlanto-occipital interval ratio is a simple and meaningful parameter that could provide prognostic value for the risk of cervical spondylotic myelopathy through the imaging examinations of the selected cases. Higher posterior atlanto-occipital interval ratios indicate a greater risk for cervical spondylotic myelopathy and cervical musculoskeletal dysfunction. A higher posterior atlanto-occipital interval ratio may manifest undetected posterior atlanto-occipital stiffness, which needs more pathological evidence in future studies.

Effects of different severities of disc degeneration on the range of motion of cervical spine

Aims and Objectives:

The human spine degenerates with age. Intervertebral disc degeneration occurs in the cervical spine. The objective of this study is to determine the effects of degenerative disc diseases on the range of motion (ROM) of the human cervical spinal column using a validated finite-element model.

Materials and Methods:

The validated intact and healthy C2–T1 finite-element model simulated the cortical shell, cancellous core, posterior elements of the vertebrae, and spinal ligaments (longitudinal, capsular, spinous and ligamentum flava, and nucleus and annulus of the discs). Three different stages of the disc disease, that is, mild, moderate, and severe, were simulated at the C5–C6, C6–C7, and C5–C6–C7 discs, respectively, and they were termed as upper single level, lower single level, and bi-level (BL) models, respectively. The material properties and geometry of the disc(s) were altered to simulate the different stages of degeneration. The external mechanical loading was applied in the sagittal mode, via flexion–extension motions and the magnitude was 2.0 Nm for each mode. They were applied to each of the healthy and disc degeneration models, and for each of the three severities of degeneration. The ROM at adjacent and index levels was extracted and normalized with respect to the healthy (baseline) spine.

Results:

A nonuniform distribution in the ROM was found for different disc degeneration states, segmental levels, and flexion–extension loading modes. The specific results for each and level are reported in the results section of the paper.

Conclusion:

Closer follow-up times may be necessary in symptomatic patients with progressive disease, especially with BL involvements.

A novel least-squares method to characterize in-vivo joint functional passive regional stiffness zones

OBJECTIVE

To present and evaluate a method to objectively quantify the functional regions of joint lumped passive stiffness.

BACKGROUND

Joint passive stiffness has an important clinical role in constraining the degrees of freedom at a given joint. Links between passive stiffness and injury, pathology and function may be better understood if joint passive stiffness can be accurately quantified. Thus, a technique was developed to objectively partition passive stiffness curves into 3 linear regions (low, transition and high stiffness).

METHODS

The passive stiffness of the lumbar spine is presented as an example. Simulated data was used to determine the sensitivity of the method to Gaussian white noise in force measurements. An experimentally determined lumbar passive flexion curve was used to demonstrate the technique on human data. Breakpoint analysis was employed on the resulting moment-angle cures to partition the curve into low, transition and high stiffness zones.

RESULTS

The proposed method was successful at discriminating between the three stiffness zones and quantifying the passive stiffness within each zone. The algorithm had difficulty determining parameters in the low-stiffness zone in the presence of noise.

CONCLUSIONS

The proposed method can be used as an objective method to investigate passive stiffness. Breakpoint Analysis can identify the three functional linear zones of passive stiffness. The slopes of these linear regions are then used as a measure of passive stiffness, which have applications in clinical populations and research studies, to assess time varying responses, or changes in stiffness following an intervention.

Development of a detailed human neck finite element model and injury risk curves under lateral impact

Advanced neck finite element modeling and development of neck injury criteria are important for the design of optimal neck protection systems in automotive and other environments. They are also important in virtual tests. The objectives of the present study were to develop a detailed finite element model (FEM) of the human neck and couple it to the existing head model, validate the model with kinematic data from legacy human volunteer and human cadaver impact datasets, and derive lateral impact neck injury risk curves using survival analysis from the upper and lower neck forces and moments. The detailed model represented the anatomy of a young adult mid-size male. It included all the cervical and first thoracic vertebrae, intervening discs, upper and lower spinal ligaments, bilateral facet joints, and passive musculature. Material properties were obtained from literature. Frontal, oblique, and lateral impacts to the distal end of the model was applied based on human volunteer and human cadaver experimental data. Corridor and cross-correlation methods were used for validation. The CORrelation and Analysis (CORA) score was used for objective assessments. Forces and moments were obtained at the occipital condyles (OC) and T1, and parametric survival analysis was used to derive injury risk curves to define human neck injury tolerance to lateral impact. The Brier Score Metric (BSM) was used to determine the hierarchical sequence among the injury metrics. The CORA scores for the lateral, frontal, and oblique impact loading conditions were 0.80, 0.91, and 0.87, respectively, for human volunteer data, and the mean score was 0.7 for human cadaver lateral impacts. Injury risk curves along with ±95% confidence intervals are given for all the four biomechanical metrics. The OC shear force was the optimal metric based on the BSM. A force of 1.5 kN was associated with the 50% probability level of AIS3+ neck injury. As a first step, the presented risk curves serve as human tolerance criteria under lateral impact, hitherto not available in published literatures, and they can be used in virtual testing and advancing restraint systems for improving human safety.

The rate of tendon failure in a collagen fibre recruitment-based model

Accurate characterization of the mechanical response of collagenous tissues is critical for investigations into mechanisms of soft tissue injury. These tissues are inherently viscoelastic, exhibiting strain-rate dependent stiffnesses, creep, and stress-relaxation. The strain-rate features of the failure portion of the stress-strain curve are less well developed. Collagen-distribution based models are improving and capable of reproducing the non-linear aspects of the elastic response of soft tissues, but still require parameterization of failure regions. Therefore, the purpose of this investigation, was to determine whether the parameters characterizing the rate of damage accumulation in a collagen-distribution model are proportional to strain rate. Fifty rat tail tendons were subjected to one of five different strain rates (0.01, 0.05, 0.1, 0.15, 0.20 s^-1) until failure in an uni-axial strain test. To test the hypothesis that the parameters associated with damage rate are proportional to strain rate, a collagen distribution model was employed with the parameters describing the rate of fibre damage being obtained by least-squares and regressed against the strain rate. The breaking function was found to be proportional to strain rate, with a proportionality constant of 60.7 s^-1. Properties characterizing the failure portion of the stress-strain curves for rat tail tendons are also reported. The Young's Modulus did not vary with strain rate and was found to be 103.3 ± 49.5 MPa. Similarly, failure stresses and strains did not vary across the strain rates tested, and were 15.6 ± 6.1 MPa and 32.2 ± 9.1%, respectively.

A detailed finite element model of a mid-sized male for the investigation of traffic pedestrian accidents

The pedestrian is one of the most vulnerable road users and comprises approximately 23% of the road crash-related fatalities in the world. To protect pedestrians during Car-to-Pedestrian Collisions (CPC), subsystem impact tests are used in regulations. These tests provide insight but cannot characterize the complex vehicle-pedestrian interaction. The main purpose of this study was to develop and validate a detailed pedestrian Finite Element (FE) model corresponding to a 50th percentile male to predict CPC induced injuries. The model geometry was reconstructed using a multi-modality protocol from medical images and exterior scan data corresponding to a mid-sized male volunteer. To investigate injury response, this model included internal organs, muscles and vessels. The lower extremity, shoulder and upper body of the model were validated against Post Mortem Human Surrogate (PMHS) test data in valgus bending, and lateral/anterior-lateral blunt impacts, respectively. The whole-body pedestrian model was validated in CPC simulations using a mid-sized sedan and simplified generic vehicles bucks and previously unpublished PMHS coronal knee angle data. In the component validations, the responses of the FE model were mostly within PMHS test corridors and in whole body validations the kinematic and injury responses predicted by the model showed similar trends to PMHS test data. Overall, the detailed model showed higher biofidelity, especially in the upper body regions, compared to a previously reported simplified pedestrian model, which recommends using it in future pedestrian automotive safety research.

An electromyography-assisted biomechanical cervical spine model: Model development and validation

BACKGROUND

In spite of the prevalence of occupational neck disorders as a result of work force fluctuating from industry to sedentary office work, most cervical spine computational models are not capable of simulating occupational and daily living activities whereas majority of cervical spine models specialized to simulate crash and impact scenarios. Therefore, estimating spine tissue loads accurately to quantify the risk of neck disorders in occupational environments within those models is not possible due to the lack of muscle models, dynamic simulation and passive spine structures.

METHODS

In this effort the structure, logic, and validation process of an electromyography-assisted cervical biomechanical model that is capable of estimating neck loading under three-dimensional complex motions is described. The developed model was designed to simulate complex dynamic motions similar to work place exposure. Curved muscle geometry, personalized muscle force parameters, and separate passive and (electromyography-driven) active muscle force components are implemented in this model.

FINDINGS

Calibration algorithms were able to reverse-engineer personalized muscle properties to calculate active and passive muscle forces of each individual.

INTERPRETATION

This electromyography-assisted cervical spine model with curved muscle model is capable to accurately predict spinal tissue loads during isometric and dynamic head and neck activities. Personalized active and passive muscle force algorithms will help to more robustly investigate person-specific muscle forces and spinal tissue loads.

A one-dimensional collagen-based biomechanical model of passive soft tissue with viscoelasticity and failure

INTRODUCTION

Strains and sprains of soft tissues, including tendons and ligaments, are frequently occurring injuries. Musculoskeletal models show great promise in prediction and prevention of these injuries. However, these models rarely account for the viscoelastic properties of ligaments and tendons, much less their failure properties. The purpose of this project was to develop, simplify, and analyze a collagen-distribution model to address these limitations.

MODEL DEVELOPMENT

A distribution-moment approximation was applied to an existing partial differential equation model to reduce its computational complexity. The resulting model was equipped with a Voigt model in series, which endowed it with viscoelastic properties in addition to failure properties.

RESULTS

The model was able to reproduce the characteristic toe, linear, and failure regions ubiquitous throughout in-vitro tests on tissue specimens. In addition, it was able to reproduce a tri-phasic creep test consisting of an initial deformation, a steady-state, and failure. Stress-relaxation and hysteresis were also reproducible by the model.

DISCUSSION AND CONCLUSION

The ability to reproduce so many characteristics of biological tissues suggests more bio-fidelity was achieved by the reduced model was other currently available models. Future work to further improve its bio-fidelity is proposed for specific tendons and ligaments.

Sex-Specific Intubation Biomechanics: Intubation Forces Are Greater in Male Than in Female Patients, Independent of Body Weight

Background

Studies of head, neck, and cervical spine morphology and tissue material properties indicate that cervical spine biomechanics differ between adult males and females. These differences result in sex-specific cervical spine kinematics and injury patterns in response to standardized loading conditions. Because direct laryngoscopy and endotracheal intubation require the application of a load to the cervical spine, intubation biomechanics should be sex-specific. The aim of this study was to determine if intubation forces during direct laryngoscopy differ between male and female patients and, if so, is the difference independent of body weight.

Methods

We pooled original data from three previously published adult clinical intubation studies that used methodologically reliable intubation force measurements (measured total laryngoscope force applied to the tongue, and force values were insensitive to or accounted for other laryngoscope blade forces). All patients had undergone direct laryngoscopy and orotracheal intubation with a Macintosh 3 blade under general anesthesia. Patient data included sex, age, height, weight, and maximum intubation force. Least squares multivariable linear regression was performed between the dependent variable (maximum intubation force) and two independent variables (patient sex and patient weight). A third term was added for the interaction between patient sex and weight.

Results

Among all patients (males n=42, females n=59), the median intubation force was 42.2 N (25^th to 75^th percentiles: 31.5 to 57.4 N). While controlling for patient body weight, intubation force differed between the sexes; P=0.011, with greater intubation force in male patients. While controlling for patient sex, there was a positive association between patient body weight and intubation force; P=0.009. In addition, there was a significant interaction between patient sex and weight; P=0.002, with intubation force in male patients having greater dependence on body weight. The difference in intubation force between male and female patients who had the same body weight exceeded 5 N when body weight exceeded 75 kg, and intubation force differences between male and female patients increased as patient body weight increased. Additional analyses using robust regression and using body mass index instead of weight provided comparable results.

Conclusion

In adult patients, the biomechanics of direct laryngoscopy and intubation are sex-specific. Our findings support the need to control for patient sex and weight in future clinical and laboratory studies of the human cervical spine and head and neck biomechanics.

Rear-Impact Neck Whiplash: Role of Head Inertial Properties and Spine Morphological Variations on Segmental Rotations

Whiplash injuries continue to be a concern in low-speed rear impact. This study was designed to investigate the role of variations in spine morphology and head inertia properties on cervical spine segmental rotation in rear-impact whiplash loading. Vertebral morphology is rarely considered as an input parameter in spine finite element (FE) models. A methodology toward considering morphological variations as input parameters and identifying the influential variations is presented in this paper. A cervical spine FE model, with its morphology parametrized using mesh morphing, was used to study the influence of disk height, anteroposterior vertebral depth, and segmental size, as well as variations in head mass, moment of inertia, and center of mass locations. The influence of these variations on the characteristic S-curve formation in whiplash response was evaluated using the peak C2–C3 flexion marking the maximum S-curve formation and time taken for the formation of maximum S-curve. The peak C2–C3 flexion in the S-curve formation was most influenced by disk height and vertebral depth, followed by anteroposterior head center of mass location. The time to maximum S-curve was most influenced by the anteroposterior location of head center of mass. The influence of gender-dependent variations, such as the vertebral depth, suggests that they contribute to the greater segmental rotations observed in females resulting in different S-curve formation from men. These results suggest that both spine morphology and head inertia properties should be considered to describe rear-impact responses.

Objective Evaluation of Whole Body Kinematics in a Simulated, Restrained Frontal Impact

The use of human body models as an additional data point in the evaluation of human-machine interaction requires quantitative validation. In this study a validation of the Global Human Body Models Consortium (GHBMC) average male occupant model (M50-O v. 4.5) in a restrained frontal sled test environment is presented. For vehicle passengers, frontal crash remains the most common mode, and the most common source of fatalities. A total of 55-time history traces of reaction loads and kinematics from the model were evaluated against corresponding PMHS data (n = 5). Further, the model's sensitivity to the belt path was studied by replicating two documented PMHS cases with prominent lateral and medial belt paths respectively. Results were quantitatively evaluated using open source CORA software. A tradeoff was observed; better correlation scores were achieved on gross measures (e.g. reaction loads), whereas better corridor scores were achieved on localized measures (rib deflections), indicating that subject specificity may dominate the comparison at localized anatomical regions. On an overall basis, the CORA scores were 0.68, 0.66 and 0.60 for force, body kinematics and chest wall kinematics. Belt force responses received the highest grouped CORA score of 0.85. Head and sternum kinematics earning a 0.8 and 0.7 score respectively. The model demonstrated high sensitivity to belt path, resulting in a 20-point increase in CORA score when the belt was routed closer to analogous location of data collection. The human model demonstrated overall reasonable biofidelity and sensitivity to countermeasures in frontal crash kinematics.

Utilization of 3D printing technology to facilitate and standardize soft tissue testing

Three-dimensional (3D) printing has become broadly available and can be utilized to customize clamping mechanisms in biomechanical experiments. This report will describe our experience using 3D printed clamps to mount soft tissues from different anatomical regions. The feasibility and potential limitations of the technology will be discussed. Tissues were sourced in a fresh condition, including human skin, ligaments and tendons. Standardized clamps and fixtures were 3D printed and used to mount specimens. In quasi-static tensile tests combined with digital image correlation and fatigue trials we characterized the applicability of the clamping technique. Scanning electron microscopy was utilized to evaluate the specimens to assess the integrity of the extracellular matrix following the mechanical tests. 3D printed clamps showed no signs of clamping-related failure during the quasi-static tests, and intact extracellular matrix was found in the clamping area, at the transition clamping area and the central area from where the strain data was obtained. In the fatigue tests, material slippage was low, allowing for cyclic tests beyond 10⁵ cycles. Comparison to other clamping techniques yields that 3D printed clamps ease and expedite specimen handling, are highly adaptable to specimen geometries and ideal for high-standardization and high-throughput experiments in soft tissue biomechanics.

A preliminary study of human model head and neck response to frontal loading in nontraditional occupant seating configurations

OBJECTIVE

Computational human body models (HBMs) are nominally omnidirectional surrogates given their structural basis in human anatomy. As a result, such models are well suited for studies related to occupant safety in anticipated highly automated vehicles (HAVs). We utilize a well-validated HBM to study the head and neck kinematics in simulations of nontraditional occupant seating configurations.

METHODS

The GHBMC M50-O v. 4.4 HBM was gravity settled into a generic seat buck and situated in a seated posture. The model was simulated in angular increments of 15 degrees clockwise from forward facing to rear facing. A pulse of 17.0 kph (NASS median) was used in each to simulate a frontal impact for each of the 13 seating configurations. Belt anchor points were rotated with the seat; the airbag was appropriately powered based on delta-V, and was not used in rear-facing orientations. Neck forces and moments were calculated.

RESULTS

The 30-degree oblique case was found to result in the maximum neck load and sagittal moment, and thus Neck Injury Criteria (NIJ). Neck loads were minimized in the rear facing condition. The moments and loads, however, were greatest in the lateral seating configuration for these frontal crash simulations.

CONCLUSIONS

In a recent policy statement on HAVs, the NHTSA indicated that vehicle manufacturers will be expected to provide countermeasures that will fully protect occupants given any planned seating or interior configurations. Furthermore, the agency indicated that virtual tests using human models could be used to demonstrate such efficacy. While the results presented are only appropriate for comparison within this study, they do indicate that human models provide reasonable biomechanical data for nontraditional occupant seating arrangements.

Failed rib region prediction in a human body model during crash events with precrash braking

OBJECTIVE

The objective of this study is 2-fold. We used a validated human body finite element model to study the predicted chest injury (focusing on rib fracture as a function of element strain) based on varying levels of simulated precrash braking. Furthermore, we compare deterministic and probabilistic methods of rib injury prediction in the computational model.

METHODS

The Global Human Body Models Consortium (GHBMC) M50-O model was gravity settled in the driver position of a generic interior equipped with an advanced 3-point belt and airbag. Twelve cases were investigated with permutations for failure, precrash braking system, and crash severity. The severities used were median (17 kph), severe (34 kph), and New Car Assessment Program (NCAP; 56.4 kph). Cases with failure enabled removed rib cortical bone elements once 1.8% effective plastic strain was exceeded. Alternatively, a probabilistic framework found in the literature was used to predict rib failure. Both the probabilistic and deterministic methods take into consideration location (anterior, lateral, and posterior). The deterministic method is based on a rubric that defines failed rib regions dependent on a threshold for contiguous failed elements. The probabilistic method depends on age-based strain and failure functions.

RESULTS

Kinematics between both methods were similar (peak max deviation: ΔX_head = 17 mm; ΔZ_head = 4 mm; ΔX_thorax = 5 mm; ΔZ_thorax = 1 mm). Seat belt forces at the time of probabilistic failed region initiation were lower than those at deterministic failed region initiation. The probabilistic method for rib fracture predicted more failed regions in the rib (an analog for fracture) than the deterministic method in all but 1 case where they were equal. The failed region patterns between models are similar; however, there are differences that arise due to stress reduced from element elimination that cause probabilistic failed regions to continue to rise after no deterministic failed region would be predicted.

CONCLUSIONS

Both the probabilistic and deterministic methods indicate similar trends with regards to the effect of precrash braking; however, there are tradeoffs. The deterministic failed region method is more spatially sensitive to failure and is more sensitive to belt loads. The probabilistic failed region method allows for increased capability in postprocessing with respect to age. The probabilistic failed region method predicted more failed regions than the deterministic failed region method due to force distribution differences.

Incorporating ligament laxity in a finite element model for the upper cervical spine

BACKGROUND CONTEXT

Predicting physiological range of motion (ROM) using a finite element (FE) model of the upper cervical spine requires the incorporation of ligament laxity. The effect of ligament laxity can be observed only on a macro level of joint motion and is lost once ligaments have been dissected and preconditioned for experimental testing. As a result, although ligament laxity values are recognized to exist, specific values are not directly available in the literature for use in FE models.

PURPOSE

The purpose of the current study is to propose an optimization process that can be used to determine a set of ligament laxity values for upper cervical spine FE models. Furthermore, an FE model that includes ligament laxity is applied, and the resulting ROM values are compared with experimental data for physiological ROM, as well as experimental data for the increase in ROM when a Type II odontoid fracture is introduced.

DESIGN/SETTING

The upper cervical spine FE model was adapted from a 50th percentile male full-body model developed with the Global Human Body Models Consortium (GHBMC). FE modeling was performed in LS-DYNA and LS-OPT (Livermore Software Technology Group) was used for ligament laxity optimization.

METHODS

Ordinate-based curve matching was used to minimize the mean squared error (MSE) between computed load-rotation curves and experimental load-rotation curves under flexion, extension, and axial rotation with pure moment loads from 0 to 3.5 Nm. Lateral bending was excluded from the optimization because the upper cervical spine was considered to be primarily responsible for flexion, extension, and axial rotation. Based on recommendations from the literature, four varying inputs representing laxity in select ligaments were optimized to minimize the MSE. Funding was provided by the Natural Sciences and Engineering Research Council of Canada as well as GHMBC. The present study was funded by the Natural Sciences and Engineering Research Council of Canada to support the work of one graduate student. There are no conflicts of interest to be reported.

RESULTS

The MSE was reduced to 0.28 in the FE model with optimized ligament laxity compared with an MSE 0f 4.16 in the FE model without laxity. In all load cases, incorporating ligament laxity improved the agreement between the ROM of the FE model and the ROM of the experimental data. The ROM for axial rotation and extension was within one standard deviation of the experimental data. The ROM for flexion and lateral bending was outside one standard deviation of the experimental data, but a compromise was required to use one set of ligament laxity values to achieve a best fit to all load cases. Atlanto-occipital motion was compared as a ratio to overall ROM, and only in extension did the inclusion of ligament laxity not improve the agreement. After a Type II odontoid fracture was incorporated into the model, the increase in ROM was consistent with experimental data from the literature.

CONCLUSIONS

The optimization approach used in this study provided values for ligament laxities that, when incorporated into the FE model, generally improved the ROM response when compared with experimental data. Successfully modeling a Type II odontoid fracture showcased the robustness of the FE model, which can now be used in future biomechanics studies.

The effect of precrash velocity reduction on occupant response using a human body finite element model

OBJECTIVE

The objective of this study is to use a validated finite element model of the human body and a certified model of an anthropomorphic test dummy (ATD) to evaluate the effect of simulated precrash braking on driver kinematics, restraint loads, body loads, and computed injury criteria in 4 commonly injured body regions.

METHODS

The Global Human Body Models Consortium (GHBMC) 50th percentile male occupant (M50-O) and the Humanetics Hybrid III 50th percentile models were gravity settled in the driver position of a generic interior equipped with an advanced 3-point belt and driver airbag. Fifteen simulations per model (30 total) were conducted, including 4 scenarios at 3 severity levels: median, severe, and the U.S. New Car Assessment Program (U.S.-NCAP) and 3 extra per model with high-intensity braking. The 4 scenarios were no precollision system (no PCS), forward collision warning (FCW), FCW with prebraking assist (FCW+PBA), and FCW and PBA with autonomous precrash braking (FCW + PBA + PB). The baseline ΔV was 17, 34, and 56.4 kph for median, severe, and U.S.-NCAP scenarios, respectively, and were based on crash reconstructions from NASS/CDS. Pulses were then developed based on the assumed precrash systems equipped. Restraint properties and the generic pulse used were based on literature.

RESULTS

In median crash severity cases, little to no risk (<10% risk for Abbreviated injury Scale [AIS] 3+) was found for all injury measures for both models. In the severe set of cases, little to no risk for AIS 3+ injury was also found for all injury measures. In NCAP cases, highest risk was typically found with No PCS and lowest with FCW + PBA + PB. In the higher intensity braking cases (1.0-1.4 g), head injury criterion (HIC), brain injury criterion (BrIC), and chest deflection injury measures increased with increased braking intensity. All other measures for these cases tended to decrease. The ATD also predicted and trended similar to the human body models predictions for both the median, severe, and NCAP cases. Forward excursion for both models decreased across median, severe, and NCAP cases and diverged from each other in cases above 1.0 g of braking intensity.

CONCLUSIONS

The addition of precrash systems simulated through reduced precrash speeds caused reductions in some injury criteria, whereas others (chest deflection, HIC, and BrIC) increased due to a modified occupant position. The human model and ATD models trended similarly in nearly all cases with greater risk indicated in the human model. These results suggest the need for integrated safety systems that have restraints that optimize the occupant's position during precrash braking and prior to impact.

A Female Ligamentous Cervical Spine Finite Element Model Validated for Physiological Loads

Mathematical cervical spine models allow for studying of impact loading that can cause whiplash associated disorders (WAD). However, existing models only cover the male anthropometry, despite the female population being at a higher risk of sustaining WAD in automotive rear-end impacts. The aim of this study is to develop and validate a ligamentous cervical spine intended for biomechanical research on the effect of automotive impacts. A female model has the potential to aid the design of better protection systems as well as improve understanding of injury mechanisms causing WAD. A finite element (FE) mesh was created from surface data of the cervical vertebrae of a 26-year old female (stature 167 cm, weight 59 kg). Soft tissues were generated from the skeletal geometry and anatomical literature descriptions. Ligaments were modeled with nonlinear elastic orthotropic membrane elements, intervertebral disks as composites of nonlinear elastic bulk elements, and orthotropic anulus fibrosus fiber layers, while cortical and trabecular bones were modeled as isotropic plastic-elastic. The model has geometrical features representative of the female cervical spine-the largest average difference compared with published anthropometric female data was the vertebral body depth being 3.4% shorter for the model. The majority the cervical segments compare well with respect to biomechanical data at physiological loads, with the best match for flexion-extension loads and less biofidelity for axial rotation. An average female FE ligamentous cervical spine model was developed and validated with respect to physiological loading. In flexion-extension simulations with the developed female model and an existing average male cervical spine model, a greater range of motion (ROM) was found in the female model.

Planar biaxial extension of the lumbar facet capsular ligament reveals significant in-plane shear forces

The lumbar facet capsular ligament (FCL) articulates with six degrees of freedom during spinal motions of flexion/extension, lateral bending, and axial rotation. The lumbar FCL is composed of highly aligned collagen fiber bundles on the posterior surface (oriented primarily laterally between the rigid articular facets) and irregularly oriented elastin on the anterior surface. Because the FCL is a capsule, it has multiple insertion sites across the lumbar facet joint, which, along with its material structure, give rise to complicated deformations in vivo. We performed planar equibiaxial mechanical tests on excised healthy cadaveric lumbar FCLs (n=6) to extract normal and shear reaction forces, and fit sample-specific two-fiber-family finite element models to the experimental force data. An eight-parameter anisotropic, hyperelastic model was used. Shear forces at maximum extension (mean values of 1.68N and 3.01N in the two directions) were of comparable magnitude to the normal forces perpendicular to the aligned collagen fiber bundles (4.67N) but smaller than normal forces in the fiber direction (16.11N). Inclusion of the experimental shear forces in the model optimization yielded fits with highly aligned fibers oriented at a specific angle across all samples, typically with one fiber population aligned nearly horizontally and the other at an oblique angle. Conversely, models fit to only the normal force data resulted in a broad range of fiber angles with low specificity. We found that shear forces generated through planar equibiaxial extension aided the model fit in describing the anisotropic nature of the FCL surface.

Out-of-Position Rear Impact Tissue-Level Investigation Using Detailed Finite Element Neck Model

OBJECTIVE

Whiplash injuries can occur in automotive crashes and may cause long-term health issues such as neck pain, headache, and visual and auditory disturbance. Evidence suggests that nonneutral head posture can significantly increase the potential for injury in a given impact scenario, but epidemiological and experimental data are limited and do not provide a quantitative assessment of the increased potential for injury. Although there have been some attempts to evaluate this important issue using finite element models, none to date have successfully addressed this complex problem.

METHODS

An existing detailed finite element neck model was evaluated in nonneutral positions and limitations were identified, including musculature implementation and attachment, upper cervical spine kinematics in axial rotation, prediction of ligament failure, and the need for repositioning the model while incorporating initial tissue strains. The model was enhanced to address these issues and an iterative procedure was used to determine the upper cervical spine ligament laxities. The neck model was revalidated using neutral position impacts and compared to an out-of-position cadaver experiment in the literature. The effects of nonneutral position (axial head rotation) coupled with muscle activation were studied at varying impact levels.

RESULTS

The laxities for the ligaments of the upper cervical spine were determined using 4 load cases and resulted in improved response and predicted failure loads relative to experimental data. The predicted head response from the model was similar to an experimental head-turned bench-top rear impact experiment. The parametric study identified specific ligaments with increased distractions due to an initial head-turned posture and the effect of active musculature leading to reduced ligament distractions.

CONCLUSIONS

The incorporation of ligament laxity in the upper cervical spine was essential to predict range of motion and traumatic response, particularly for repositioning of the neck model prior to impact. The results of this study identify a higher potential for injury in out-of-position rear collisions and identified at-risk locations based on ligament distractions. The model predicted higher potential for injury by as much as 50% based on ligament distraction for the out-of-position posture and reduced potential for injury with muscle activation. Importantly, this study demonstrated that the location of injury or pain depends on the initial occupant posture, so that both the location of injury and kinematic threshold may vary when considering common head positions while driving.

Development of a computationally efficient full human body finite element model

INTRODUCTION

A simplified and computationally efficient human body finite element model is presented. The model complements the Global Human Body Models Consortium (GHBMC) detailed 50th percentile occupant (M50-O) by providing kinematic and kinetic data with a significantly reduced run time using the same body habitus.

METHODS

The simplified occupant model (M50-OS) was developed using the same source geometry as the M50-O. Though some meshed components were preserved, the total element count was reduced by remeshing, homogenizing, or in some cases omitting structures that are explicitly contained in the M50-O. Bones are included as rigid bodies, with the exception of the ribs, which are deformable but were remeshed to a coarser element density than the M50-O. Material models for all deformable components were drawn from the biomechanics literature. Kinematic joints were implemented at major articulations (shoulder, elbow, wrist, hip, knee, and ankle) with moment vs. angle relationships from the literature included for the knee and ankle. The brain of the detailed model was inserted within the skull of the simplified model, and kinematics and strain patterns are compared.

RESULTS

The M50-OS model has 11 contacts and 354,000 elements; in contrast, the M50-O model has 447 contacts and 2.2 million elements. The model can be repositioned without requiring simulation. Thirteen validation and robustness simulations were completed. This included denuded rib compression at 7 discrete sites, 5 rigid body impacts, and one sled simulation. Denuded tests showed a good match to the experimental data of force vs. deflection slopes. The frontal rigid chest impact simulation produced a peak force and deflection within the corridor of 4.63 kN and 31.2%, respectively. Similar results vs. experimental data (peak forces of 5.19 and 8.71 kN) were found for an abdominal bar impact and lateral sled test, respectively. A lateral plate impact at 12 m/s exhibited a peak of roughly 20 kN (due to stiff foam used around the shoulder) but a more biofidelic response immediately afterward, plateauing at 9 kN at 12 ms. Results from a frontal sled simulation showed that reaction forces and kinematic trends matched experimental results well. The robustness test demonstrated that peak femur loads were nearly identical to the M50-O model. Use of the detailed model brain within the simplified model demonstrated a paradigm for using the M50-OS to leverage aspects of the M50-O. Strain patterns for the 2 models showed consistent patterns but greater strains in the detailed model, with deviations thought to be the result of slightly different kinematics between models. The M50-OS with the deformable skull and brain exhibited a run time 4.75 faster than the M50-O on the same hardware.

CONCLUSIONS

The simplified GHBMC model is intended to complement rather than replace the detailed M50-O model. It exhibited, on average, a 35-fold reduction in run time for a set of rigid impacts. The model can be used in a modular fashion with the M50-O and more broadly can be used as a platform for parametric studies or studies focused on specific body regions.

Application of Radial Basis Function Methods in the Development of a 95th Percentile Male Seated FEA Model

Human body finite element models (FEMs) are a valuable tool in the study of injury biomechanics. However, the traditional model development process can be time-consuming. Scaling and morphing an existing FEM is an attractive alternative for generating morphologically distinct models for further study. The objective of this work is to use a radial basis function to morph the Global Human Body Models Consortium (GHBMC) average male model (M50) to the body habitus of a 95th percentile male (M95) and to perform validation tests on the resulting model. The GHBMC M50 model (v. 4.3) was created using anthropometric and imaging data from a living subject representing a 50th percentile male. A similar dataset was collected from a 95th percentile male (22,067 total images) and was used in the morphing process. Homologous landmarks on the reference (M50) and target (M95) geometries, with the existing FE node locations (M50 model), were inputs to the morphing algorithm. The radial basis function was applied to morph the FE model. The model represented a mass of 103.3 kg and contained 2.2 million elements with 1.3 million nodes. Simulations of the M95 in seven loading scenarios were presented ranging from a chest pendulum impact to a lateral sled test. The morphed model matched anthropometric data to within a rootmean square difference of 4.4% while maintaining element quality commensurate to the M50 model and matching other anatomical ranges and targets. The simulation validation data matched experimental data well in most cases.