The effects of different articulate curvature of artificial disc on loading distribution

Purpose: Deeper insights into the mechanical behavior of lumbar disc prostheses are required. Prior studies on the biomechanical performance of artificial discs were mostly performed with finite element analyses, but this has never been analyzed with altering articulate curvature. This study aimed to ascertain the influence of the geometry of a ball-and-socket disc prosthesis for the lumbar spine. Materials and Methods: Three-dimensional finite element model of human L4-L5 was reconstructed. Convex, concave, and elliptic artificial disc models were also established with Computer-Aided-Design software. Simulations included: (1) three articulate types of polyethylene (PE) insert were implanted inferiorly and (2) concave and convex PE inserts were implanted on the superior or inferior sides in flexion/extension, lateral bending, and axial rotation in the lumbar spine. Shear stresses and von Mises stresses on PE insert were assessed for their loading distributions. Results: High shear stresses of all articulate types occurred in flexion, and convex PE insert performed the maximum stress of 23.81 MPa. Under all conditions, stresses on concave PE inserts were distributed more evenly and lower than those on the convex type. Elliptic geometry enabled confining the rotation of the motion unit. The shear force on the convex PE insert on the inferior side could induce transverse crack because the shear stress exceeded yielding shear stress. Conclusions: The concave PE insert on the inferior side not only decreased loading concentration but had relatively low stress. Such a design may be applicable for artificial discs.     

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.     

Measurement of geometric deformation of lumbar intervertebral discs under in-vivo weightbearing condition

Quantitative data of spinal intervertebral disc deformation is instrumental for investigation of spinal disc pathology. In this study, we employed a combined dual fluoroscopic imaging system and the MR imaging technique to determine the lumbar disc deformation in living human subjects. Discs at L2-3, L3-4 and L4-5 levels were investigated in 8 normal subjects. The geometric deformation of the discs under full body weight loading condition (upright standing) was determined using the supine, non-weightbearing condition as a reference. The average maximum tensile deformation was ?21% in compression and 24% in tension, and maximum shear deformation on the disc surface reached 26%. The data indicated that different portions of the disc are under different tensile and shear deformation. Further, discs of L2-3, L3-4 and L4-5 have different deformation behavior under the physiological weightbearing condition. In general, the higher level discs have higher deformation values. The technique used in this study can be used to investigate the deformation behaviors of diseased discs as well as the efficacy of different surgical modalities at restoring normal disc deformation patterns.     

Load-displacement properties of lower cervical spine motion segments

The load-displacement behavior of 35 fresh adult cervical spine motion segments was measured in compression, shear, flexion, extension, lateral bending and axial torsion tests. Motion segments were tested both intact and with posterior elements removed. Applied forces ranged to 73.6 N in compression and to 39 N in shear, while applied moments ranged to 2.16 Nm. For each mode of loading, principal and coupled motions were measured and stiffnesses were calculated. The effect of disc degeneration on motion segment stiffnesses and the moments required for motion segment failure were also measured. In compression, the stiffnesses of the cervical motion segments were similar to those of thoracic and lumbar motion segments. In other modes of loading, cervical stiffnesses were considerably smaller than thoracic or lumbar stiffnesses. Removal of the posterior elements decreased cervical motion segment stiffnesses by as much as 50%. Degenerated cervical discs were less stiff in compression and stiffer in shear than less degenerated discs, but in bending or axial torsion, no statistically significant differences were evident. Bending moments causing failure were an order of magnitude lower than those for lumbar segments.     

Biomechanical comparison between fusion of two vertebrae and implantation of an artificial intervertebral disc

Surgical treatments for lower back pain can be distributed into two main groups: fusion (arthrodesis) and disc replacement (arthroplasty). The objective of this study was to compare, under severe loading conditions, the biomechanics of the lumbar spine treated either by fusion or total disc replacement (TDR). A three-dimensional model of a two-level ligamentous lumbar segment was created and simulated through static analyses with the finite-element method (FEM) software ABAQUS. The model was validated by comparing mobility, pressure on the facets, force in the ligaments, maximum stresses, disc bulge, and endplate deflection with measured data given in the literature. The FEM analysis predicted that the mobility of the model after arthrodesis on the upper level was reduced in all rotational degrees of freedom by an average of approximately 44%, relative to healthy normal discs. Conversely, the mobility of the model after TDR on the upper level was increased in all rotational degrees of freedom by an average of approximately 52%. The level implanted with the artificial disc showed excessive ligament tensions (greater than 500 N), high facet pressures (greater than 3 MPa), and a high risk of instability. The mobility and the stresses in the level adjacent to the arthroplasty were also increased. In conclusion, the model for an implanted movable artificial disc illustrated complications common to spinal arthroplasty and showed greater risk of instability and further degeneration than predicted for the fused model. This modeling technique provides an accurate means for assessing potential biomechanical risks and can be used to improve the design of future artificial intervertebral discs.     

Influence of single-level lumbar degenerative disc disease on the behavior of the adjacent segments—A finite element model study

The current study investigated mechanical predictors for the development of adjacent disc degeneration. A 3-D finite element model of a lumbar spine was modified to simulate two grades of degeneration at the L4–L5 disc. Degeneration was modeled by changes in geometry and material properties. All models were subjected to follower preloads of 800 N and moment loads in the three principal directions of motion using a hybrid protocol. Degeneration caused changes in the loading and motion patterns of the segments above and below the degenerated disc. At the level (L3–L4) above the degenerated disc, the motion increased due to moderate degeneration by 21% under lateral bending, 26% under axial rotation and 28% under flexion/extension. At the level (L5-S1) below the degenerated disc, motion increased only during lateral bending by 20% due to moderate degeneration. Both the L3–L4 and L5-S1 segment showed a monotonic increase in both the maximum von Mises stress and shear stress in the annulus as degeneration progressed for all loading directions, expect extension at L3–L4. The most significant increase in stress was observed at the L5-S1 level during axial rotation with nearly a ten-fold increase in the maximum shear stress and 103% increase in the maximum von Mises stress. The L5-S1 segment also showed a progressive increase in facet contact force for all loading directions with degeneration. Nucleus pressure did not increase significantly for any loading direction at either the caudal or cephalic adjacent segment. Results suggest that single-level degeneration can increase the risk for injury at the adjacent levels.     

Response analysis of the lumbar spine during regular daily activities—A finite element analysis

A non-linear poroelastic finite element model of the lumbar spine was developed to investigate spinal response during daily dynamic physiological activities. Swelling was simulated by imposing a boundary pore pressure of 0.25 MPa at all external surfaces. Partial saturation of the disc was introduced to circumvent the negative pressures otherwise computed upon unloading. The loading conditions represented a pre-conditioning full day followed by another day of loading: 8 h rest under a constant compressive load of 350 N, followed by 16 h loading phase under constant or cyclic compressive load varying in between 1000 and 1600 N. In addition, the effect of one or two short resting periods in the latter loading phase was studied.The model yielded fairly good agreement with in-vivo and in-vitro measurements. Taking the partial saturation of the disc into account, no negative pore pressures were generated during unloading and recovery phase. Recovery phase was faster than the loading period with equilibrium reached in only ～3 h. With time and during the day, the axial displacement, fluid loss, axial stress and disc radial strain increased whereas the pore pressure and disc collagen fiber strains decreased. The fluid pressurization and collagen fiber stiffening were noticeable early in the morning, which gave way to greater compression stresses and radial strains in the annulus bulk as time went by. The rest periods dampened foregoing differences between the early morning and late in the afternoon periods. The forgoing diurnal variations have profound effects on lumbar spine biomechanics and risk of injury.     

High pressures and asymmetrical stresses in the scoliotic disc in the absence of muscle loading

Background

Loads acting on scoliotic spines are thought to be asymmetric and involved in progression of the scoliotic deformity; abnormal loading patterns lead to changes in bone and disc cell activity and hence to vertebral body and disc wedging. At present however there are no direct measurements of intradiscal stresses or pressures in scoliotic spines. The aim of this study was to obtain quantitative measurements of the intradiscal stress environment in scoliotic intervertebral discs and to determine if loads acting across the scoliotic spine are asymmetric. We performed in vivo measurements of stresses across the intervertebral disc in patients with scoliosis, both parallel (termed horizontal) and perpendicular (termed vertical) to the end plate, using a side mounted pressure transducer (stress profilometry)

Methods

Stress profilometry was used to measure horizontal and vertical stresses at 5 mm intervals across 25 intervertebral discs of 7 scoliotic patients during anterior reconstructive surgery. A state of hydrostatic pressure was defined by identical horizontal and vertical stresses for at least two consecutive readings. Results were compared with similar stress profiles measured during surgery across 10 discs of 4 spines with no lateral curvature and with data from the literature.

Results

Profiles across scoliotic discs were very different from those of normal, young, healthy discs of equivalent age previously presented in the literature. Hydrostatic pressure regions were only seen in 14/25 discs, extended only over a short distance. Non-scoliotic discs of equivalent age would be expected to show large centrally placed hydrostatic nuclear regions in all discs. Mean pressures were significantly greater (0.25 MPa) than those measured in other anaesthetised patients (<0.07 MPa). A stress peak was seen in the concave annulus in 13/25 discs. Stresses in the concave annulus were greater than in the convex annulus indicating asymmetric loading in these anaesthetised, recumbent patients.

Conclusion

Intradiscal pressures and stresses in scoliotic discs are abnormal, asymmetrical and high in magnitude even in the absence of significant applied muscle loading. The origin of these abnormal stresses is unclear.     

Viscoelastic properties of the central region of porcine temporomandibular joint disc in shear stress-relaxation

In this study, shear relaxation properties of the porcine temporomandibular joint (TMJ) disc are investigated. Previous studies have shown that, in fatigue failure and damage of cartilage and fibrocartilage, shear loads could be one of the biggest contributors to the failure. The aim of the present study is to develop an evaluation method to study shear properties of the disc and to do a mathematical characterization of it. For the experiments, twelve porcine discs were used. Each disc was dissected from the TMJ and, then, static strain control tests were carried out to obtain the shear relaxation modulus for the central region of the discs. From the results, it was found that the disc presents a viscoelastic behavior under shear loads. Relaxation modulus decreased with time. Shear relaxation was 10% of the instantaneous stress, which implies that the viscous properties of the disc cannot be neglected. The present results lead to a better understanding of the discs mechanical behavior under realistic TMJ working conditions.     

Direct measurement of intervertebral disc maximum shear strain in six degrees of freedom: motions that place disc tissue at risk of injury

Human intervertebral disc specimens were tested to determine the regions of largest maximum shear strain (MSS) experienced by disc tissues in each of three principal displacements and three rotations, and to identify the physiological rotations and displacements that may place the disc at greatest risk for large tissue strains and injury. Tearing of disc annulus may be initiated by large interlamellar shear strains. Nine human lumbar discs were tagged with radiographic markers on the endplates, disc periphery and with a grid of wires in the mid-transverse plane and subjected to each of the six principal displacements and rotations. Stereo-radiographs were taken in each position and digitized for reconstruction of the three-dimensional position of each marker. Maximum tissue shear strains were calculated from relative marker displacements and normalized by the input displacement or rotation. Lateral shear, compression, and lateral bending were the motions that produced the mean (95% confidence interval) largest mean MSS of 9.6 (0.7)%/mm, 9.0 (0.5)%/mm, and 5.8 (1.6)%/ degrees , respectively, and which occurred in the posterior, posterolateral and lateral peripheral regions of the disc. After taking into account the reported maximum physiological range of motion for each degree of freedom, motions producing the highest physiological MSS were lateral bending (57.8 (16.2)%) and flexion (38.3 (3.3)%), followed by lateral shear (14.4 (1.1)%) and compression (12.6 (0.7)%).     

Effects of several temporomandibular disorders on the stress distributions of temporomandibular joint: a finite element analysis

The aim of this study was to evaluate stress distributions in the temporomandibular joints (TMJs) with temporomandibular disorders (TMDs) for comparison with healthy TMJs. A model of mandible and normal TMJs was developed according to CT images. The interfaces between the discs and the articular cartilages were treated as contact elements. Nonlinear cable elements were used to simulate disc attachments. Based on this model, seven models of various TMDs were established. The maximum stresses of the discs with anterior, posterior, medial and lateral disc displacement (ADD, PDD, MDD and LDD) were 12.09, 9.33, 10.71 and 6.07 times magnitude of the identically normal disc, respectively. The maximum stresses of the posterior articular eminences in ADD, PDD, MDD, LDD, relaxation of posterior attachments and disc perforation models were 21, 59, 46, 21, 13 and 15 times greater than the normal model, respectively. TMDs could cause increased stresses in the discs and posterior articular eminences.     

C4–C5 segment finite element model development,validation, and load-sharing investigation

Detailed cervical spine models are necessary to better understand cervical spine response to loading, improve our understanding of injury mechanisms, and specifically for predicting occupant response and injury in auto crash scenarios. The focus of this study was to develop a C4–C5 finite element model with accurate representations of each tissue within the segment. This model incorporates more than double the number of elements of existing models, required for accurate prediction of response. The most advanced material data available were then incorporated using appropriate nonlinear constitutive models to provide accurate predictions of response at physiological levels of loading. This tissue-scale segment model was validated against a wide variety of experimental data including different modes of loading (axial rotation, flexion, extension, lateral bending, and translation), and different load levels. In general, the predicted response of the model was within the single standard deviation response corridors for both low and high load levels. Importantly, this model demonstrates that appropriate refinement of the finite element mesh, representation at the tissue level, and sufficiently detailed material properties and constitutive models provide excellent response predictions without calibration of the model to experimental data. Load sharing between the disc, ligaments, and facet joints was investigated for various modes of loading, and the dominant load-bearing structure was found to correlate with typical anatomical injury sites for these modes of loading. The C4–C5 model forms the basis for the development of a full cervical spine model. Future studies will focus on tissue-level injury prediction and dynamic response.     

Morphometric analysis of vertebrae and intervertebral discs as a basis of disc replacement

In this study measurements were obtained from 359 vertebrae and 215 intervertebral discs in an attempt to classify discs by their size. At the cervical and thoracic levels, this attempt was unproductive because of extensive variations. In the lumbar spine, discs were allocated to six size-matched groups and to two height-matched groups. The breadth of marginal rims were also measured with a view to provide surgeons operating on the spine with precise data on disc morphology to facilitate disc replacement.     

Investigation of whiplash injuries in the upper cervical spine using a detailed neck model

Whiplash injuries continue to have significant societal cost; however, the mechanism and location of whiplash injury is still under investigation. Recently, the upper cervical spine ligaments, particularly the alar ligament, have been identified as a potential whiplash injury location. In this study, a detailed and validated explicit finite element model of a 50th percentile male cervical spine in a seated posture was used to investigate upper cervical spine response and the potential for whiplash injury resulting from vehicle crash scenarios. This model was previously validated at the segment and whole spine levels for both kinematics and soft tissue strains in frontal and rear impact scenarios. The model predicted increasing upper cervical spine ligament strain with increasing impact severity. Considering all upper cervical spine ligaments, the distractions in the apical and alar ligaments were the largest relative to their failure strains, in agreement with the clinical findings. The model predicted the potential for injury to the apical ligament for 15.2 g frontal or 11.7 g rear impacts, and to the alar ligament for a 20.7 g frontal or 14.4 g rear impact based on the ligament distractions. Future studies should consider the effect of initial occupant position on ligament distraction.     

Inertial properties and loading rates affect buckling modes and injury mechanisms in the cervical spine

Cervical spine injuries continue to be a costly societal problem. Future advancements in injury prevention depend on improved physical and computational models which, in turn, are predicated on a better understanding of the responses of the neck during dynamic loading. Previous studies have shown that the tolerance of the neck is dependent on its initial position and its buckling behavior. This study uses a computational model to examine the mechanical factors influencing buckling behavior during impact to the neck. It was hypothesized that the inertial properties of the cervical spine influence the dynamics during compressive axial loading. The hypothesis was tested by performing parametric analyses of vertebral mass, mass moments of inertia, motion segment stiffness, and loading rate. Increases in vertebral mass resulted in increasingly complex kinematics and larger peak loads and impulses. Similar results were observed for increases in stiffness. Faster loading rates were associated with higher peak loads and higher-order buckling modes. The results demonstrate that mass has a great deal of influence on the buckling behavior of the neck, particularly with respect to the expression of higher-order modes. Injury types and mechanisms may be substantially altered by loading rate because inertial effects may influence whether the cervical spine fails in a compressive mode, or a bending mode.     

Spine biomechanics

Current trends in spine research are reviewed in order to suggest future opportunities for biomechanics. Recent studies show that psychosocial factors influence back pain behaviour but are not important causes of pain itself. Severe back pain most often arises from intervertebral discs, apophyseal joints and sacroiliac joints, and physical disruption of these structures is strongly but variably linked to pain. Typical forms of structural disruption can be reproduced by severe mechanical loading in-vitro, with genetic and age-related weakening sometimes leading to injury under moderate loading. Biomechanics can be used to quantify spinal loading and movements, to analyse load distributions and injury mechanisms, and to develop therapeutic interventions. The authors suggest that techniques for quantifying spinal loading should be capable of measurement "in the field" so that they can be used in epidemiological surveys and ergonomic interventions. Great accuracy is not required for this task, because injury risk depends on tissue weakness as much as peak loading. Biomechanical tissue testing and finite-element modelling should complement each other, with experiments establishing proof of concept, and models supplying detail and optimising designs. Suggested priority areas for future research include: understanding interactions between intervertebral discs and adjacent vertebrae; developing prosthetic and tissue-engineered discs; and quantifying spinal function during rehabilitation. "Mechanobiology" has perhaps the greatest future potential, because spinal degeneration and healing are both mediated by the activity of cells which are acutely sensitive to their local mechanical environment. Precise characterisation and manipulation of this environment will be a major challenge for spine biomechanics.     

Contact analysis of a posterior cervical spine plate using a three-dimensional canine finite element model

The development of a three-dimensional finite element model of a posteriorly plated canine cervical spine (C3-C6) including contact nonlinearities is described. The model was created from axial CT scans and the material properties were derived from the literature. The model demonstrated sufficient accuracy from the results of a mesh convergence test. Significant steps were taken toward establishing model validation by comparison of plate surface strains with a posteriorly plated canine cervical spine under three-point bending. This model was developed to better characterize the contact pressures at the various interfaces under average physiologic canine loading. The analysis showed that the screw-plate interfaces had the highest values of all the mechanical parameters evaluated.     

Search for critical loading condition of the spine–A meta analysis of a nonlinear viscoelastic finite element model

The relative vulnerability of spinal motion segments to different loading combinations remains unknown. The meta-analysis described here using the results of a validated L2–L3 nonlinear viscoelastic finite element model was designed to investigate the critical loading and its effect on the internal mechanics of the human lumbar spine. A Box-Behnken experimental design was used to design the magnitude of seven independent variables associated with loads, rotations and velocity of motion. Subsequently, an optimization method was used to find the primary and secondary variables that influence spine mechanical output related to facet forces, disc pressure, ligament forces, annulus matrix compressive/shear stresses and anulus fibers strain. The mechanical responses with respect to the two most-relevant variables were then regressed linearly using the response surface quadratic model. Axial force and sagittal rotation were identified as the most-relevant variables for mechanical responses. The procedure developed can be used to find the critical loading for finite element models with multi input variables. The derived meta-models can be used to predict the risk associated with various loading parameters and in setting safer load limits.     

Search for critical loading condition of the spine--a meta analysis of a nonlinear viscoelastic finite element model

The relative vulnerability of spinal motion segments to different loading combinations remains unknown. The meta-analysis described here using the results of a validated L2-L3 nonlinear viscoelastic finite element model was designed to investigate the critical loading and its effect on the internal mechanics of the human lumbar spine. A Box-Behnken experimental design was used to design the magnitude of seven independent variables associated with loads, rotations and velocity of motion. Subsequently, an optimization method was used to find the primary and secondary variables that influence spine mechanical output related to facet forces, disc pressure, ligament forces, annulus matrix compressive/shear stresses and anulus fibers strain. The mechanical responses with respect to the two most-relevant variables were then regressed linearly using the response surface quadratic model. Axial force and sagittal rotation were identified as the most-relevant variables for mechanical responses. The procedure developed can be used to find the critical loading for finite element models with multi input variables. The derived meta-models can be used to predict the risk associated with various loading parameters and in setting safer load limits.     

Biphasic finite element simulation of the TMJ disc from in vivo kinematic and geometric measurements

Understanding the biomechanical nature of the degeneration of the temporomandibular joint requires a coupling between experimental measurements and numerical simulation. In this study, geometry measured from MRI, and motion obtained from a specially designed optoelectronic system are fed into a three-dimensional biphasic finite element analysis to generate the spatial and temporal mechanical response of the disc. This study demonstrates how this coupling can effectively predict the biomechanical response of the temporomandibular joint disc to physiological loading. For small jaw opening movements, asymmetries in the load of the disc are found, with especially higher shear stresses in the lateral portion.