Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Bioeng. Biotechnol., 03 September 2025

Sec. Biomechanics

Volume 13 - 2025 | https://doi.org/10.3389/fbioe.2025.1597267

A biomechanical study of neck strength and impact dynamics on head and neck injury parameters

Rahid ZamanRahid ZamanMohammad Ibrahim HossainMohammad Ibrahim HossainAhmed Zubayer RaiyanAhmed Zubayer RaiyanShuvo ChowdhuryShuvo ChowdhuryAaron JacksonAaron JacksonArthur Thomas KosterArthur Thomas KosterAshfaq Adnan
Ashfaq Adnan*
  • Department of Mechanical and Aerospace Engineering, University of Texas at Arlington, Arlington, TX, United States

Introduction: Head and neck injuries, including traumatic brain injuries (TBI), are a leading cause of disability and death worldwide. It affects millions of people worldwide, from automobiles to sports to military personnel. This study investigates the influence of impact locations, severities, and neck strength on head and neck injury parameters using a musculoskeletal head-neck model in OpenSim software.

Methods: We hypothesize that eccentric impacts, particularly those on the anterolateral side, increase GAMBIT and Neck Injury Criteria (NIC) due to elevated rotational accelerations, and that higher neck strength mitigates GAMBIT and NIC under these impacts. To test our hypotheses, we investigated a total of 63 cases in which seven impact locations (two from the anterior side, two from the posterior side, and three from lateral sides), three neck strengths (low, mid, high strength capacity), and three impact severities (low, moderate, and high) were explored. Seven output parameters were analyzed: linear and rotational accelerations, the Generalized Acceleration Model for Brain Injury Threshold (GAMBIT), neck force, neck moment, and Neck Injury Criteria (NIC) and neck muscle strain.

Results: Results reveal that anterolateral eccentric impacts pose the greatest risk, with rotational acceleration reaching 4,176 that is 4.75 times higher than anterior central impacts (879 rad/s2). GAMBIT values for moderate and high severity impacts are 1.44 and 1.54 times greater than low severity impacts, respectively. Head and neck injury parameters vary minimally (10) with neck strength.

Discussion: In summary, the severities and location of impacts had a significant role in GAMBIT and NIC, and the anterolateral eccentric impact had a higher probability of head and neck injury than the other six impact locations. These findings underscore the critical role of impact location and severity in injury risk and suggest helmet padding in lateral and anterolateral zones with energy-absorbing materials to reduce rotational acceleration.

1 Introduction

An estimated 1.6 to 3.8 million sport and recreation-related traumatic brain injuries (TBIs) occur in the USA annually (Borich et al., 2013). Over 6 million passenger car accidents take place every year, and more than 42,514 people lost their lives in 2022 in the USA (Insurance Institute for Highway Safety and Highway Loss Data Institute, 2024). According to the CDC, in 2019, 60,611 TBI-related deaths occurred in the USA (Centers for Disease Control and Prevention, 2019). Along with the TBI, neck injuries are also common in high-energy events like blasts, automobile accidents, and sports impacts. It is reported that half a million people faced cervical spine injuries in high-energy accidents in the USA from 2005 to 2013, and the trend of incidents has increased from 4.1% to 5.4% (Johnson et al., 2020). A major portion of these head and neck injury events has been widely underreported. Even in fully equipped settings, like sports, 66% of TBI go unreported (McCrea et al., 2004). According to the Defense Health Agency, 505,896 service members faced mild to severe head and neck injuries from 2000 to 2024 (Defense Health Agency, 2024). In a study, it was found that 57% of service members who faced head and neck injuries did not seek medical care (Escolas et al., 2020). So, the percentage of the civil population that ignores head and neck injuries is really a concern. Moreover, the complications from TBI and neck injuries have long-term effects and can alter a person’s thinking, sensation, language, or emotions, which may not be immediately apparent. That’s why TBI and neck injuries are frequently referred to as a “silent epidemic” (Takhounts et al., 2011).

TBI is a complex pathophysiological process induced by mechanical loading of the brain. TBI can occur because of a direct impact or blow to the head or an indirect inertial response due to a non-cranial impact on the human body. Most of the life-threatening head and neck injuries from contact sports, automobiles, and blasts involve the torso, spine, and head. The rotational response of the head has been correlated to the risk of TBIs such as subdural hematomas and Diffuse Axonal Injury (DAI) (Browne et al., 2011). The effect of neck strength on concussions is still under investigation. Neck resistance and strength have an effect on the human head response during and after an impact (Farmer et al., 2022). Additionally, certain brain structures exhibit directional response dependencies, such as the falx and tentorium cerebri, particularly in the corpus callosum and brainstem, which undergo significant strain during axial rotations (Ho et al., 2017). Most studies indicate that stronger and stiffer necks reduce the risk of concussion (Fanta et al., 2013; Eckner et al., 2014; Hasegawa et al., 2014; Mang et al., 2015). However, some studies disagree with this (Mihalik et al., 2011). To gain a comprehensive understanding of the effects of neck strength and impact dynamics on head and neck injuries, three major parameters need to be carefully investigated. First, a computational or experimental head-neck model that can mimic head-neck response with high fidelity. Second, the input data for the analysis should be free from bias. Third, impacts from all around the head should be considered, as they can cause head and neck injuries.

A very detailed model of the human head and neck is a vital tool to assess any injury. The cervical spine is a pivotal structure of head dynamics and controls the head’s response after external loading. It establishes the mechanical linkage between the torso and head, and stabilizes the head using cervical disc and neck muscle tissue resistance. The most popular experimental and computational tools for reconstructing head and neck injury events are anthropomorphic test devices (ATD), multibody models, finite element models (FEM), and musculoskeletal models. Since all ATD use pin joints to connect the vertebrae, it allows limited rotation and bending in a single plane (Farmer et al., 2022). ATDs have 2-10 times stiffer necks and are mainly used for frontal, rear, and near-side automotive impact (Farmer et al., 2022). It lacks the natural compliance and multi-planar motion of human intervertebral discs and ligaments. ATDs can not mimic the intervertebral rotation with high biofidelity. Therefore, the evaluation of the rotational acceleration-based head injury criterion and Neck Injury Criterion (NIC) is not fully reliable with ATD experiments. Moreover, ATDs are expensive for repetitive blast and severe impact experiments. For these reasons, there is a growing interest in using fast and high bio-fidelity computational models. At the initial stages, lumped parameter models (Bumberger et al., 2025; Merrill et al., 1984; Deng and Goldsmith, 1987) and multibody models (Garcia and Ravani, 2003; Huo and Wu, 2017) of head and neck were developed. Most of those models are based on simplified rigid linkages. With the increase in computational efficiency, finite element models (FEM) became popular (Kleiven, 2007; Pramudita et al., 2009; Zhang et al., 2011; Gabler et al., 2019; Wu et al., 2019; Dymek et al., 2021; Liu et al., 2022; Putra et al., 2022; Singh et al., 2022; Tang et al., 2023; Wang et al., 2023). Finite element models were widely used to study the head and neck motion for whiplash (Zhang et al., 2011; Putra et al., 2022). However, FEM models are computationally very heavy and require predefined boundary conditions. Also, it is very difficult for the FEM model to exactly mimic the nature of muscle tissues and their interactions with the surrounding anatomical structures, such as muscle wrapping, dynamics, muscle excitations and contractions, and interactions with tendons and skeletons. The reliability of muscle force-generating capacity is a contributing factor in determining the head-neck model response. Not only does the shape and alignment of the cervical spine stabilize the head-neck joints, but the ligaments and voluntarily controlled muscles also play a contributing factor (Fanta et al., 2013). To understand the response of the head and neck under external loadings, it is essential to integrate all these physiological and biomechanical properties into the model to estimate a model’s muscle force-generating capacity. However, FEM model lacks those properties, hindering its ability to accurately simulate anatomical interactions.

Computational musculoskeletal models offer enhanced biofidelity compared to FEM and multibody models, as they incorporate dynamic inputs, thereby reducing reliance on static assumptions. These steps are critical for accurately capturing the biomechanical response to impact (Halloran et al., 2010). Musculoskeletal models incorporate time-dependent dynamic inputs such as muscle activations, contractions, wrapping, tendon compliances, and joint movement. Furthermore, computational musculoskeletal models are easily scalable to subject-specific anthropometric details and can provide segmental acceleration, velocity, and force information for a specific subject. These features reduce the kinematic errors at the early stage that can propagate throughout a biomechanical analysis. Musculoskeletal models can help us study the cause-and-effect relationship between various muscle and joint physiological properties and head and neck injuries, which is a vital part of this study. Additionally, computational musculoskeletal models are less time-consuming and computationally less expensive than FEM, cadaveric, or dummy tests. Considering all these factors and hypotheses of this study, we employed a musculoskeletal model to investigate the effects of neck strength and impact dynamics on head and neck injury criteria. In our model, we considered three neck strengths: low force capacity (sedentary or untrained), mid force capacity (recreationally active), and high force capacity (Athletic).

To perform a detailed injury analysis, sources of kinematics and kinetics data should be free from bias and cover a broad spectrum of events. The acceleration loading curve on the head has a significant effect on TBI (Post et al., 2014). Musculoskeletal models have been widely used for low-speed human movements like weight lifting, walking, running, etc. Recently, there have been some studies where the musculoskeletal model has been used for moderate severity impacts, such as NFL (Kleiven, 2007; Jin et al., 2017; Mortensen et al., 2018; Kuo et al., 2019; Mortensen et al., 2020) and rugby impacts (Silvestros et al., 2024). Head and neck injuries, because of low and high severity impacts, were not studied much using computational musculoskeletal models. In all those studies, sports impacts were the only source of real-life input data to study the effect of impact dynamics and neck strength on head and neck injury. However, sports data has visual perception and risk compensation phenomenon error (Schmidt et al., 2014). That means players with stronger and larger cervical musculature may be confident that they are more protected from head and neck injuries and may engage in more violent collisions. On the other hand, weaker subjects tend to avoid violent collisions. Furthermore, it was also reported that without visual perception, muscle activation takes 0.027s after an impact, while with visual perception, muscle activation takes 0.127s before the impact (Fanta et al., 2013). It is also reported in the same study that head injury criteria show about 30% decrease with visual perception. For these reasons, the studies and decisions about the effect of impact dynamics and neck strength on head injuries, based only on sports data, may mislead the injury prediction. To fill up these gaps, in our study, we consider three types of impact data: automobile impact test as low severity impact (LSI), NFL as moderate severty impact (MSI), and blast as high severity impact (HSI). Thus, we make our data independent of visual perception error. Details about the severity classification of impact are mentioned in Section 2.3.1.

Another research gap in the literature is that the impacts from all around the head are not widely studied for both head and neck injuries in terms of neck strength. Impact locations influence the head and neck injury metric (Mychasiuk et al., 2016; Yu et al., 2022). Impacts on the frontal or occipital region increase the risk of subdural hematoma, whereas parenchymal contusions are more likely due to the side impacts on the head (Post et al., 2015). In a study, the effects of anterior central, posterior central, lateral central, and posterolateral impacts were analyzed for sports-related concussive and sub-concussive impacts (Mortensen et al., 2020). They found that neck strength has statistical significance in terms of head injury criteria. In another study, they considered cranial anterior, cranial posterior, lateral mid-posterior, mid-anterior, and lateral interior impacts from rugby players (Silvestros et al., 2024). However, they studied only neck injuries; head injuries were not included. Anterior eccentric and posterior eccentric impacts were not studied for both head and neck injuries. However, our hypothesis is that eccentric impact may create higher accelerations and momentum around the head than concentric impacts. To fix this research gap, we considered seven impact locations for all the above-mentioned impact severities and neck strengths: anterior concentric, anterior eccentric, posterior concentric, posterior eccentric, lateral concentric, posterolateral eccentric, and anterolateral eccentric.

To investigate the effect of neck strength and impact dynamics on head and neck injury, a comprehensive study is necessary to examine different neck strengths and consider varying impact severities from all directions to reach a conclusive decision. We hypothesize that eccentric impacts will result in higher GAMBIT and NIC values compared to concentric impacts due to the increased rotational accelerations associated with eccentric loading. In this study, we designed 63 cases guided by these hypotheses (3 impact severities × 3 neck strengths × 7 impact locations). In Section 2, we will provide details about the head-neck dynamics framework, injury criteria, impact dynamics, and simulation workflow. In Section 3, we will present the results, and in Section 4, we will discuss the influence of impact locations, neck strengths, and impact severities on head and neck injury criteria.

2 Methods

The musculoskeletal head-neck model is a muscle-driven model. Neck muscles generate forces based on complex head-neck muscle activation and musculo-tendon dynamics. Under an accelerative environment, all neck muscles collectively generate a muscle-driven head-neck response (head acceleration, muscle strain, and joint reaction load of cervical spine), which is crucial in head and neck injury prediction with high bio-fidelity. The head-neck dynamics framework is provided below in Figure 1.

Figure 1
Diagram of a head-neck muscle dynamics model illustrating the interaction between external forces, neck muscle dynamics, and forward dynamics. It shows muscle activation and contraction processes, joint geometry, and moment arms. Inset includes a detailed schematic of muscle-tendon dynamics with force equations. Graphs display muscle force-length, force-velocity relations, and tendon strain. Annotations define variables such as muscle excitations, activations, forces, joint torques, and joint accelerations, velocities, and angles. The illustration aims to model physiological responses to impacts.

Figure 1. Head-neck dynamics framework uses a muscle-driven simulation to determine the head accelerations. Impact direction (θ), impact severiy (Γ), and Eccentric multiplier (EM) are considered as design variables of this simulation (yellow color).

2.1 Head-neck dynamics framework

2.1.1 Head-neck dynamics

The dynamics of the human musculoskeletal system can be formulated with the Euler-Lagrange equations (Zuo et al., 2024) as per Equation 1:

Mqq¨=Cq,q˙+Gq+RqFTu+Eq,q˙(1)

Where, q, q˙, and q¨ are the positions, velocities, and accelerations of the joints, Mq is the mass distribution matrix which contains masses and inertial properties of the body segments, Cq,q˙ is the Coriolis and centrifugal force vector which arises when Newton’s laws of motion are applied in reference frames that are fixed to rotating bodies, Gq is the gravitational force vector, Rq is muscle moment arm matrix, FTu is the tendon force vectors, which is a function of muscle excitations (u), Eq,q˙ is the external forces vector that represents the interactions between the body and environment, including the impact direction (θ) and impact severities (Γ).

The acceleration of the body in response to muscle forces and other loads can be computed using the equations of motion for the body as per Equation 2.

q¨=M1qCq,q˙+Gq+RqFTu+Eq,q˙(2)

The process to get the muscle force vectors is shown in Equation 4.

2.1.2 Muscle-force generating parameters

Muscle activation dynamics takes muscle excitation as the input and muscle activation as the output. Muscle excitation ut represents the strength of the excitation signal from the nerve to the muscle, whereas muscle activation at represents the availability of calcium ions within intracellular space (Uchida and Delp, 2021). The relationship between muscle activation at and muscle excitation ut can be expressed as a first-order ordinary differential equation, as shown in Equation 3:

a˙=utatτa,u(3)
whereτa,u=τA0.5+1.5at,ifut>atτD0.5+1.5at,otherwise

τA and τD represent the activation and deactivation time constants. τA is smaller that τD. Although the values vary based on age, muscle, composition and other factors, typical values are 10 and 40 ms, respectively (Uchida and Delp, 2021). This relationship shows that the rate of activation slows as the activation level increases due to less amount of calcium release and diffusion. Similarly, the rate of deactivation slows as the muscle activation level decreases. The force-length-velocity relationship of a muscle largely depends on the muscle activation dynamics.

Neck muscle forces can be computed from the activation at, normalized fiber length lM), and normalized muscle velocity (vM). Muscle force at particular muscle length, velocity and activation is the product of the corresponding values on the force-length and force-velocity curves, as shown in Equation 4:

fMt=foMatfLlMtfvvMt+fPElMt(4)

Where fMt,foM,anda are fiber force, maximum isometric force, and muscle activation, respectively.

The total generated neck muscle force is a function of muscle fiber length and muscle velocity. The force-length curve, fLlMt provides the relationship between neck muscle length and generated force at a specific time. The force-velocity curve fvvMt provides the relationship between neck muscle velocity and generated force at a specific time. As our impact severities are classified based on force and duration, we considered neck strength to be a function of the force-velocity curve fvvMt. In musculoskeletal modeling, maximum force output can be expressed as an Eccentric Multiplier (EM) which represents the maximum eccentric contraction velocity, as shown in Figure 2. This parameter ranges between 1.1 and 1.8 and varies from subject to subject.

Figure 2
Graph showing the relationship between normalized force and normalized velocity. Three curves represent different force capacities: high (green), intermediate (blue), and low (red). The graph includes illustrations of muscle shortening for negative velocity and muscle lengthening for positive velocity.

Figure 2. Force-velocity relationship. A higher eccentric multiplier results in higher forces during eccentric (lengthening) contractions. Positive velocity indicates lengthening, and negative velocity indicates shortening.

Based on the eccentric multiplier, we consider three subjects (neck strengths): (i) Low force capacity (Sedentary or untrained): EM = 1.2, (ii) Mid force capacity (Recreationally active): EM = 1.4, and (iii) High force capacity (Athletic): EM 1.8.

2.1.3 Forward dynamics

Forward dynamics takes the musculoskeletal model, initial conditions, muscle excitations, and external forces as the input and provides kinematic data over time, muscle states, and joint torques as the output. It computes muscle force using the Hill-type muscle model (Hill, 1938) based on current states and excitations, solves Equation 2 to get the accelerations, and integrates accelerations over time intervals Δt to update velocities (q˙t+1) and positions qt+1 as shown in Equation 5:

q˙t+1=q˙t+q¨Δt(5)
qt+1=qt+q˙t+1Δt

2.1.4 Joint reaction load

Joint reaction load analysis provides the forces and moments between two consecutive bodies based on the output of forward dynamics. This analysis determines the forces and moments acting on each cervical joint to maintain equilibrium, accounting for external forces, inertial forces, and internal forces, such as force generated from muscles and tendons. It is calculated by solving Newton-Euler equations where all translational and rotational dynamics are presented between those two consecutive bodies. The Newton-Euler Equations for estimating the atlantoaxial (C1-C2) joint reaction load is shown in Equation 6

RSL=M1q1q¨1Fm+Fgq1+Roccc1(6)

where M1q1 is the 6 × 6 mass matrix of cervical spine C1, q1 is the vector of linear and angular displacement of the C1, q¨1 is the vector of the linear and angular accelerations of the C1, Fm are the required muscle forces and moments to follow the given kinematics, Fgq1 is the gravitational loading, Roccc1 is the atlanto-occipital (O-C1) joint reaction force and moment, which can be found from the similar step for upper bodies, and RSL is the atlantoaxial (C1-C2) joint reaction force and moment. Rc and Rs are the compressive and shear forces on the atlantoaxial (C1-C2) joint in Figure 3.

Figure 3
Diagram of a human head and neck showing the brain and cervical spine labeled with anatomical structures: skull, cortex, corpus callosum, thalamus, brain stem, cerebellum, and cervical spine segments C1-C7. Enlarged insets highlight biomechanical forces and vectors on the cervical spine with labels: Occ, C1, C2, F_m, R_SL, F_g(q1), R_occ-c1, R_C, and R_S.

Figure 3. Joint reaction load analysis (simplified C1-C2 joint reaction calculation).

2.2 Integration of head, neck, and muscle injury criteria with the framework

2.2.1 Head injury criteria

Acceleration is measured in multiples of the acceleration of gravity (g), and time is measured in seconds. Head Injury Criteria (HIC) is proposed by the US National Highway Traffic Safety Administration (NHTSA), as shown in Equation 7

HIC=max1t2t1 t1t2atdt 2.5t2t1(7)

Where, t1 and t2 represent two arbitrary time points during the acceleration pulse. Acceleration is measured in gravitational force (g), and time is measured in seconds. HIC36 means that HIC is measured for each 36 ms apart (t2 and t1 is less than 36 ms apart). The HIC model is developed based on animal and cadaver experiments. Suggested human tolerance for 50th percentile males are: HIC36<1000 and HIC15<700. According to a study, HIC36<1000 represents 18% probability of a severe head injury, 55% probability of serious injury, and 90% probability of a moderate head injury to the average adult (Mackay, 2007). In another study, concussions were found at HIC15=250 in most athletes (Viano, 2005; Viano and Pellman, 2005).

However, rotational acceleration is responsible for shear stress that damages brain tissue. HIC does not consider rotational acceleration of the head in its injury evaluation. The Generalized Acceleration model Brain Injury Threshold (GAMBIT) model was developed by Newman (Newman, 1986) by combining translational and rotational components of head acceleration as shown in Equation 8.

GAMBIT=atacn+ϕ¨tϕc¨m1k(8)

Where at and ϕ¨t denote the translational and rotational acceleration, respectively. ac and ϕ¨c represent critical tolerance levels of those accelerations. It is reported that GAMBIT is one of the strongest predictors of human mTBI (Hernandez et al., 2015). Newman (1986) proposed 250 g and 10,000 rad/s2 as their critical thresholds. n,m, and k are constants (Zhan et al., 2021). A GAMBIT value of G = 1 represents a 50% probability of serious brain injury (Feist et al., 2009). Recent studies propose 80–100 g for linear accelerations, and 5,000-6,500 rad/s2 for 50% risk of mTBI (Funk et al., 2007; Rowson et al., 2012). However, to the best of authors’ knowledge, there is no contemporary studies that tried to connect the revised threshold based GAMBIT rating with Abbreviated Injury Scale (AIS) or probability of TBI risk. For this reason, in this study, we used the legacy thresholds 250 g and 10,000 rad/s2 for linear and rotational accelerations, respectively, as per the literature data we found (Chinn et al., 2001; Fernandes and Sousa, 2015; Schmitt et al., 2019). In this study, we assess the effect of impact locations. The eccentric impacts create significant rotational accelerations. As HIC does not consider rotational acceleration, we used GAMBIT to assess the head injury risk. However, as the GAMBIT threshold varies from study to study, we used HIC values to validate our data in comparison with the literature data.

Margulies and Thibault (1992) reported that the threshold value of shear strain has been reported to be about 5%–10% in case of DAI. It was observed that 50 percent concussion probability is observed within the corpus callosum at 21 percent strain (Kleiven, 2007). Deck and willinager proposed a 31% for maxmimum principal strain threshold and a 25% shear strain to indicate the onset of axonal damage. Rika carlsen adopted tissue level threshold 18% as the axonal strain injury tolerance. Bain et al. considered 18% strain as for the onset of electrophysiological impairment and 21% as morphological damage of axon. Kimapra et al. found out that at 63 g and 4,267 rad/s^2 accelerations, the brain tissue strain is about 29% for 25% mTBI. However, the macro level paramters and tissue level paraterms interactions requuiires further researche to devleop the bridge.

2.2.2 Neck injury criteria

The existing understanding of neck injury criteria is based on load combinations of axial force and flexion-extension moments that are measured with ATDs. Studies have shown that a compressive force of 1750 N–4800 N can cause cervical disc injury (Myers and Winkelstein, 1995; Whyte et al., 2020). One of the most utilized neck injury criteria is Nij (Prasad and Daniel, 1984) which is adopted by the US National Highway Traffic Safety Administration (NHTSA) for assessing automobile injuries. Nij is based on linear combination of axial force and bending moment. However, the neck injury criteria are based on force and moments, and the injuries are completely ignored due to coronal plane loads that cause lateral bending. As external impacts for this study are horizontal and not compressive forces, we considered relative velocity and acceleration based neck injury criteria (NIC) (Boström and Kullgren, 2007). The Neck Injury Criteria (NIC) is shown in Equation 9:

NIC t=H*arelt+vrelt2(9)

arelt,vrelt = Relative acceleration and velocity between the first thoracic vertebra (T1) and first cervical spine (C1), H is the neck length (Boström and Kullgren, 2007). Injury prediction is done considering the pressure gradient caused by a sudden change of the fluid flow inside the fluid compartment of the cervical spine. Bohman and Haland (2000) suggested NIC value 15 m2/s2 for AIS1 (minor injury). Nij is based on compressive force and moments on cervical spine and considers only flexion and extension related injury. Therefore, we considered NIC as we have lateral impacts, and it creates rotational movement along with flexion and extension.

2.2.3 Muscle-tendon strain

Muscle-tendon strain is expressed as a percentage increase of its current length compared to the length when it is at rest and developing no force. The resting length is called the slack length. The strain of muscle-tendon at any given instant can be expressed in Equation 10.

ϵMT=lMTlsMTlsMT(10)

Here, lMT is muscle-tendon current length, and lsMT is muscle-tendon slack length. In this study, we considered three hyoid muscles and three superficial multifidus muscles strain to assess the muscle injury probability. Three hyoid muscles: sternohyoid, sternothyroid, and omohyoid. The superficial multifidus muscles are superficial multifidus (C5/6-C2), superficial multifidus (C6/7-C2), and superficial multifidus (T1-C4). For more than 10% strain, tendon begins to experience mechanical failure, and there is a high risk of injury (Uchida and Delp, 2021).

2.3 Impact parameters

2.3.1 Impact severities (Γ)

We adopted the classification of impact severity from Beeman et al. (2012). Based on impact durations and force, we classified external forces into three categories: (i) Low-severity impact (LSI), Γ1: impact force <1500 N and impact durations >30ms, (ii) Medium severity impact (MSI), Γ2: (impact force 15003000 N and impact durations <30ms), (iii) High severity impact (HSI), Γ3: impact force >3500 N and impact durations <5ms. Average low velocity (<35 mph) automobile barrier impact tests conducted by the National Highway Traffic Safety Administration (NHTSA) were considered for LSI (National Highway Traffic Safety Administration (NHTSA)), as the impact duration is longer and the change of velocity and impact force are low during the test. American football head impacts data were used for the MSI as the impact forces are higher, the impact duration is shorter, and severity is moderate (Mortensen et al., 2020). Blast data were used as HSI as the force is very high, the impact duration is very short, and the severity is highest (Kulkarni et al., 2014). However, there is no standard classification of external forces. The impact severity is classified to test the hypotheses of this study.

2.3.2 Impact locations (θ)

The selection of impact locations was driven by both clinical injury prevalence and biomechanical considerations. These locations encompass a comprehensive range of impact directions observed in high-energy events such as automotive crashes, sports collisions, and blasts (Pellman et al., 2003; Post et al., 2015; National Highway Traffic Safety Administration, 2020). Post et al. (2015) reported that subdural hematoma occurred from impacts to the frontal and occipital regions, whereas impacts to the sides of the head produced parenchymal contusions than the impacts to the front and occipital regions of the head. Pellman et al. (2003) indicated that 78% of NFL impacts happen from the frontal left and right sides.

We categorized the location of head contact into seven impact locations (θ). Two anterior impacts, two posterior impacts, and three lateral impacts. We chose two impact locations for frontal and rear impacts, as head and neck anatomy are symmetric with respect to the sagittal plane. On the other hand, we choose three lateral impacts as head anatomical properties and neck muscles are asymmetric with respect to coronal plane. Figure 4a shows the impact locations. The two anterior impacts (θ1,θ2) are anterior central (AC) and anterior eccentric (AE). The two posterior impacts (θ3,θ4) are posterior central (PC) and posterior eccentric (PE). The lateral impacts (θ5,θ6, θ7) are lateral central (LC), anterolateral eccentric (ALE), and posterolateral eccentric (PLE). Seven impact locations cover key anatomical planes. We choose ALE and PLE lateral eccentric directions as it is reported that 57% of the NFL impacts involve these regions (Pellman et al., 2003).

Figure 4
Illustration showing (a) a schematic view of a head with labeled impact locations: AE, AC, PC, PE, PLE, LC, ALE. (b) Three diagrams depicting head accelerations with arrows: Translational (linear forward), Rotational (sideways turning), and Angular (tilting backward).

Figure 4. (a) Impact locations (AC = Anterior Concentric, AE = Anterior Eccentric, PC = Posterior Concentric, PE = Posterior Eccentric, LC = Lateral Concentric, PLE = Posterolateral Eccentric, ALE = Anterolateral Eccentric), (b) types of head acceleration.

There are three different types of acceleration that affect the skull and the head: translational, rotational, and angular acceleration, as shown in Figure 4b. For the translational type, also known as linear, the center of gravity of the head will move in a straight line without rotation, and this usually results in focal injury, while in rotational acceleration, there is no movement of the center, but the head will rotate around it, resulting in diffuse shear strain. Lastly, the angular type is a combination of both rotational and linear properties, and the center of gravity moves in an angular manner. The angular type is the most common type and is responsible for the bulk of the TBI cases.

2.4 Simulation workflow

2.4.1 Data collection

We collected three datasets to analyze the relationship of neck strength, impact locations, and impact severities on head and neck injury parameters. For the low-severity impacts (LSI), we collected an automotive crash test from the NHTSA database. The test number was 11,149 (National Highway Traffic Safety Administration, 2020). A frontal rigid barrier impact test was conducted on a 2020 Toyota Highlander SUV. The impact velocity of the vehicle was 56.35 km/h (35 mph). A 50th percentile male Hybrid III ATD was placed in the driver sitting position. The impact duration was about 75 ms, and the impact force was 1210N. For the moderate-severity impacts (MSI), we used extrapolated American Football data (Mortensen et al., 2018). The impact duration was 13 ms, and the maximum impact force was 1510N. For the high-severity impacts (HSI), we used simulated blast data (Kulkarni et al., 2014). The impact duration was 5 ms and the impact force was 3455N. The exact force and time duration data were also used to validate the model.

2.4.2 Musculoskeletal model simulation

A modified head-neck musculoskeletal model was employed in OpenSim software (Delp et al., 2007; Seth et al., 2018). The model was verified using 10 subjects’ data for moderate severity impacts. The head neck model is based on Mortensen (Mortensen et al., 2018) and Vasavada’s model (Vasavada et al., 1998). The model has a total of 72 Millard muscle actuators (Millard et al., 2013; Kuo et al., 2019). Neck muscles such as multifidus and Semispinalis Cervicis are short and highly pennated. The Millard muscles better capture muscle-specific behaviors, such as tendon compliance and fiber pennation angles. Moreover, the Millard model has detailed non-linear tendon dynamics. Therefore, it reflects neck muscle activation dynamics better than traditional hill-type muscle. The musculoskeletal model incorporates soft-tissue damping and rate-dependent ligaments. It also includes passive cervical spine ligaments and viscoelastic properties of muscles. These passive components stabilize the cervical spine during dynamic loading, ensuring biofidelic GAMBIT and NIC predictions for lateral and flexion-extension impacts. These properties are crucial for accurately simulating energy dissipation during rapid stretching in impacts. The details about the model can be found in (Kuo et al., 2018; Mortensen et al., 2018; Kuo et al., 2019; Mortensen et al., 2020).

Muscle activation data for neck stiffness were collected from the literature (Kuo et al., 2018; Kuo et al., 2019). Three neck strengths, three impact severities, and seven impact locations, a total of 63 cases, were performed. A MATLAB script was used to run these 63 cases in OpenSim forward dynamics simulations. We used OpenSim’s adaptive time step size for the forward dynamics. The maximum step size was set to 1, and the minimum step size was set to 1 × 10−7. The integrator tolerance was set to 1 × 10−5. During the forward dynamics, the integrator dynamically chooses a step size between these bounds to satisfy the integrator error tolerance. For slower and steadier motion, it will increase the step size to save computation time. For highly dynamic simulations, the step size shrinks to capture fast-changing dynamics. The workflow diagram is shown in Figure 5.

Figure 5
Flowchart showing the analysis of head and neck injuries using a musculoskeletal model. It includes sections for impact locations, neck strengths, and impact types. The chart outlines 63 cases and their forward dynamics, leading to calculations for head, cervical spine, and neck muscle injury criteria using specific equations. Sections are labeled with different subsections for detailed examination.

Figure 5. Musculoskeletal model simulation workflow for 63 cases.

2.4.3 Statistical analyses

Statistical analyses were performed using IBM SPSS Statistics version 29.0 (IBM Corp., Armonk, NY, USA). Initially, normality of the data was assessed with the Shapiro-Wilk test. Parametric one-way analysis of variance (ANOVA) was employed for normally distributed data, while the nonparametric Kruskal–Wallis test was utilized for datasets deviating from normality. For pairwise comparisons, the Mann-Whitney U test was conducted.

A post hoc power analysis for GAMBIT, conducted in SPSS using an effect size of f = 10.46 for the observed means (LSI: 0.6179, MSI: 0.8934, HSI: 0.9542; SDs: 0.0252, 0.0735, 0.0288; Table 3), indicates that n = 3 per neck strength category (N = 27 total) achieves a power of 1.000 for one-way ANOVA (α = 0.05, two-tailed). For a 20% difference in GAMBIT (f = 1.53), power is 1.000, confirming robustness for significant findings reported in Section 4.3.

3 Results

3.1 Head injury risk analysis based on linear and rotational accelerations

The linear accelerations are presented in Figure 6a. The rotational accelerations are presented in Figure 6b. The linear and rotational accelerations of the head are calculated at its center of mass.

Figure 6
Top image shows a diagram of a schematic view of a head with labeled impact locations: AE, AC, PC, PE, PLE, LC, ALE. Below are two sets of 3D bar graphs for three subjects (EM 1.2, EM 1.4, and EM 1.8). The left set represents linear accelerations measured in g, while the right set shows rotational accelerations in radians per second squared. Bars are color-coded from blue (low) to red (high).

Figure 6. Accelerations from 63 cases, (a) Linear accelerations and (b) rotational accelerations for three subjects, three impact severities, and seven impact locations.

3.2 Head injury risk analysis based on GAMBIT

The GAMBIT results, based on Equation 7 are presented in Figure 7.

Figure 7
Three 3D bar charts compare GAMBIT scores for three subjects with EM values 1.2, 1.4, and 1.8. Each chart includes vertical bars labeled with impact severity (HSI, MSI, LSI), and impact locations (AC, AE, PC, PE, LC, PLE, and ALE). Bars are color-coded from blue (low) to red (high).

Figure 7. GAMBIT for three subjects, three impact severities, and seven impact locations.

3.3 Neck injury risk assessment based on cervical spine force and moment

The calculated neck force and moment data are presented in Figure 8.

Figure 8
Six 3D bar charts display force and moment measurements for three subjects for different neck strengths (EM 1.2, 1.4, and 1.8). The left column shows force in newtons, while the right shows moment in newton meters. Bars are color-coded from blue (low) to red (high).

Figure 8. (a) neck force, (b) neck moment for three subjects, three impact severities, and seven impact locations.

3.4 Neck injury assessment based on NIC

The NIC data, based on Equation 8 are presented in Figure 9.

Figure 9
Three 3D bar charts compare NIC scores for three subjects with EM values 1.2, 1.4, and 1.8. Each chart includes vertical bars labeled with impact severity (HSI, MSI, LSI), and impact locations (AC, AE, PC, PE, LC, PLE, and ALE). Bars are color-coded from blue (low) to red (high).

Figure 9. NIC for three subjects, three impact severities, and seven impact locations.

3.5 Muscle strain injury

The strains for six different muscles, that take part in resistive action of the neck are presented in Figure 10.

Figure 10
(a) A 3D bar graph shows muscle strain percentages relative to impact locations, with a color gradient from blue (low strain) to red (high strain). (b) An anatomical illustration of neck muscles, including the omohyoid and sternohyoid, with labels. (c) A 3D anatomical model highlights the superficial multifidus muscles along the cervical spine.

Figure 10. (a) Muscle Strain for three subjects, three impact severities, and seven impact locations. (b) Hyoid muscles, (c) superficial multifidus muscles.

4 Discussion

In this study, we developed a computational musculoskeletal head-neck model using OpenSim software to investigate the effects of impact location, neck strength, and impact severity on head and neck injury parameters. We designed a total of 63 simulation cases, varying across seven impact locations, three neck strength levels, and three impact severities.

4.1 Effect of impact location

The box plots in Figure 11 illustrate the linear and rotational accelerations, GAMBIT, neck force, neck moment, and NIC across various impact locations, with data taken from three subjects. The boxes represent the interquartile range (IQR), the red horizontal lines indicate the medians, and the red circles show the mean values, while the whiskers display the minimum and maximum values within the data range. Table 1 shows the mean value of each data for different locations. It can be seen from 9 (a) that linear accelerations do not vary much based on impact directions and locations and stay within 194 g–218 g. Eccentric impacts create higher linear acceleration than central impacts however, the increment is less than 10%. The lateral impacts generate about 10% more linear accelerations than anterior and posterior impacts. The hyoid muscles are highly stretched (>10%) during frontal and ALE impacts, whereas the multifidus muscles were highly stretched during posterior and lateral impacts.

Figure 11
Six box plots, labeled (a) to (f), show data across various impact locations: AC, AE, PC, PE, LC, PLE, and ALE. (a) shows linear acceleration, (b) rotational acceleration, (c) GAMBIT values, (d) force, (e) moment, and (f) NIC. Each plot includes median lines, whiskers, and overlaid red line graphs indicating trends.

Figure 11. Effect of impact locations on (a) linear acceleration, (b) rotational acceleration, (c) GAMBIT, (d) neck force, (e) neck moment, and (f) NIC.

Table 1
www.frontiersin.org

Table 1. Effect of impact locations.

To quantify the effect of impact location on injury parameters, we conducted Kruskal–Wallis tests due to normality violations (Shapiro-Wilk p < 0.05 for most variables). Significant main effects were found for NIC (χ26 = 16.650, p = 0.011), rotational acceleration (χ26 = 49.357, p < 0.001), neck force (χ26 = 45.771, p < 0.001), and neck moment (χ26 = 48.720, p < 0.001), but not for GAMBIT (p = 0.266) or linear acceleration (p = 0.289). Pairwise Mann-Whitney U tests (Bonferroni-corrected, p < 0.025) comparing anterolateral eccentric (ALE) to anterior central (AC) and anterior eccentric (AE) impacts confirmed significant differences for NIC (ALE: 21.28 m2/s2 vs. AC: 15.64 m2/s2, p = 0.004; vs. AE: 15.80 m2/s2, p = 0.002), rotational acceleration (ALE: 4,554.66 rad/s2 vs. AC: 893.20 rad/s2, p < 0.001; vs. AE: 1,148.54 rad/s2, p < 0.001), neck force (ALE: 201.10 N vs. AC: 602.64 N, p < 0.001; vs. AE: 638.88 N, p < 0.001), and neck moment (ALE: 7.27 Nm vs. AC: 1.57 Nm, p < 0.001; vs. AE: 1.98 Nm, p < 0.001). These results confirm that ALE impacts produce significantly higher rotational acceleration and neck moment, supporting the claim of increased injury risk for anterolateral impacts. No significant differences were found for GAMBIT (ALE: 0.98 vs. AC: 0.88, p = 0.046; vs. AE: 0.88, p = 0.050) or linear acceleration (ALE: 241.98 m2/s2 vs. AC: 218.52 m2/s2, p = 0.059; vs. AE: 219.92 m2/s2, p = 0.063). The non-significant results for GAMBIT suggest that GAMBIT is more sensitive to the magnitude of acceleration (severity) than its directional application (location) (1.44x and 1.54x higher for MSI and HSI vs. LSI, as mentioned in Section 4.3). Also, it indicates that head injuries are more sensitive to rotational accelerations than linear accelerattions which aligns with previous studies (Wright, 2012; Wright et al., 2013). Anterolateral eccentric (ALE) impacts exhibited greater consistency in neck force values, suggesting a stable biomechanical response that may contribute to their classification as the most risky impact location compared to other locations.

The rotational accelerations in Figure 11b vary significantly based on impact directions and locations. The eccentric impacts generated about 1.51, 2.61, and 4.11 times more rotational acceleration than central impacts, for the anterior, posterior, and lateral sides, respectively. Posterior and lateral impacts create 1.58 and 2.51 times more rotational acceleration than anterior impacts, respectively. The maximum average rotation acceleration is for anterolateral eccentric impact (4,176 rad/s2) which is 4.75 times more than the average anterior central impact (879 rad/s2). The rotational acceleration of 4,267 rad/s2 has 25% chance of mTBI (King, 2018; Zaman et al., 2024). Translational or linear accelerations cause the focal injuries to the brain, such as contusion and hematoma, whereas rotational acceleration causes diffuse injuries, such as DAI and subdural hematoma. Rotational accelerations create multidirectional strain fields that primarily affect the white matter structures such as axons. Margulies and Thibault (Margulies and Thibault, 1992) reported that the threshold for shear strain associated with diffuse axonal injury (DAI) ranges from 5% to 10%. Kleiven (2007) mentioned that a 50% concussion probability occurs within the corpus callosum at a 21% strain level. Deck and Willinger (2008) proposed a maximum principal strain threshold of 31% and a shear strain threshold of 25% to indicate the onset of axonal damage. Write et al. (2013) adopted an axonal strain injury tolerance at 18%, while Bain and Meaney (2000) identified 18% strain as the threshold for electrophysiological impairment and 21% for morphological axonal damage. Kimpara and Iwamoto (2012) determined that brain tissue strain reaches approximately 29% under accelerations of 63g and 4,267 rad/s2, corresponding to a 25% mTBI risk. However, further research is needed to clarify the interplay between macro-level parameters and tissue-level metrics to establish a comprehensive link between these scales.

The GAMBIT values in Figure 11c do not vary very much with the impact directions and locations similar to linear accelerations. It was within 0.79–0.83 for all impact locations except for anterolateral eccentric, which was 0.89. In summary, impact locations have a significant effect on rotational acceleration, and anterolateral impact is the riskiest impact location for both head and neck.

The average neck forces for the anterior impacts were about 2.5 times higher than for posterior and lateral impacts. The maximum average neck force was for the anterior eccentric impacts (676N). Neck moments in Figure 11e vary significantly based on impact locations and directions. The highest average neck moment was for anterolateral eccentric (9.43Nm), which was 5.85 times more than the posterior central neck moment (1.61Nm). Also, eccentric impacts generate about 157% more neck moments than central impacts. NIC values in Figure 11f were within 15–24 m2/s2. The lateral impact values were higher than the anterior and posterior impacts.

4.2 Effect of neck strength

The effects of neck strength are presented in Figure 12 and Table 2. Average linear accelerations in Figure 12a do not vary much with neck muscle strength and were within 202 g to 205 g. The average rotational accelerations in Figure 12b also do not vary much and stay between 1960 rad/s2 to 2037 rad/s2. The average GAMBIT in Figure 12c ranges from 0.81-0.83.

Figure 12
Six box plots (a-f) depict the effects of neck strength (EM 1.2, 1.4, 1.8) on various metrics: (a) Linear Acceleration (g), (b) Rotational Acceleration (rad/s²), (c) GAMBIT, (d) Force (N), (e) Moment (Nm), and (f) NIC. Each plot shows medians and variability with red circles indicating means.

Figure 12. Effect of neck strength on (a) linear acceleration, (b) rotational acceleration, (c) GAMBIT, (d) neck force, (e) neck moment, and (f) NIC.

Table 2
www.frontiersin.org

Table 2. Effect of neck strength.

Kruskal–Wallis tests showed no significant effects of neck strength (EM 1.2, 1.4, 1.8) on any injury parameters (p > 0.184 for linear acceleration, rotational acceleration, GAMBIT, neck force, neck moment, and NIC), consistent with the minimal variation (< 10%) in medians and means across neck strength levels.

The average neck force in Figure 12d ranges from 334N to 417N. The average neck moment in Figure 12e ranges from 4.5Nm to 5.6Nm. The average NIC value ranges from 14.7 to 15.5. In summary, the average head and neck injury parameters do not vary by more than 10% based on neck strength, which agrees with the literature (Mortensen, Vasavada et al., and statistical analysis.

4.3 Effect of impact severities

The effects of impact severities are presented in Figure 13 and Table 3. Kruskal–Wallis tests revealed significant effects of impact severities (LSI, MSI, HSI) on NIC (p < 0.001), GAMBIT (p < 0.001), and linear acceleration (p < 0.001), but not on rotational acceleration (p = 0.183), neck force (p = 0.064), or neck moment (p = 0.048). Pairwise Mann-Whitney U tests (Bonferroni-corrected, p < 0.0167) showed HSI significantly exceeded LSI for NIC (19.69 vs. 11.82 m2/s2, p < 0.001), GAMBIT (0.95 vs. 0.61, p < 0.001), and linear acceleration (236.60 vs. 153.09 m/s2, p < 0.001), and exceeded MSI for NIC (19.69 vs. 16.21 m2/s2, p < 0.001), GAMBIT (0.95 vs. 0.90, p < 0.001), and linear acceleration (236.60 vs. 223.65 m/s2, p < 0.001). MSI also exceeded LSI for GAMBIT (0.90 vs. 0.61, p < 0.001) and linear acceleration (223.65 vs. 153.09 m/s2, p < 0.001), but not for NIC (p = 0.030, non-significant at 0.0167).

Figure 13
Six box plots compare different impact severity measures: (a) Linear Acceleration, (b) Rotational Acceleration, (c) GAMBIT, (d) Force, (e) Moment, and (f) NIC across three severity levels: LSI, ISI, and HSI. Each plot shows data distribution, median, and trend lines indicating changes with increasing severity, highlighting higher values in HSI compared to LSI and ISI for most metrics.

Figure 13. Effect of impact severities on (a) linear acceleration, (b) rotational acceleration, (c) GAMBIT, (d) neck force, (e) neck moment, and (f) NIC.

Table 3
www.frontiersin.org

Table 3. Effect of impact severities.

Our model incorporates both passive and active muscle responses, as detailed in Section 2.1.2 (Equation 3), using muscle activation data from literature (Kuo et al., 2018; Kuo et al., 2019). However, in HSI, the rapid onset limits active muscle activation, making passive responses dominant. For LSI and MSI impacts, active muscle responses contribute significantly, reflecting anticipatory stiffening in automotive or sports scenarios. NIC and GAMBIT may be conservatively high due to limited active response time for HSI. However, we ignored pre-activation in all 63 cases, as it is reported that pre-activation can reduce head acceleration by ∼30% in sports contexts (Fanta et al., 2013). This was omitted to standardize responses across all impact severities. It ensures consistency across the 63 cases and avoids bias from sports-specific pre-activation data.

The average linear accelerations in Figure 13a for MSI and HSI are 1.45 and 1.54 times higher than LSI, and range from 153 g to 237 g. The maximum linear accelerations were 251 g for HSI. The average rotational accelerations in Figure 13b for MSI and HSI are respectively 1.18 and 1.35 times higher than LSI, and range from 1704 rad/s2 to 2,290 rad/s2. The maximum rotational accelerations were 4859 rad/s2 for HSI. The average GAMBIT in Figure 13c for MSI and HSI were 1.44 and 1.54 times higher than LSI. The maximum GAMBIT was 1.02 for HSI.

The neck forces and moments did not vary much based on the impact severities, as in Figures 13d,e. The NIC in Figure 13f for MSI and HSI were 1.24 and 1.68 times more than LSI. The maximum NIC value was 30.28 for HSI. In summary, the risk probability of head and neck injuries is higher for HSI and MSI than LSI.

4.4 Validation

The musculoskeletal model was validated in previous studies for MSI using NFL (Mortensen et al., 2018; Kuo et al., 2019; Mortensen et al., 2020) and rugby (Silvestros et al., 2024) data. However, it was not validated for LSI and HSI. To validate the model for LSI and HSI, we ran the simulation using the literature data for LSI and HSI for all neck strength. The literature only has the HIC value and max resultant linear acceleration, which we compared with our result for validation. However, the neck strengths of the experimental subjects were not mentioned in the literature. So, we ran the simulation for all neck strength. We can see from Table 4, for LSI, the experimental data indicates that the subject’s (ATD’s) neck strength was between EM 1.2 and EM 1.4. For HSI, the experimental data indicates that the subject’s (ATD’s) neck strength was between 1.4 and 1.8. The simulation results were consistent with the literature data.

Table 4
www.frontiersin.org

Table 4. Validation for LSI and HSI data.

We also compared our average results with Hybrid III ATD results (Viano, 2005) as shown in Table 5. Although the experiment design and impact forces differ, the means of most results are within one standard deviation for rotational accelerations, neck force, and moments.

Table 5
www.frontiersin.org

Table 5. Validation of model data with ATD experiments.

From our results, the rotational accelerations of the anterolateral eccentric and posterolateral eccentric are 4.7 and 3.7 times higher than the anterior central impacts, respectively. It is mentioned in the literature that spin rotational acceleration (anterolateral eccentric or posterolateral eccentric) of the head is 3.5 times higher than whiplash acceleration (anterior impact) (Razzaghi, 2019). Our acceleration data is consistent with the literature. The distance ratio from the front of the head to the center of mass of the head and from back of the skull to the center of mass of the skull is 0.52 and 0.48. That means the anterolateral impact will create more acceleration and moment than the posterolateral impact. Viano et al. (2007) studied twenty-five helmet collisions involving NFL players and found that the maximum rotation occurred around the superior-inferior axis. The maximum acceleration was 9,678 rad/s2 and happened on the anterolateral locations. This finding is consistent with our results.

4.5 Limitation

There are some limitations to this study. First, we ignored torso velocity and accelerations on the impact of head and neck injuries. We considered the torso as a fixed body. Secondly, we did not consider the whole blast impact on the head and neck injuries. We considered an idealized blast scenario, which did not account for multidimensional waveforms. This does not perfectly reproduce the experimental environment and head impacts. As a result, it may underestimate the rotational kinematics. Future models should incorporate these complexities to improve the realism. Thirdly, we consider height, weight, and anthropometric data for all tests to be the same, in order to consider the effects of neck strength, impact locations, and impact severities on head and neck injuries. However, with different anthropometric and muscle activation data, head-neck responses after an incident may vary. Also, we used the MSI muscle activation data for all cases. The muscles’ activations data were collected from (Mortensen et al., 2018; Mortensen et al., 2020).

Fourth, while n = 3 per neck strength category provides high power (1.000 for observed GAMBIT differences, 1.000 for 20% differences, Section 2.4.3), it increases the risk of Type II errors for smaller effect sizes, particularly for non-significant parameters like rotational acceleration (p = 0.183) and neck moment (p = 0.048, Section 2.4). Fifth, we considered that the left and right sides of the human head and brain are symmetric. The brain has different lobes and different functional areas, as shown in Figure 14a. Also, the internal structure of brain tissue is heterogeneous and consists of grey matter, white matter, cerebrospinal fluid, and blood vessels as shown in Figure 14b. The outer layer of the brain (skull) can be considered as symmetric; however, the brain is not symmetric. Also, the functions of left-side brain and right-side brain are different. The right hemisphere is responsible for analytical thinking and logical expressions, whereas the left hemisphere is responsible for creativity and intuition. The symptoms caused by the left side impacts and right side impacts may be different. Sixth, the subjects in this study were scaled as an adult. The anthropometric data of infants, children, and adults are different. Infants are not miniature adults. Infants and children have greater head-mass to body-mass ratio. Also, their head’s center of gravity with respect to the neck is higher than that of adults (Burdi et al., 1969). All these factors may increase their TBI risk compared to adults for similar incidents (Eckner et al., 2014). The results on a scaled musculoskeletal model based on children may be different than this study and a scope for future researchers.

Figure 14
Illustration of the human brain with labeled regions. Left image shows the side view highlighting frontal, parietal, temporal, and occipital lobes, each associated with specific functions like motor, sensory, auditory, and visual. Right image shows a cross-section displaying internal structures including grey matter, white matter, thalamus, and hippocampus.

Figure 14. (a) Lobes and functional areas of the brain, (b) cross-section of the brain. Images modified and used with permission from Servier Medical Art, Creative Commons Attribution 4.0 Unreported License.

Also, the biomechanical model utilized in this study is based on a 50th percentile male, which may limit its generalizability to female populations. Sex-based differences in neck musculature and TBI susceptibility suggest that females may exhibit distinct injury responses due to variations in anatomy and physiology (Putra et al., 2022). This study does not quantify the differences in inter-subject variability, such as age, sex, and muscle composition. Moreover, in this study, we did not quantify the spatial distribution of impact eccentricity, such as moment arm distances. This limitation may restrict a detailed dose-response analysis that future research could explore, if relevant relationships are identified. Careful consideration should be given to the limitations outlined in this study during interpretation of the presented data and decision making. This study provides valuable insights. However, the limitations of those study should be kept in mind for any conclusion in a broader context. Further works are necessary to enhance the robustness and generalizability of the outcomes.

5 Conclusion

In this study, we investigated the varying relationship of neck strength, impact locations, and impact severities with head and neck injury criteria. We employed a musculoskeletal model and used forward dynamics to run 63 cases. The results were validated using different literature data. It was found that impact locations have a significant effect on head and neck injury parameters, and anterolateral eccentric impact is the riskiest impact location for both head and neck. We also found that average head and neck injury parameters do not vary more than 10% based on neck strength. Finally, the risk probability of head and neck injuries is higher for HSI and MSI than for LSI. These findings indicate that people who face an anterolateral eccentric impact are at higher risk than the person who faces an impact from other directions and should seek medical attention. The biomechanical findings of this study suggest a need for enhanced lateral impact protection in helmet design, potentially through energy-absorbing materials. These recommendations are preliminary and require clinical validation to ensure efficacy. To our knowledge, this is the first musculoskeletal model based study, where head and neck injury parameters were assessed for multiple impact locations and for all three types of severity impacts. Although the musculoskeletal-model-based framework has limitations, it provides vital information about head and neck injuries for future researchers. The high biofidelity data from this study can be utilized as the input for tissue-level and molecular-level head and neck injury studies. That will enhance our understanding about the cause-and-effects of upper extremity organ-level injuries.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

Ethical approval was not required for the study involving humans in accordance with the local legislation and institutional requirements. Written informed consent to participate in this study was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and the institutional requirements.

Author contributions

RZ: Data curation, Formal Analysis, Investigation, Methodology, Writing – original draft. MH: Data curation, Formal Analysis, Writing – original draft. AR: Formal Analysis, Visualization, Writing – original draft. SC: Formal Analysis, Visualization, Writing – original draft. AJ: Writing – review and editing. AK: Writing – review and editing. AA: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing – review and editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. AA acknowledges the financial support from Office of Naval Research (ONR)’s Force Health Protection (FHP) program (through the Award # ONR: N00014-21-1-2051 and ONR:N000142412324 Dr. Timothy Bentley, Program Manager).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Bain, A. C., and Meaney, D. F. (2000). Tissue-level thresholds for axonal damage in an experimental model of central nervous system white matter injury. J. Biomech. Eng. 122 (6), 615–622. doi:10.1115/1.1324667

PubMed Abstract | CrossRef Full Text | Google Scholar

Beeman, S. M., Kemper, A. R., Madigan, M. L., Franck, C. T., and Loftus, S. C. (2012). Occupant kinematics in low-speed frontal sled tests: human volunteers, hybrid III ATD, and PMHS. Accid. Analysis and Prev. 47, 128–139. doi:10.1016/j.aap.2012.01.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Bohman, K., Haland, Y., H Land, Y., and Kullgren, A. (2000). A study of AIS1 neck injury parameters in 168 frontal collisions using a restrained hybrid III dummy. SAE Tech. Pap. 44, 103–116. doi:10.4271/2000-01-sc08

CrossRef Full Text | Google Scholar

Borich, M. R., Cheung, K. L., Jones, P., Khramova, V., Gavrailoff, L., Boyd, L. A., et al. (2013). Concussion: current concepts in diagnosis and management. J. Neurologic Phys. Ther. 37 (3), 133–139. doi:10.1097/npt.0b013e31829f7460

PubMed Abstract | CrossRef Full Text | Google Scholar

Boström, O., and Kullgren, A. (2007). Characteristics of anti-whiplash seat designs with good real-life performance. Proc. IRCOBI Conf.

Google Scholar

Browne, K. D., Chen, X.-H., Meaney, D. F., and Smith, D. H. (2011). Mild traumatic brain injury and diffuse axonal injury in swine. J. neurotrauma 28 (9), 1747–1755. doi:10.1089/neu.2011.1913

PubMed Abstract | CrossRef Full Text | Google Scholar

Bumberger, R., Acar, M., and Bouazza-Marouf, K. (2025). Parametric computational head-and-neck model and anthropometric effects on whiplash loading.

Google Scholar

Burdi, A. R., Huelke, D. F., Snyder, R. G., and Lowrey, G. H. (1969). Infants and children in the adult world of automobile safety design: pediatric and anatomical considerations for design of child restraints. J. Biomechanics 2 (3), 267–280. doi:10.1016/0021-9290(69)90083-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Centers for Disease Control and Prevention (CDC) (2019). TBI surveillance report US 2018 2019.

Google Scholar

Chinn, B., Canaple, B., Derler, S., Doyle, D., Otte, D., Schuller, E., et al. (2001). COST 327 motorcycle safety helmets, 3. Belgium: European Commission, Directorate General for Energy and Transport, 1135–1148.

Google Scholar

Deck, C., and Willinger, R. (2008). Improved head injury criteria based on head FE model. Int. J. Crashworthiness 13 (6), 667–678. doi:10.1080/13588260802411523

CrossRef Full Text | Google Scholar

Defense Health Agency (DHA) (2024). DOD numbers for traumatic brain injury worldwide. Def. Health Agency (DHA).

Google Scholar

Delp, S. L., Anderson, F. C., Arnold, A. S., Loan, P., Habib, A., John, C. T., et al. (2007). OpenSim: open-source software to create and analyze dynamic simulations of movement. IEEE Trans. Biomed. Eng. 54 (11), 1940–1950. doi:10.1109/tbme.2007.901024

PubMed Abstract | CrossRef Full Text | Google Scholar

Deng, Y.-C., and Goldsmith, W. (1987). Response of a human head/neck/upper-torso replica to dynamic loading—II. analytical/Numerical model. J. biomechanics 20 (5), 487–497. doi:10.1016/0021-9290(87)90249-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Dymek, M., Ptak, M., Ratajczak, M., Fernandes, F. A., Kwiatkowski, A., and Wilhelm, J. (2021). Analysis of HIC and hydrostatic pressure in the human head during NOCSAE tests of American football helmets. Brain Sci. 11 (3), 287. doi:10.3390/brainsci11030287

PubMed Abstract | CrossRef Full Text | Google Scholar

Eckner, J. T., Oh, Y. K., Joshi, M. S., Richardson, J. K., and Ashton-Miller, J. A. (2014). Effect of neck muscle strength and anticipatory cervical muscle activation on the kinematic response of the head to impulsive loads. Am. J. sports Med. 42 (3), 566–576. doi:10.1177/0363546513517869

PubMed Abstract | CrossRef Full Text | Google Scholar

Escolas, S. M., Luton, M., Ferdosi, H., Chavez, B. D., and Engel, S. D. (2020). Traumatic brain injuries: unreported and untreated in an army population. Mil. Med. 185 (Suppl. ment_1), 154–160. doi:10.1093/milmed/usz259

PubMed Abstract | CrossRef Full Text | Google Scholar

Fanta, O., Hadraba, D., Lopot, F., Kubový, P., Bouček, J., and Jelen, K. (2013). Pre-activation and muscle activity during frontal impact in relation to whiplash associated disorders. Neuroendocrinol. Lett. 34 (7), 708–716.

PubMed Abstract | Google Scholar

Farmer, J., Mitchell, S., Sherratt, P., and Miyazaki, Y. (2022). A human surrogate neck for traumatic brain injury research. Front. Bioeng. Biotechnol. 10, 854405. doi:10.3389/fbioe.2022.854405

PubMed Abstract | CrossRef Full Text | Google Scholar

Feist, F., Gugler, J., Arregui-Dalmases, C., Del Pozo de Dios, E., López-Valdés, F., Deck, D., et al. (2009). “Pedestrian collisions with flat-fronted vehicles: injury patterns and importance of rotational accelerations as a predictor for traumatic brain injury (TBI),” in 21st international technical conference on the enhanced safety of vehicles (ESV).

Google Scholar

Fernandes, F. A., and Sousa, R. J. A. d. (2015). Head injury predictors in sports trauma–a state-of-the-art review. Proc. Institution Mech. Eng. Part H J. Eng. Med. 229 (8), 592–608. doi:10.1177/0954411915592906

CrossRef Full Text | Google Scholar

Funk, J. R., Duma, S., Manoogian, S., and Rowson, S. (2007). “Biomechanical risk estimates for mild traumatic brain injury,” in Annual proceedings/association for the advancement of automotive medicine (Melbourne, VIC: Association for the Advancement of Automotive Medicine).

Google Scholar

Gabler, L. F., Crandall, J. R., and Panzer, M. B. (2019). Development of a second-order system for rapid estimation of maximum brain strain. Ann. Biomed. Eng. 47, 1971–1981. doi:10.1007/s10439-018-02179-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Garcia, T., and Ravani, B. (2003). A biomechanical evaluation of whiplash using a multi-body dynamic model. J. Biomech. Eng. 125 (2), 254–265. doi:10.1115/1.1556856

PubMed Abstract | CrossRef Full Text | Google Scholar

Halloran, J. P., Ackermann, M., Erdemir, A., and Van den Bogert, A. J. (2010). Concurrent musculoskeletal dynamics and finite element analysis predicts altered gait patterns to reduce foot tissue loading. J. biomechanics 43 (14), 2810–2815. doi:10.1016/j.jbiomech.2010.05.036

PubMed Abstract | CrossRef Full Text | Google Scholar

Hasegawa, K., Takeda, T., Nakajima, K., Ozawa, T., Ishigami, K., Narimatsu, K., et al. (2014). Does clenching reduce indirect head acceleration during rugby contact? Dent. Traumatol. 30 (4), 259–264. doi:10.1111/edt.12082

PubMed Abstract | CrossRef Full Text | Google Scholar

Hernandez, F., Shull, P. B., and Camarillo, D. B. (2015). Evaluation of a laboratory model of human head impact biomechanics. J. biomechanics 48 (12), 3469–3477. doi:10.1016/j.jbiomech.2015.05.034

PubMed Abstract | CrossRef Full Text | Google Scholar

Hill, A. V. (1938). The heat of shortening and the dynamic constants of muscle. Proc. R. Soc. Lond. Ser. B-Biological Sci. 126 (843), 136–195.

Google Scholar

Ho, J., Zhou, Z., Li, X., and Kleiven, S. (2017). The peculiar properties of the falx and tentorium in brain injury biomechanics. J. biomechanics 60, 243–247. doi:10.1016/j.jbiomech.2017.06.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Huo, F., and Wu, H. (2017). Predicting the head-neck posture and muscle force of the driver based on the combination of biomechanics with multibody dynamics. SAE Tech. Pap. doi:10.4271/2017-01-0407

CrossRef Full Text | Google Scholar

Insurance Institute for Highway Safety and Highway Loss Data Institute (IIHS-HLDI) (2024). Fatality facts 2022 yearly snapshot.

Google Scholar

Jin, X., Feng, Z., Mika, V., Li, H., Viano, D. C., and Yang, K. H. (2017). The role of neck muscle activities on the risk of mild traumatic brain injury in American football. J. biomechanical Eng. 139 (10), 101002. doi:10.1115/1.4037399

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnson, D., Koya, B., and Gayzik, F. S. (2020). Comparison of neck injury criteria values across human body models of varying complexity. Front. Bioeng. Biotechnol. 8, 985. doi:10.3389/fbioe.2020.00985

PubMed Abstract | CrossRef Full Text | Google Scholar

Kimpara, H., and Iwamoto, M. (2012). Mild traumatic brain injury predictors based on angular accelerations during impacts. Ann. Biomed. Eng. 40 (1), 114–126. doi:10.1007/s10439-011-0414-2

PubMed Abstract | CrossRef Full Text | Google Scholar

King, A. I. (2018). The biomechanics of impact injury.

Google Scholar

Kleiven, S. (2007). Predictors for traumatic brain injuries evaluated through accident reconstructions. SAE Tech. Pap. 51, 81–114. doi:10.4271/2007-22-0003

PubMed Abstract | CrossRef Full Text | Google Scholar

Kulkarni, K. B., Ramalingam, J., and Thyagarajan, R. (2014). Assessment of the accuracy of certain reduced order models used in the prediction of occupant injury during under-body blast events. SAE Int. J. Transp. Saf. 2 (2), 307–319. doi:10.4271/2014-01-0752

CrossRef Full Text | Google Scholar

Kuo, C., Fanton, M., Wu, L., and Camarillo, D. (2018). Spinal constraint modulates head instantaneous center of rotation and dictates head angular motion. J. biomechanics 76, 220–228. doi:10.1016/j.jbiomech.2018.05.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuo, C., Sheffels, J., Fanton, M., Yu, I. B., Hamalainen, R., and Camarillo, D. (2019). Passive cervical spine ligaments provide stability during head impacts. J. R. Soc. Interface 16 (154), 20190086. doi:10.1098/rsif.2019.0086

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J., Jin, J. J., Eckner, J. T., Ji, S., and Hu, J. (2022). Influence of morphological variation on brain impact responses among youth and young adults. J. biomechanics 135, 111036. doi:10.1016/j.jbiomech.2022.111036

PubMed Abstract | CrossRef Full Text | Google Scholar

Mackay, M. (2007). The increasing importance of the biomechanics of impact trauma. Sadhana 32 (4), 397–408. doi:10.1007/s12046-007-0031-9

CrossRef Full Text | Google Scholar

Mang, D. W., Siegmund, G. P., Brown, H. J., Goonetilleke, S. C., and Blouin, J.-S. (2015). Loud preimpact tones reduce the cervical multifidus muscle response during rear-end collisions: a potential method for reducing whiplash injuries. spine J. 15 (1), 153–161. doi:10.1016/j.spinee.2014.08.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Margulies, S. S., and Thibault, L. E. (1992). A proposed tolerance criterion for diffuse axonal injury in man. J. biomechanics 25 (8), 917–923. doi:10.1016/0021-9290(92)90231-o

PubMed Abstract | CrossRef Full Text | Google Scholar

McCrea, M., Hammeke, T., Olsen, G., Leo, P., and Guskiewicz, K. (2004). Unreported concussion in high school football players: implications for prevention. Clin. J. sport Med. 14 (1), 13–17. doi:10.1097/00042752-200401000-00003

PubMed Abstract | CrossRef Full Text | Google Scholar

Merrill, T., Goldsmith, W., and Deng, Y. (1984). Three-dimensional response of a lumped parameter head-neck model due to impact and impulsive loading. J. biomechanics 17 (2), 81–95. doi:10.1016/0021-9290(84)90126-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Mihalik, J. P., Guskiewicz, K. M., Marshall, S. W., Greenwald, R. M., Blackburn, J. T., and Cantu, R. C. (2011). Does cervical muscle strength in youth ice hockey players affect head impact biomechanics? Clin. J. sport Med. 21 (5), 416–421. doi:10.1097/jsm.0b013e31822c8a5c

PubMed Abstract | CrossRef Full Text | Google Scholar

Millard, M., Uchida, T., Seth, A., and Delp, S. L. (2013). Flexing computational muscle: modeling and simulation of musculotendon dynamics. J. Biomechanical Eng. 135 (2), 021005. doi:10.1115/1.4023390

PubMed Abstract | CrossRef Full Text | Google Scholar

Mortensen, J. D., Vasavada, A. N., and Merryweather, A. S. (2018). The inclusion of hyoid muscles improve moment generating capacity and dynamic simulations in musculoskeletal models of the head and neck. PloS one 13 (6), e0199912. doi:10.1371/journal.pone.0199912

PubMed Abstract | CrossRef Full Text | Google Scholar

Mortensen, J. D., Vasavada, A. N., and Merryweather, A. S. (2020). Sensitivity analysis of muscle properties and impact parameters on head injury risk in American football. J. biomechanics 100, 109411. doi:10.1016/j.jbiomech.2019.109411

PubMed Abstract | CrossRef Full Text | Google Scholar

Mychasiuk, R., Hehar, H., Candy, S., Ma, I., and Esser, M. J. (2016). The direction of the acceleration and rotational forces associated with mild traumatic brain injury in rodents effect behavioural and molecular outcomes. J. Neurosci. methods 257, 168–178. doi:10.1016/j.jneumeth.2015.10.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Myers, B. S., and Winkelstein, B. A. (1995). Epidemiology, classification, mechanism, and tolerance of human cervical spine injuries. Crit. Reviews™ Biomed. Eng. 23 (5-6), 307–409. doi:10.1615/critrevbiomedeng.v23.i5-6.10

PubMed Abstract | CrossRef Full Text | Google Scholar

National Highway Traffic Safety Administration (2020). “New car assessment program (NCAP), frontal barrier impact test, Toyota Motor manufacturing, Indiana, INC, 2020 Toyota highlander,” in Five door suv, nhtsa no: m20205100, u. s. department of transportation, national highway traffic safety administration, office of crashworthiness standards.

Google Scholar

Newman, J. (1986). “A generalized acceleration model for brain injury threshold (GAMBIT),” in Proceedings of international IRCOBI conference, 1986.

Google Scholar

Pellman, E. J., Viano, D. C., Tucker, A. M., and Casson, I. R. (2003). Concussion in professional football: location and direction of helmet impacts—Part 2. Neurosurgery 53 (6), 1328–1341. doi:10.1227/01.neu.0000093499.20604.21

PubMed Abstract | CrossRef Full Text | Google Scholar

Post, A., Hoshizaki, T. B., Gilchrist, M. D., Brien, S., Cusimano, M. D., and Marshall, S. (2014). The influence of acceleration loading curve characteristics on traumatic brain injury. J. biomechanics 47 (5), 1074–1081. doi:10.1016/j.jbiomech.2013.12.026

CrossRef Full Text | Google Scholar

Post, A., Hoshizaki, T. B., Gilchrist, M. D., Brien, S., Cusimano, M., and Marshall, S. (2015). Traumatic brain injuries: the influence of the direction of impact. Neurosurgery 76 (1), 81–91. doi:10.1227/neu.0000000000000554

PubMed Abstract | CrossRef Full Text | Google Scholar

Pramudita, J. A., Ujihashi, S., Ono, K., Ejima, S., Sato, F., Yamazaki, K., et al. (2009). Development of a head-neck finite element model and analysis of intervertebral strain response during rear impact. Trans. Jpn. Soc. Mech. Eng. Ser. A 75 (759), 1549–1555. doi:10.1299/kikaia.75.1549

CrossRef Full Text | Google Scholar

Prasad, P., and Daniel, R. P. (1984). A biomechanical analysis of head, neck, and torso injuries to child surrogates due to sudden torso acceleration. SAE Trans., 784–799.

Google Scholar

Putra, I. P. A., Iraeus, J., Sato, F., Svensson, M. Y., and Thomson, R. (2022). Finite element human body models with active reflexive muscles suitable for sex based whiplash injury prediction. Front. Bioeng. Biotechnol. 10, 968939. doi:10.3389/fbioe.2022.968939

PubMed Abstract | CrossRef Full Text | Google Scholar

Razzaghi, M. (2019). Mechanical analysis of head impact. Open Access J. Neurology and Neurosurg. 10 (3), 47–51. doi:10.19080/oajnn.2019.10.555787

CrossRef Full Text | Google Scholar

Rowson, S., Duma, S. M., Beckwith, J. G., Chu, J. J., Greenwald, R. M., Crisco, J. J., et al. (2012). Rotational head kinematics in football impacts: an injury risk function for concussion. Ann. Biomed. Eng. 40 (1), 1–13. doi:10.1007/s10439-011-0392-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmidt, J. D., Guskiewicz, K. M., Blackburn, J. T., Mihalik, J. P., Siegmund, G. P., and Marshall, S. W. (2014). The influence of cervical muscle characteristics on head impact biomechanics in football. Am. J. sports Med. 42 (9), 2056–2066. doi:10.1177/0363546514536685

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmitt, K.-U., Niederer, P. F., Cronin, D. S., Morrison, B., Muser, M. H., and Walz, F. (2019). (Book) trauma biomechanics: an introduction to injury biomechanics. Springer.

Google Scholar

Seth, A., Hicks, J. L., Uchida, T. K., Habib, A., Dembia, C. L., Dunne, J. J., et al. (2018). OpenSim: simulating musculoskeletal dynamics and neuromuscular control to study human and animal movement. PLoS Comput. Biol. 14 (7), e1006223. doi:10.1371/journal.pcbi.1006223

PubMed Abstract | CrossRef Full Text | Google Scholar

Silvestros, P., Quarrington, R. D., Preatoni, E., Gill, H. S., Jones, C. F., and Cazzola, D. (2024). An extended neck position is likely to produce cervical spine injuries through buckling in accidental head-first impacts during rugby tackling. Ann. Biomed. Eng. 52, 3125–3139. doi:10.1007/s10439-024-03576-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Singh, A., Harmukh, A., and Ganpule, S. (2022). Investigation of role of falx and tentorium on brain simulant strain under impact loading. J. Biomechanics 144, 111347. doi:10.1016/j.jbiomech.2022.111347

PubMed Abstract | CrossRef Full Text | Google Scholar

Takhounts, E. G., Hasija, V., Ridella, S. A., Rowson, S., and Duma, S. M. (2011). “Kinematic rotational brain injury criterion (BRIC),” in Proceedings of the 22nd enhanced safety of vehicles conference. Paper.

Google Scholar

Tang, J., Zhou, Q., Shen, W., Chen, W., and Tan, P. (2023). Can we reposition finite element human body model like dummies? Front. Bioeng. Biotechnol. 11, 1176818. doi:10.3389/fbioe.2023.1176818

PubMed Abstract | CrossRef Full Text | Google Scholar

Uchida, T. K., and Delp, S. L. (2021). “Biomechanics of movement: the science of sports,” in Robotics, and rehabilitation (MIT Press).

Google Scholar

Vasavada, A. N., Li, S., and Delp, S. L. (1998). Influence of muscle morphometry and moment arms on the moment-generating capacity of human neck muscles. Spine 23 (4), 412–422. doi:10.1097/00007632-199802150-00002

PubMed Abstract | CrossRef Full Text | Google Scholar

Viano, D. C. (2005). “Head impact biomechanics in sport: implications to head and neck injury tolerances,” in IUTAM symposium on impact biomechanics: from fundamental insights to applications (Springer).

Google Scholar

Viano, D. C., and Pellman, E. J. (2005). Concussion in professional football: biomechanics of the striking player—part 8. Neurosurgery 56 (2), 266–280. doi:10.1227/01.neu.0000150035.54230.3c

PubMed Abstract | CrossRef Full Text | Google Scholar

Viano, D. C., Casson, I. R., Pellman, E. J., Bir, C. A., Zhang, L., Sherman, D. C., et al. (2005). Concussion in professional football: comparison with boxing head impacts—part 10. Neurosurgery 57 (6), 1154–1172. doi:10.1227/01.neu.0000187541.87937.d9

PubMed Abstract | CrossRef Full Text | Google Scholar

Viano, D. C., Casson, I. R., and Pellman, E. J. (2007). Concussion in professional football: biomechanics of the struck player—part 14. Neurosurgery 61 (2), 313–328. doi:10.1227/01.neu.0000279969.02685.d0

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Jiang, H., Teo, E. C., and Gu, Y. (2023). Finite element analysis of head–neck kinematics in rear-end impact conditions with headrest. Bioengineering 10 (9), 1059. doi:10.3390/bioengineering10091059

PubMed Abstract | CrossRef Full Text | Google Scholar

Whyte, T., Melnyk, A. D., Van Toen, C., Yamamoto, S., Street, J., Oxland, T. R., et al. (2020). A neck compression injury criterion incorporating lateral eccentricity. Sci. Rep. 10 (1), 7114. doi:10.1038/s41598-020-63974-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Wright, R. M. (2012). A computational model for traumatic brain injury based on an axonal injury criterion. Baltimore, Maryland. The Johns Hopkins University.

Google Scholar

Wright, R. M., Post, A., Hoshizaki, B., and Ramesh, K. T. (2013). A multiscale computational approach to estimating axonal damage under inertial loading of the head. J. neurotrauma 30 (2), 102–118. doi:10.1089/neu.2012.2418

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, T., Alshareef, A., Giudice, J. S., and Panzer, M. B. (2019). Explicit modeling of white matter axonal fiber tracts in a finite element brain model. Ann. Biomed. Eng. 47, 1908–1922. doi:10.1007/s10439-019-02239-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, X., Halldin, P., and Ghajari, M. (2022). Oblique impact responses of hybrid III and a new headform with more biofidelic coefficient of friction and moments of inertia. Front. Bioeng. Biotechnol. 10, 860435. doi:10.3389/fbioe.2022.860435

PubMed Abstract | CrossRef Full Text | Google Scholar

Zaman, R., Rifat, M. N. I., Maliha, F., Hossain, M. N. B., Akhtaruzzaman, R., and Adnan, A. (2024). Multiscale structure of brain and challenges in traumatic brain injury risk prediction. Multiscale Sci. Eng. 6, 124–146. doi:10.1007/s42493-024-00117-7

CrossRef Full Text | Google Scholar

Zhan, X., Li, Y., Liu, Y., Domel, A. G., Alizadeh, H. V., Raymond, S. J., et al. (2021). The relationship between brain injury criteria and brain strain across different types of head impacts can be different. J. R. Soc. Interface 18 (179), 20210260. doi:10.1098/rsif.2021.0260

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J.-G., Wang, F., Zhou, R., and Xue, Q. (2011). A three-dimensional finite element model of the cervical spine: an investigation of whiplash injury. Med. and Biol. Eng. and Comput. 49, 193–201. doi:10.1007/s11517-010-0708-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Zuo, C., He, K., Shao, J., and Sui, Y. (2024). Self model for embodied intelligence: modeling full-body human musculoskeletal system and locomotion control with hierarchical low-dimensional representation. IEEE: IEEE International Conference on Robotics and Automation ICRA.

CrossRef Full Text | Google Scholar

Keywords: traumatic brain injury, neck injury, injury criteria, head impact, injury biomechanics, musculoskeletal modeling, neck strength, impact biomechanics

Citation: Zaman R, Hossain MI, Raiyan AZ, Chowdhury S, Jackson A, Koster AT and Adnan A (2025) A biomechanical study of neck strength and impact dynamics on head and neck injury parameters. Front. Bioeng. Biotechnol. 13:1597267. doi: 10.3389/fbioe.2025.1597267

Received: 20 March 2025; Accepted: 14 August 2025;
Published: 03 September 2025.

Edited by:

Daniele Cafolla, Swansea University, United Kingdom

Reviewed by:

Di Xu, Fujian Provincial Hospital, China
Yan Bin, Xinjiang Medical University, China

Copyright © 2025 Zaman, Hossain, Raiyan, Chowdhury, Jackson, Koster and Adnan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Ashfaq Adnan, YWFkbmFuQHV0YS5lZHU=, YXNoZmFxYWRuYW5AZ21haWwuY29t

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.