Biomechanical contributions of upper cervical ligamentous structures in Type II odontoid fractures

Fractures of the odontoid present frequently in spinal trauma, and Type II odontoid fractures, occurring at the junction of the odontoid process and C2 vertebrae, represent the bulk of all traumatic odontoid fractures. It is currently unclear what soft-tissue stabilizers contribute to upper cervical motion in the setting of a Type II odontoid fracture, and evaluation of how concomitant injury contributes to cervical stability may inform surgical decision-making as well as allow for the creation of future, accurate, biomechanical models of the upper cervical spine. The objective of the current study was to determine the contribution of soft-tissue stabilizers in the upper cervical spine following a Type II odontoid fracture. Eight cadaveric C0-C2 specimens were evaluated using a robotic testing system with motion tracking. The unilateral facet capsule (UFC) and anterior longitudinal ligament (ALL) were serially resected to determine their biomechanical role following odontoid fracture. Range of motion (ROM) and moment at the end of intact specimen replay were the primary outcomes. We determined that fracture of the odontoid significantly increases motion and decreases resistance to intact motion for flexion–extension (FE), axial rotation (AR), and lateral bending (LB). Injury to the UFC increased AR by 3.2° and FE by 3.2°. ALL resection did not significantly increase ROM or decrease end-point moment. The UFC was determined to contribute to 19% of intact flexion resistance and 24% of intact AR resistance. Overall, we determined that Type II fracture of the odontoid is a significant biomechanical destabilizer and that concurrent injury to the UFC further increases upper cervical ROM and decreases resistance to motion in a cadaveric model of traumatic Type II odontoid fractures.     

Sensitivity of lumbar spine response to follower load and flexion moment: finite element study

The follower load (FL) combined with moments is commonly used to approximate flexed/extended posture of the lumbar spine in absence of muscles in biomechanical studies. There is a lack of consensus as to what magnitudes simulate better the physiological conditions. Considering the in-vivo measured values of the intradiscal pressure (IDP), intervertebral rotations (IVRs) and the disc loads, sensitivity of these spinal responses to different FL and flexion moment magnitudes was investigated using a 3D nonlinear finite element (FE) model of ligamentous lumbosacral spine. Optimal magnitudes of FL and moment that minimize deviation of the model predictions from in-vivo data were determined. Results revealed that the spinal parameters i.e. the IVRs, disc moment, and the increase in disc force and moment from neutral to flexed posture were more sensitive to moment magnitude than FL magnitude in case of flexion. The disc force and IDP were more sensitive to the FL magnitude than moment magnitude. The optimal ranges of FL and flexion moment magnitudes were 900–1100 N and 9.9–11.2 Nm, respectively. The FL magnitude had reverse effect on the IDP and disc force. Thus, magnitude for FL or flexion that minimizes the deviation of all the spinal parameters together from the in-vivo data can vary. To obtain reasonable compromise between the IDP and disc force, our findings recommend that FL of low magnitude must be combined with flexion moment of high intensity and vice versa.     

Statistical factorial analysis approach for parameter calibration on material nonlinearity of intervertebral disc finite element model

In the biomechanics field, material parameters calibration is significant for finite element (FE) model to ensure a legit estimation of biomechanical response. Determining an appropriate combination of calibration factors is challenging as each constitutive component responds differently. This study proposes a statistical factorial analysis approach using L₁₆(4⁵) orthogonal array to evaluate material nonlinearity and applicable calibration factor of the intervertebral disc FE model in pure moment. The calibrated model exhibits improved agreement to the experimental findings for all directions. Appropriate combination of calibration parameter reduces the estimation gap to the experimental findings, ensuring agreeable biomechanical responses.     

Estimation of in vivo inter-vertebral loading during motion using fluoroscopic and magnetic resonance image informed finite element models

Finite element (FE) models driven by medical image data can be used to estimate subject-specific spinal biomechanics. This study aimed to combine magnetic resonance (MR) imaging and quantitative fluoroscopy (QF) in subject-specific FE models of upright standing, flexion and extension. Supine MR images of the lumbar spine were acquired from healthy participants using a 0.5 T MR scanner. Nine 3D quasi-static linear FE models of L3 to L5 were created with an elastic nucleus and orthotropic annulus. QF data was acquired from the same participants who performed trunk flexion to 60° and trunk extension to 20°. The displacements and rotations of the vertebrae were calculated and applied to the FE model. Stresses were averaged across the nucleus region and transformed to the disc co-ordinate system (S1 = mediolateral, S2 = anteroposterior, S3 = axial). In upright standing S3 was predicted to be −0.7 ± 0.6 MPa (L3L4) and −0.6 ± 0.5 MPa (L4L5). S3 increased to −2.0 ± 1.3 MPa (L3L4) and −1.2 ± 0.6 MPa (L4L5) in full flexion and to −1.1 ± 0.8 MPa (L3L4) and −0.7 ± 0.5 MPa (L4L5) in full extension. S1 and S2 followed similar patterns; shear was small apart from S23. Disc stresses correlated to disc orientation and wedging. The results demonstrate that MR and QF data can be combined in a participant-specific FE model to investigate spinal biomechanics in vivo and that predicted stresses are within ranges reported in the literature.     

Effect of lumbar fasciae on the stability of the lower lumbar spine

The biomechanical effect of tensioning the lumbar fasciae (LF) on the stability of the spine during sagittal plane motion was analysed using a validated finite element model of the normal lumbosacral spine (L4-S1). To apply the tension in the LF along the direction of the fibres, a local coordinate was allocated using dummy rigid beam elements that originated from the spinous process. Up to 10 Nm of flexion and 7.5 Nm of extension moment was applied with and without 20 N of lateral tension in the LF. A follower load of 400 N was additionally applied along the curvature of the spine. To identify how the magnitude of LF tension related to the stability of the spine, the tensioning on the fasciae was increased up to 40 N with an interval of 10 N under 7.5 Nm of flexion/extension moment. A fascial tension of 20 N produced a 59% decrease in angular motion at 2.5 Nm of flexion moment while there was a 12.3% decrease at 10 Nm in the L5-S1 segment. Its decrement was 53 and 9.6% at 2.5 Nm and 10 Nm, respectively, in the L4-L5 segment. Anterior translation was reduced by 12.1 and 39.0% at the L4-L5 and L5-S1 segments under 10 Nm of flexion moment, respectively. The flexion stiffness shows an almost linear increment with the increase in fascial tension. The results of this study showed that the effect of the LF on the stability of the spine is significant.     

Resultant lower extremity joint moments in below-knee amputees during running stance

Resultant flexion/extension lower extremity joint moments of four below-knee amputees running between 2.5 and 5.7 m s-1 were computed during stance on their intact and prosthetic limbs. All subjects wore patellar tendon-bearing prostheses with either a SACH or Greissinger foot component. During stance on the prosthesis, the resultant hip extensor moment on the amputated side was greater in magnitude and duration than its counterpart on the intact limb during its corresponding stance period. Since the artificial foot was planted on the ground, such a moment may help control knee flexion and promote knee extension of the residual limb. For the three subjects whose knees continued to flex at the beginning of stance, there was a dominant extensor moment about the knee joint during stance on the prosthesis. By contrast, for the fourth subject whose knee remained straight or hyperextended throughout stance on the prosthesis, a flexor moment was dominant.     

An accurate finite element model of the cervical spine under quasi-static loading

Cervical disc injury due to impact has been observed in clinical and biomechanical investigations; however, there is a lack of data that helps to elucidate the mechanisms of disc injury during these collisions. Therefore, it is necessary to understand the behavior of the cervical spine under different types of loading situations. A three dimensional finite element (FE) model for the multi-level cervical spine segment (C0-C7) was developed using computed tomography (CT) data and applied to study the internal stresses and strains of the intervertebral discs under quasi-static loading conditions. The intervertebral discs were treated as nonlinear, anisotropic and incompressible subjected to large deformations. The model accuracy was validated by comparing it with previously published experimental and numerical results for different movements. It was shown that the use of a fiber reinforced model to describe the behavior of the annulus of the discs would predict higher maximum shear strains than an isotropic one, being therefore important the use of complex constitutive models in order to be able to detect the appearance of injured zones, since those strains and stresses are supposed to be related with damage to soft tissues. Several movements were analyzed: flexion, extension and axial rotation, obtaining that the maximum shear stresses in the disc were higher for a flexo-extension movement.     

Numerical analysis of the influence of nucleus pulposus removal on the biomechanical behavior of a lumbar motion segment

Nucleus replacement was deemed to have therapeutic potential for patients with intervertebral disc herniation. However, whether a patient would benefit from nucleus replacement is technically unclear. This study aimed to investigate the influence of nucleus pulposus (NP) removal on the biomechanical behavior of a lumbar motion segment and to further explore a computational method of biomechanical characteristics of NP removal, which can evaluate the mechanical stability of pulposus replacement. We, respectively, reconstructed three types of models for a mildly herniated disc and three types of models for a severely herniated disc based on a L4–L5 segment finite element model with computed tomography image data from a healthy adult. First, the NP was removed from the herniated disc models, and the biomechanical behavior of NP removal was simulated. Second, the NP cavities were filled with an experimental material (Poisson's ratio = 0.3; elastic modulus = 3 MPa), and the biomechanical behavior of pulposus replacement was simulated. The simulations were carried out under the five loadings of axial compression, flexion, lateral bending, extension, and axial rotation. The changes of the four biomechanical characteristics, i.e. the rotation degree, the maximum stress in the annulus fibrosus (AF), joint facet contact forces, and the maximum disc deformation, were computed for all models. Experimental results showed that the rotation range, the maximum AF stress, and joint facet contact forces increased, and the maximum disc deformation decreased after NP removal, while they changed in the opposite way after the nucleus cavities were filled with the experimental material.     

Finite element analysis of moment-rotation relationships for human cervical spine

A comprehensive, geometrically accurate, nonlinear C0-C7 FE model of head and cervical spine based on the actual geometry of a human cadaver specimen was developed. The motions of each cervical vertebral level under pure moment loading of 1.0 Nm applied incrementally on the skull to simulate the movements of the head and cervical spine under flexion, tension, axial rotation and lateral bending with the inferior surface of the C7 vertebral body fully constrained were analysed. The predicted range of motion (ROM) for each motion segment were computed and compared with published experimental data. The model predicted the nonlinear moment-rotation relationship of human cervical spine. Under the same loading magnitude, the model predicted the largest rotation in extension, followed by flexion and axial rotation, and least ROM in lateral bending. The upper cervical spines are more flexible than the lower cervical levels. The motions of the two uppermost motion segments account for half (or even higher) of the whole cervical spine motion under rotational loadings. The differences in the ROMs among the lower cervical spines (C3-C7) were relatively small. The FE predicted segmental motions effectively reflect the behavior of human cervical spine and were in agreement with the experimental data. The C0-C7 FE model offers potentials for biomedical and injury studies.     

Contact pressure in the facet joint during sagittal bending of the cadaveric cervical spine

The facet joint contributes to the normal biomechanical function of the spine by transmitting loads and limiting motions via articular contact. However, little is known about the contact pressure response for this joint. Such information can provide a quantitative measure of the facet joint's local environment. The objective of this study was to measure facet pressure during physiologic bending in the cervical spine, using a joint capsule-sparing technique. Flexion and extension bending moments were applied to six human cadaveric cervical spines. Global motions (C2-T1) were defined using infra-red cameras to track markers on each vertebra. Contact pressure in the C5-C6 facet was also measured using a tip-mounted pressure transducer inserted into the joint space through a hole in the postero-inferior region of the C5 lateral mass. Facet contact pressure increased by 67.6 ± 26.9 kPa under a 2.4 Nm extension moment and decreased by 10.3 ± 9.7 kPa under a 2.7 Nm flexion moment. The mean rotation of the overall cervical specimen motion segments was 9.6 ± 0.8° and was 1.6 ± 0.7° for the C5-C6 joint, respectively, for extension. The change in pressure during extension was linearly related to both the change in moment (51.4 ± 42.6 kPa/Nm) and the change in C5-C6 angle (18.0 ± 108.9 kPa/deg). Contact pressure in the inferior region of the cervical facet joint increases during extension as the articular surfaces come in contact, and decreases in flexion as the joint opens, similar to reports in the lumbar spine despite the difference in facet orientation in those spinal regions. Joint contact pressure is linearly related to both sagittal moment and spinal rotation. Cartilage degeneration and the presence of meniscoids may account for the variation in the pressure profiles measured during physiologic sagittal bending. This study shows that cervical facet contact pressure can be directly measured with minimal disruption to the joint and is the first to provide local pressure values for the cervical joint in a cadaveric model.     

Effect of alignment changes on socket reaction moments during gait in transfemoral and knee-disarticulation prostheses: Case series

The alignment of a lower-limb prosthesis is critical to the successful prosthetic fitting and utilization by the wearer. Loads generated by the socket applied to the residual limb while walking are thought to be different in transfemoral and knee-disarticulation prostheses. The aim of this case series was to compare the socket reaction moments between transfemoral and knee-disarticulation prostheses and to investigate the effect of alignment changes on them. Two amputees, one with a transfemoral prosthesis and another with a knee-disarticulation prosthesis, participated in this study. A Smart Pyramid™ was used to measure socket reaction moments while walking under 9 selected alignment conditions; including nominally aligned, angle malalignments of 6° (flexion, extension, abduction and adduction) and translation malalignments of 15 mm (anterior, posterior, medial and lateral) of the socket relative to the foot. This study found that the pattern of the socket reaction moments was similar between transfemoral and knee-disarticulation prostheses. An extension moment in the sagittal plane and a varus moment in the coronal plane were dominant during stance under the nominally aligned condition. This study also demonstrated that alignment changes might have consistent effects on the socket reaction moments in transfemoral and knee-disarticulation prostheses. Extension and posterior translation of the socket resulted in increases in an extension moment, while abduction and lateral translation of the socket resulted in increases in a varus moment. The socket reaction moments may potentially serve as useful biomechanical parameters to evaluate alignment in transfemoral and knee-disarticulation prostheses.     

Influence of cervical disc degeneration after posterior surgical techniques in combined flexion-extension--a nonlinear analytical study

Laminectomy and facetectomy are surgical techniques used for decompression of the cervical spinal stenosis. Recent in vitro and finite element studies have shown significant cervical spinal instability after performing these surgical techniques. However, the influence of degenerated cervical disk on the biomechanical responses of the cervical spine after these surgical techniques remains unknown. Therefore, a three-dimensional nonlinear finite element model of the human cervical spine (C2-C7) was created. Two types of disk degeneration grades were simulated. For each grade of disk degeneration, the intact as well as the two surgically altered models simulating C5 laminectomy with or without C5-C6 total facetectomies were exercised under flexion and extension. Intersegmental rotational motions, internal disk annulus, cancellous and cortical bone stresses were obtained and compared to the normal intact model. Results showed that the cervical rotational motion decreases with progressive disk degeneration. Decreases in the rotational motion due to disk degeneration were accompanied by higher cancellous and cortical bone stress. The surgically altered model showed significant increases in the rotational motions after laminectomies and facetectomies when compared to the intact model. However, the percentage increases in the rotational motions after various surgical techniques were reduced with progressive disk degeneration.     

Relationships between individual isometric muscle forces, EMG activity and joint torque in monkeys

The aim of the present work was to determine the EMG activity and the moment of force developed by the main elbow flexor muscles, and to establish on this basis the degree of their participation in isometric contractions performed at various positions of the elbow. This was achieved by recording the following biomechanical parameters: EMG and tensile stress (or force) from biceps brachii (BB) and brachioradialis (BR); EMG from brachialis; external resultant force (FE). There was: a linear or quadratic relationship between the integrated EMG from each muscle and FE; a linear relationship between the force produced by BB or BR and FE. The slope of these relationships depended on the elbow angle, except for that between BB force and FE. It is proposed that iEMG changes compensate for those of the force lever arm. It has been calculated that the contribution of BR to external torque decreased from the extension to flexion while that of BB increased from 70 degrees to 90 degrees and then decreased. How far these data can be extrapolated to man is a matter of discussion based on iEMG and anthropometrical data.     

Effect of annular incision type on the change in biomechanical properties in a herniated lumbar intervertebral disc

The technique used to incise the disc during discectomy may play a role in the subsequent healing and change in biomechanical stiffness of the disc. Several techniques of lumbar disc annulotomy have been described in clinical reports. The purpose of this paper was to study the influence of annulotomy technique on motion segment stiffness using a finite element model. Four incision methods (square, circular, cross, and slit) were compared. The analyses showed that each of the annular incisions produced increase in motions under axial moment loadings with circular incision producing the largest change in the corresponding rotational motion. Under shear loading mode, cross and slit-type annular incisions produced slightly larger changes in the principal motions of the disc than square and circular incisions. All other incision types considered in the current study produced negligibly small increase in motion under rest of the loading conditions. In addition to annulotomy, when nucleotomy was also included in the analyses, once again cross and slit incisions produced larger change in motion under shear loading mode as compared to the other two incision types. A comparison between the four types of annular incisions showed that cross incision produced an increase in motion larger than those produced by the other three incisions under flexion/extension and lateral moment loading and both shear force loadings. Circular incision produced the largest increase in motion under axial moment load in comparison to those produced by square, cross, and slit incisions. Sagittal plane symmetry was influenced by the incision injury to the motion segment leading to coupled motions as well as increased facet loads. From the study it can be concluded that the increase inflexibility of the disc due to annulotomy depends on the type of annulotomy and the annulotomy also produce asymmetrical deformations leading to increased facet loading.     

Biomechanical effect after Coflex and Coflex rivet implantation for segmental instability at surgical and adjacent segments: a finite element analysis

The Coflex device may provide stability to the surgical segment in extension but does not restore stability in other motion. Recently, a modified version called the Coflex rivet has been developed. The effects of Coflex and Coflex rivet implantation on the adjacent segments are still not clear; therefore, the purpose of this study was to investigate the biomechanical differences between Coflex and Coflex rivet implantation by using finite element analyses. The results show that the Coflex implantation can provide stability in extension, lateral bending, and axial rotation at the surgical segment, and it had no influence at adjacent segments except for extension. The Coflex rivet implantation can provide stability in all motions and reduce disc annulus stress at the surgical segment. Therefore, the higher range of motion and stress induced by the Coflex rivet at both adjacent discs may result in adjacent segment degeneration in flexion and extension.     

Biomechanical effect after Coflex and Coflex rivet implantation for segmental instability at surgical and adjacent segments: a finite element analysis

The Coflex device may provide stability to the surgical segment in extension but does not restore stability in other motion. Recently, a modified version called the Coflex rivet has been developed. The effects of Coflex and Coflex rivet implantation on the adjacent segments are still not clear; therefore, the purpose of this study was to investigate the biomechanical differences between Coflex and Coflex rivet implantation by using finite element analyses. The results show that the Coflex implantation can provide stability in extension, lateral bending, and axial rotation at the surgical segment, and it had no influence at adjacent segments except for extension. The Coflex rivet implantation can provide stability in all motions and reduce disc annulus stress at the surgical segment. Therefore, the higher range of motion and stress induced by the Coflex rivet at both adjacent discs may result in adjacent segment degeneration in flexion and extension.     

Development of a model based method for investigating facet articulation

Reported investigations of facet articulation in the human spine have often been conducted through the insertion of pressure sensitive film into the joint space, which requires incision of the facet capsule and may alter the characteristics of interaction between the facet surfaces. Load transmission through the facet has also been measured using strain gauges bonded to the articular processes. While this method allows for preservation of the facet capsule, it requires extensive instrumentation of the spine, as well as strain-gauge calibration, and is highly sensitive to placement and location of the strain gauges. The inherently invasive nature of these techniques makes it difficult to translate them into medical practice. A method has been developed to investigate facet articulation through the application of test kinematics to a specimen-specific rigid-body model of each vertebra within a lumbar spine segment. Rigid-body models of each vertebral body were developed from CT scans of each specimen. The distances between nearest-neighboring points on each facet surface were calculated for specific time frames of each specimen's flexion/extension test. A metric describing the proportion of each facet surface within a distance (2 mm) from the neighboring surface, the contact area ratio (CAR), was calculated at each of these time frames. A statistically significant difference (p<0.037) was found in the CAR between the time frames corresponding to full flexion and full extension in every level of the lumbar spine (L1-L5) using the data obtained from the seven specimens evaluated in this study. The finding that the contact area of the facet is greater in extension than flexion corresponds to other findings in the literature, as well as the generally accepted role of the facets in extension. Thus, a biomechanical method with a sufficiently sensitive metric is presented as a means to evaluate differences in facet articulation between intact and treated or between healthy and pathologic spines.     

The effect of degenerative morphological changes of the intervertebral disc on the lumbar spine biomechanics: a poroelastic finite element investigation

Intervertebral disc degeneration involves changes in the spinal anatomical structures. The mechanical relevance of the following changes was investigated: disc height, endplate sclerosis, disc water content, permeability and depressurisation. A poroelastic nonlinear finite element model of the L4–L5 human spine segments was employed. Loads represented a daily cycle (500 N compression combined with flexion–extension motion for 16 h followed by 200 N compression for 8 h). In non-degenerative conditions, the model predicted a diurnal axial displacement of 1.32 mm and a peak intradiscal pressure of 0.47 MPa. Axial displacement, facet force and range of motion in flexion–extension are decreased by decreasing disc height. By decreasing the initial water content, axial displacement, facet force and fluid loss were all reduced. Endplate sclerosis did not have a significant influence on the calculated results. Depressurisation determined an increase of the disc effective stress, possibly inducing failure. Degenerative instability was not calculated in any simulations.     

A multibody modelling approach to determine load sharing between passive elements of the lumbar spine

The human spinal segment is an inherently complex structure, a combination of flexible and semi-rigid articulating elements stabilised by seven principal ligaments. An understanding of how mechanical loading is shared among these passive elements of the segment is required to estimate tissue failure stresses. A 3D rigid body model of the complete lumbar spine has been developed to facilitate the prediction of load sharing across the passive elements. In contrast to previous multibody models, this model includes a non-linear, six degrees of freedom intervertebral disc, facet bony articulations and all spinal ligaments. Predictions of segmental kinematics and facet joint forces, in response to pure moment loading (flexion–extension), were compared to published in vitro data. On inclusion of detailed representation of the disc and facets, the multibody model fully captures the non-linear flexibility response of the spinal segment, i.e. coupled motions and a mobile instantaneous centre of rotation. Predicted facet joint forces corresponded well with reported values. For the loading case considered, the model predicted that the ligaments are the main stabilising elements within the physiological motion range; however, the disc resists a greater proportion of the applied load as the spine is fully flexed. In extension, the facets and capsular ligaments provide the principal resistance. Overall patterns of load distribution to the spinal ligaments are in agreement with previous predictions; however, the current model highlights the important role of the intraspinous ligament in flexion and the potentially high risk of failure. Several important refinements to the multibody modelling of the passive elements of the spine have been described, and such an enhanced passive model can be easily integrated into a full musculoskeletal model for the prediction of spinal loading for a variety of daily activities.     

Realistic loading conditions for upper body bending

Different modes of load applications are used to simulate flexion and extension of the upper body. It is not clear which loading modes deliver realistic results and allow the comparison of different studies.In a numerical study, a validated finite element model of the lumbar spine, ranging from the vertebra L1 to the disc L5–S1 was employed. Each of six different loading modes was studied for simulating flexion and extension, including pure moments, an eccentric axial force, using a wedged fixture, and applying upper body weight plus follower load plus muscle forces. Intersegmental rotations, intradiscal pressures and facet joint contact forces were calculated. Where possible, results were compared to data measured in vivo.The results of the loading modes studied show a large variance for some values. Outcome measures such as flexion angle and intradiscal pressure differed at a segment by up to 44% and 88%, respectively, related to their maximum values. Intradiscal pressure is mainly determined by the magnitude of the applied compressive force. For flexion maximum contact forces between 0 and 69 N are predicted in each facet joint for different loading modes. For both flexion and extension, applying upper body weight plus follower load plus muscle forces as well as a follower load together with a bending moment delivers results which agreed well with in vivo data from the literature.Choosing an adequate loading mode is important in spine biomechanics when realistic results are required for intersegmental rotations, intradiscal pressure and facet joint contact forces. Only then will results of different studies be comparable.