Can extra-articular strains be used to measure facet contact forces in the lumbar spine? An in-vitro biomechanical study

Empirical joint contact mechanics: A comprehensive review

Empirical joint contact mechanics measurement (EJCM; e.g. contact area or force, surface velocities) enables critical investigations of the relationship between changing joint mechanics and the impact on surface-to-surface interactions. In orthopedic biomechanics, understanding the changes to cartilage contact mechanics following joint pathology or aging is critical due to its suggested role in the increased risk of osteoarthritis (OA), which might be due to changed kinematics and kinetics that alter the contact patterns within a joint. This article reviews and discusses EJCM approaches that have been applied to articulating joints such that readers across different disciplines will be informed of the various measurement and analysis techniques used in this field. The approaches reviewed include classical measurement approaches (radiographic and sectioning, dye staining, casting, surface proximity, and pressure measurement), stereophotogrammetry/motion analysis, computed tomography (CT), magnetic resonance imaging (MRI), and high-speed videoradiography. Perspectives on approaches to advance this field of EJCM are provided, including the value of considering relative velocity in joints, tractional stress, quantification of joint contact area shape, consideration of normalization techniques, net response (superposition) of multiple input variables, and establishing linkages to regional cartilage health status. EJCM measures continue to provide insights to advance our understanding of cartilage health and degeneration and provide avenues to assess the efficacy and guide future directions of developing interventions (e.g. surgical, biological, rehabilitative) to optimize joint's health and function long term.

Comparison of the biomechanical effects of lumbar disc degeneration on normal patients and osteoporotic patients: A finite element analysis

BACKGROUND

Some older patients who suffered from both conditions (disc degeneration and osteoporosis) have higher surgical risks and longer postoperative recovery times. Understanding the relation between disc degeneration and osteoporosis is fundamental to know the mechanisms of orthopedic disorders and improve clinical treatment. However, there is a lack of finite element (FE) studies to predict the combined effects of disc degeneration and osteoporosis. So the aim of the present study is to explore the differences of biomechanical effects of lumbar disc degeneration on normal patients and osteoporotic patients.

METHODS

A normal lumbar spine finite element model (FEM) was developed based on the geometric information of a healthy male subject (age 35 years; height 178 cm; weight 65 kg). This normal lumbar spine FEM was modified to build three lumbar spine degeneration models simulating mild, moderate and severe grades of disc degeneration at the L4-L5 segment. Then the degenerative lumbar spine models for osteoporotic patients were constructed on the basis of the above-mentioned degeneration models. Firstly, the normal model (flexion: 8 Nm; extension: 6 Nm; lateral bending: 6 Nm; torsion: 4 Nm) and degenerative models (10 Nm) were calibrated under pure moment load, respectively. Secondly, under a 400 N follower load, the 7.5 Nm moments of different directions were applied on all models to simulate different motion postures. Finally, under the above loading conditions, we calculated and analyzed the range of motion (ROM), Mises stress in cortical (MSC1), Mises stress in endplate (MSE), Mises stress in cancellous (MSC2), and Mises stress in post (MSP).

RESULTS

Compared with disc degeneration patients without osteoporosis, the ROM, MSC1, and MSE of osteoporosis patients with various disc degeneration decreased in all postures, while the MSC2 and MSP increased. With increase in the degree of disc degeneration, the reduction proportions of ROM and MSE in osteoporotic patients gradually increased, while the reduction percentages in MSC1 of osteoporotic patients gradually decreased. The increase percentages of MSC2 in osteoporotic patients gradually increased. Given the progressive changes of disc degeneration, the changes in MSP in osteoporosis patients were uneven.

CONCLUSION

In summary, the effect of disc degeneration on flexibility in the two kinds of patients (osteoporosis and non-osteoporosis patients) was nearly same. By comparing the remaining biomechanical parameters (MSC1, MSE, MSC2, and MSP), we found that degenerated intervertebral discs caused changes in loading patterns of osteoporosis patients. Disc degeneration reduced the Mises stress in the cortical and endplate, which increased the Mises stress in the cancellous and post. That is to say, in order to cope with the changes in bone stresses caused by disc degeneration and osteoporosis, clinicians should be more careful in choosing the surgical option for osteoporotic patients with disc degeneration.

Yeditepe spine mesh: Finite element modeling and validation of a parametric CAD model of lumbar spine

Finite element analysis is a powerful tool that is often used to study the biomechanical response of the spine. The primary objective of this study was to illustrate the mechanical behavior of a previously proposed parametric CAD spine model in comparison with a segmented FSU model and the literature. In this study, two finite element models of the L4-L5 spinal level were developed from the same patient's CT scan data. The first was developed using well-known segmentation methods, whereas the second was developed from the new by using a novel parametric CAD model. Both models were subjected to the same loading and boundary conditions to perform flexion, extension, lateral bending and axial rotation motions. The segmented finite element model was observed to be in good agreement with the literature. The parametric finite element model results were also observed to be in good agreement with the segmented finite element model and with the literature except under extension.

Biomechanical modelling of the facet joints: a review of methods and validation processes in finite element analysis

There is an increased interest in studying the biomechanics of the facet joints. For in silico studies, it is therefore important to understand the level of reliability of models for outputs of interest related to the facet joints. In this work, a systematic review of finite element models of multi-level spinal section with facet joints output of interest was performed. The review focused on the methodology used to model the facet joints and its associated validation. From the 110 papers analysed, 18 presented some validation of the facet joints outputs. Validation was done by comparing outputs to literature data, either computational or experimental values; with the major drawback that, when comparing to computational values, the baseline data was rarely validated. Analysis of the modelling methodology showed that there seems to be a compromise made between accuracy of the geometry and nonlinearity of the cartilage behaviour in compression. Most models either used a soft contact representation of the cartilage layer at the joint or included a cartilage layer which was linear elastic. Most concerning, soft contact models usually did not contain much information on the pressure-overclosure law. This review shows that to increase the reliability of in silico model of the spine for facet joints outputs, more needs to be done regarding the description of the methods used to model the facet joints, and the validation for specific outputs of interest needs to be more thorough, with recommendation to systematically share input and output data of validation studies.

The Effects of Orientation of Lumbar Facet Joints on the Facet Joint Contact Forces: An In Vitro Biomechanical Study

STUDY DESIGN

A biomechanical human cadaveric study.

OBJECTIVE

The aim of this study was to measure L2-L3 facet joint contact forces in a flexibility test using thin film electroresistive sensors, and facet joint orientation on computed tomographic (CT) scan images, to examine the effects of orientation of lumbar facet joint on the facet joint contact forces.

SUMMARY OF BACKGROUND DATA

Biomechanically, the bilateral facet joints play a critical role in maintaining stability of the lumbar spine. The effect of orientation of lumbar facet joints on the contact forces remains unknown.

METHODS

Eight human cadaveric lumbar spine specimens (L2-L3) were tested by applying a pure moment of ±7.5 Nm in three directions of loading (flexion-extension, lateral bending, and axial rotation) with and without a follower preload of 300 N. The orientation of the lumbar facet joints at the L2-L3 was measured on axial CT scans. Bilateral facet contact forces were measured during flexibility tests using thin film electroresistive sensors (Tekscan 6900).

RESULTS

The average total peak facet loads was 66 N in axial rotation, 27 N in extension, and 20 N in lateral bending under a pure moment. Under a pure moment with a follower preload of 300 N, the average total peak facet loads was 53 N in axial rotation, 43 N in extension, and 24 N in lateral bending. The facet joint forces were correlated positively and significantly with the orientation in all directions with and without a compressive follower preload (P < 0.05). In addition, the facet joint contact forces at neutral position with a follower preload were correlated positively with the orientation (rs = 0.759, P = 0.001).

CONCLUSION

This study identified that the greater coronal orientation of lumbar facet joints is, the higher the facet joint contact forces are.

LEVEL OF EVIDENCE

3.

Impact of material and morphological parameters on the mechanical response of the lumbar spine - A finite element sensitivity study

Finite element models are frequently used to study lumbar spinal biomechanics. Deterministic models are used to reflect a certain configuration, including the means of geometrical and material properties, while probabilistic models account for the inherent variability in the population. Because model parameters are generally uncertain, their predictive power is frequently questioned. In the present study, we determined the sensitivities of spinal forces and motions to material parameters of intervertebral discs, vertebrae, and ligaments and to lumbar morphology. We performed 1200 model simulations using a generic model of the human lumbar spine loaded under pure moments. Coefficients of determination and of variation were determined for all parameter and response combinations. Material properties of the vertebrae displayed the least impact on results, whereas those of the discs and morphology impacted most. The most affected results were the axial compression forces in the vertebral body and in several ligaments during flexion and the facet-joint forces during extension. Intervertebral rotations were considerably affected only when several parameters were varied simultaneously. Results can be used to decide which model parameters require careful consideration in deterministic models and which parameters might be omitted in probabilistic studies. Findings allow quantitative estimation of a model׳s precision.

Biomechanical response of lumbar facet joints under follower preload: a finite element study

Background

Facet joints play a significant role in providing stability to the spine and they have been associated with low back pain symptoms and other spinal disorders. The influence of a follower load on biomechanics of facet joints is unknown. A comprehensive research on the biomechanical role of facets may provide insight into facet joint instability and degeneration.

Method

A nonlinear finite element (FE) model of lumbar spine (L1-S1) was developed and validated to study the biomechanical response of facets, with different values of follower preload (0 N,500 N,800 N,1200 N), under loadings in the three anatomic planes. In this model, special attention was paid to the modeling of facet joints, including cartilage layer. The asymmetry in the biomechanical response of facets was also discussed. A rate of change (ROC) and an average asymmetry factor (AAF) were introduced to explore and evaluate the preload effect on these facet contact parameters and on the asymmetry under different loading conditions.

Results

The biomechanical response of facets changed according to the loading condition. The preload amplified the facet force, contact area and contact pressure in flexion-extension; the same effect was observed on the ipsilateral facet while an opposite effect could be seen on the contralateral facet during lateral bending. For torsion loading, the preload increased contact area, decreased the mean contact pressure, but had almost no effect on facet force. However, all the effects of follower load on facet response became weaker with the increase of preload. The greatest asymmetry of facet response could be found on the ipsilateral side during lateral bending, followed by flexion, bending (contralateral side), extension and torsion. This asymmetry could be amplified by preload in the bending (ipsilateral), torsion loading group, while being reduced in the flexion group.

Conclusions

An analysis combining patterns of contact pressure distribution, facet load, contact area and contact pressure can provide more insight into the biomechanical role of facets under various moment loadings and follower loads. The effect of asymmetry on facet joint response should be fully considered in biomechanical studies of lumbar spine, especially in post structures subjected to physiological loadings.

Electronic supplementary material

The online version of this article (doi:10.1186/s12891-016-0980-4) contains supplementary material, which is available to authorized users.

A computational spinal motion segment model incorporating a matrix composition-based model of the intervertebral disc

The extracellular matrix of the intervertebral disc is subjected to changes with age and degeneration, affecting the biomechanical behaviour of the spine. In this study, a finite element model of a generic spinal motion segment that links spinal biomechanics and intervertebral disc biochemical composition was developed. The local mechanical properties of the tissue were described by the local matrix composition, i.e. fixed charge density, amount of water and collagen and their organisation. The constitutive properties of the biochemical constituents were determined by fitting numerical responses to experimental measurements derived from literature. This general multi-scale model of the disc provides the possibility to evaluate the relation between local disc biochemical composition and spinal biomechanics.

Characterization and prediction of rate-dependent flexibility in lumbar spine biomechanics at room and body temperature

BACKGROUND CONTEXT

The soft tissues of the spine exhibit sensitivity to strain-rate and temperature, yet current knowledge of spine biomechanics is derived from cadaveric testing conducted at room temperature at very slow, quasi-static rates.

PURPOSE

The primary objective of this study was to characterize the change in segmental flexibility of cadaveric lumbar spine segments with respect to multiple loading rates within the range of physiologic motion by using specimens at body or room temperature. The secondary objective was to develop a predictive model of spine flexibility across the voluntary range of loading rates.

STUDY DESIGN

This in vitro study examines rate- and temperature-dependent viscoelasticity of the human lumbar cadaveric spine.

METHODS

Repeated flexibility tests were performed on 21 lumbar function spinal units (FSUs) in flexion-extension with the use of 11 distinct voluntary loading rates at body or room temperature. Furthermore, six lumbar FSUs were loaded in axial rotation, flexion-extension, and lateral bending at both body and room temperature via a stepwise, quasi-static loading protocol. All FSUs were also loaded using a control loading test with a continuous-speed loading-rate of 1-deg/sec. The viscoelastic torque-rotation response for each spinal segment was recorded. A predictive model was developed to accurately estimate spine segment flexibility at any voluntary loading rate based on measured flexibility at a single loading rate.

RESULTS

Stepwise loading exhibited the greatest segmental range of motion (ROM) in all loading directions. As loading rate increased, segmental ROM decreased, whereas segmental stiffness and hysteresis both increased; however, the neutral zone remained constant. Continuous-speed tests showed that segmental stiffness and hysteresis are dependent variables to ROM at voluntary loading rates in flexion-extension. To predict the torque-rotation response at different loading rates, the model requires knowledge of the segmental flexibility at a single rate and specified temperature, and a scaling parameter. A Bland-Altman analysis showed high coefficients of determination for the predictive model.

CONCLUSIONS

The present work demonstrates significant changes in spine segment flexibility as a result of loading rate and testing temperature. Loading rate effects can be accounted for using the predictive model, which accurately estimated ROM, neutral zone, stiffness, and hysteresis within the range of voluntary motion.

Increase in facet joint loading after nucleotomy in the human lumbar spine

Low-back pain has been related to degenerative changes after nucleotomy. Although several etiologies for pain after nucleotomy have been proposed, there is evidence of pain arising in the facet joints in general, which may be related to changes in load transfer. This study addresses the effect of nucleotomy on facet joint loading. Nine human lumbar motion segments (age: 40-59 years) were loaded in axial compression and extension-flexion. Reaction forces were compared with soft tissue structures sequentially removed. After nucleotomy the facets supported significantly greater load, almost doubling from a median of 8.6% of the applied external force to 15.8%. Force transmission related to the capsular ligament increased significantly from an intact median of 1.2-5.1% after nucleotomy. No correlation was observed between force increase on the facets and the proportion of disc nucleus removed. Even a small quantity of nucleus removal (range: 0.7-1.7g) increased the forces transmitted over the facet joints, both with and without capsular ligaments. This suggests that the proportion of material removed might not be important clinically with regard to facet joint degeneration and pain.

Lumbar facet joint and intervertebral disc loading during simulated pelvic obliquity

BACKGROUND CONTEXT

Intervertebral disc and facet joints are the two primary load-bearing structures of the lumbar spine, and altered loading to these structures may be associated with frontal plane spinal deviations.

PURPOSE

To determine the load on the lumbar facet joint and intervertebral disc under simulated frontal plane pelvic obliquity combined loading, an in vitro biomechanical study was conducted.

STUDY DESIGN/SETTING

An in vitro biomechanical study using a repeated-measures design was used to compare L4-L5 facet joint and intervertebral disc loading across pure moment and combined loading conditions.

METHODS

Eight fresh-frozen lumbosacral specimens were tested under five loading conditions: flexion/extension, lateral bending, axial rotation using pure moment bending (±10 Nm), and two additional tests investigating frontal plane pelvic obliquity and axial rotation (sacrum tilted left 5° and at 10° followed by a ±10-Nm rotation moment). Three-dimensional kinematics, facet load, and intradiscal pressures were recorded from the L4-L5 functional spinal unit.

RESULTS

Sagittal and frontal plane loading resulted in significantly smaller facet joint forces compared with conditions implementing a rotation moment (p<.05). The facet joint had the highest peak load during the 10° combined loading condition (124.0±30.2 N) and the lowest peak load in flexion (26.8±16.1 N). Intradiscal pressure was high in lateral flexion (495.6±280.9 kPa) and flexion (429.0±212.9 kPa), whereas intradiscal pressures measured in rotation (253.2±135.0 kPa) and 5° and 10° combined loading conditions were low (255.5±132.7 and 267.1±127.1 kPa, respectively).

CONCLUSIONS

Facet loading increased during simulated pelvic obliquity in frontal and transverse planes, whereas intradiscal pressures were decreased compared with sagittal and frontal plane motions alone. Altered spinopelvic alignment may increase the loads experienced by spinal tissue, especially the facet joints.

The effect of nucleotomy on facet joint loading - a porcine in vitro study

BACKGROUND

Lumbar facet joints have been cited as a possible origin of low-back pain. A relationship between disc height decrease and facet joint degeneration has been reported. Facet joint degeneration may also be triggered by nucleotomy, performed on prolapsed discs, which might change the natural load sharing between the anterior and posterior structures of the spine. In this study load bearing of the facet joints was compared between natural and nucleotomised spinal segments.

METHODS

Nine porcine lumbar motion segments were tested quasi-statically in ±1.5° extension-flexion under 700 N constant compression loading. The kinematics of the spinal segments were recorded as a response to the applied load. These kinematics were subsequently applied to the segments with the ligaments and disc sequentially removed and the reaction forces measured. This was performed in samples with and without nucleotomy. Comparison of the reaction forces allowed a direct comparison between healthy and pathological force transmission over the facet joints. Load sharing was related to the proportion of removed nucleus.

FINDINGS

The proportion of applied compression force supported by the facets increased from a mean of 40.7% (standard deviation, SD 10.0%) to 82.0% (SD 7.2%) after nucleotomy averaged over the entire extension-flexion regime. No correlation was observed between facet loading and the proportion of the nucleus removed.

INTERPRETATION

Increased facet loading after nucleotomy might cause greater cartilage wear, which may be related to facet joint degeneration. The independence of facet loading on the proportion of nucleus removed might be due to a complete pressure loss once the annulus is incised.

Effect of centers of rotation on spinal loads and muscle forces in total disk replacement of lumbar spine

The placement of artificial disks can alter the center of rotation and kinematic pattern; therefore, forces in the spine during the motion will be affected as a result. The relationship between the location of joint center of artificial disks and forces in the spinal components is not investigated. A musculoskeletal model of the spine was developed, and three location cases of center of rotation were investigated varying 5 mm anteriorly and posteriorly from the default center. Resultant joint forces, ligament forces, facet forces, and muscle forces for each case were predicted during sagittal motion. No considerable difference was observed for joint force (maximum 14%). Anterior shift of center of rotation induced the most ligament forces (200 N) and facet forces (130 N) among the three cases. Posterior and anterior shifts of centers of rotation from the default location caused considerable changes in muscle forces, respectively: 108% and 70% of increase in multifidi muscle and 157% and 187% of increase in short segmental muscle. This study showed that the centers of rotation due to the design and the surgical placement of artificial disk can affect the kinetic results in the spine.

Facet joint contact pressure is not significantly affected by ProDisc cervical disc arthroplasty in sagittal bending: a single-level cadaveric study

BACKGROUND CONTEXT

Total disc arthroplasty is a motion-preserving spinal procedure that has been investigated for its impact on spinal motions and adjacent-level degeneration. However, the effects of disc arthroplasty on facet joint biomechanics remain undefined despite the critical role of these posterior elements on guiding and limiting spinal motion.

PURPOSE

The goal was to measure the pressure in the facet joint in cadaveric human cervical spines subjected to sagittal bending before and after implantation of the ProDisc-C (Synthes Spine Company, L.P, West Chester, PA, USA).

STUDY DESIGN

A biomechanical study was performed using cadaveric human cervical spines during sagittal bending in the intact and implanted conditions.

METHODS

Seven C2-T1 osteoligamentous cadaveric cervical spines were instrumented with a transducer to measure the C5-C6 facet pressure profiles during physiological sagittal bending, before and after implantation of a ProDisc-C at that level. Rotations of the index segment and global cervical spine were also quantified.

RESULTS

The mean C5-C6 range of motion significantly increased (p=.009) from 9.6°±5.1° in the intact condition to 16.2°±3.6° after implantation. However, despite such changes in rotation, there was no significant difference in the facet contact pressure during extension between the intact (64±30 kPa) and implanted (44±55 kPa) conditions. Similarly, there was no difference in facet pressure developed during flexion.

CONCLUSIONS

Although implantation of a ProDisc-C arthroplasty device at the C5-C6 level increases angular rotations, it does not significantly alter the local facet pressure at the index level in flexion or extension. Using a technique that preserves the capsular ligament, this study provides the first direct measurement of cervical facet pressure in a disc arthroplasty condition.

Kinematic evaluation of one- and two-level Maverick lumbar total disc replacement caudal to a long thoracolumbar spinal fusion

PURPOSE

Adjacent level degeneration that occurs above and/or below long fusion constructs is a documented clinical problem that is widely believed to be associated with the considerable change in stiffness caused by the fusion. Some researchers have suggested that early degeneration at spinal joints adjacent to a fusion could be treated by implanting total disc replacements at these levels. It is thought that further degeneration could be prevented through the disc replacement's design aims to reproduce normal disc heights, kinematics and tissue loading. For this reason, there is a clinical need to evaluate if a total disc replacement can maintain both the quantity of motion (i.e. range) and the quality of motion (i.e. center of rotation and coupling) at segments adjacent to a long spinal fusion. The purpose of this study was to experimentally evaluate range of motion (ROM-the intervertebral motion measured) and helical axis of motion (HAM) changes due to one- and two-level Maverick total disc replacement (TDR) adjacent to a long spinal fusion.

METHODS

Seven spine specimens (T8-S1) were used in this study (66 ± 19 years old, 3F/4 M). A continuous pure moment of ±5.0 Nm was applied to the specimen in flexion-extension (FE), lateral bending (LB) and axial rotation (AR), with a compressive follower preload of 400 N. The 5.0 Nm data were analyzed to evaluate the operated segment biomechanics at the level of the disc replacements. The data were also analyzed at lower moments using a modified version of Panjabi's proposed "hybrid" method to evaluate adjacent segment kinematics (intervertebral motion at the segments adjacent to the fusion) under identical overall (T8-S1) specimen rotations. The motion of each vertebra was monitored with an optoelectronic camera system. The biomechanical test was completed for (1) the intact condition and repeated after each surgical technique was applied to the specimen, (2) capsulotomy at L4-L5 and L5-S1, (3) T8-L4 fusion and capsulotomy at L4-L5 and L5-S1, (4) Maverick at L4-L5, and (5) Maverick at L5-S1. The capsulotomy was performed to allow measurement of facet joint loads in a companion study. Paired t tests were used to determine if differences in the kinematic parameters measured were significant. Holm-Sidak corrections for multiple comparisons were applied where appropriate.

RESULTS

Under the 5.0 Nm loads, L4-L5 ROMs tended to decrease in all directions following L4-L5 Maverick replacement (mean = 22 %, compared to the fused condition). Two-level Maverick implantation also tended to reduce L4-S1 ROM (mean 18, 7 and 31 % in FE, LB and AR, respectively, compared to the fused condition without TDR). Following TDR replacement, the HAM location tended to shift posteriorly in FE (at L5-S1), anteriorly in AR, and inferiorly in LB. However, although the above-mentioned trends were observed, neither one- nor two-level TDR replacement showed statistically significant ROM or HAM change in any of the three directions. At the identical T8-S1 posture identified by the modified hybrid analysis, the L4-L5 and L5-S1 levels underwent significant larger motions, relative to the overall specimen rotation, after fusion. In the hybrid analysis, there were no significant differences between the ROM after fusion with intact natural discs at L4-L5 and L5-S1 and the motions at those levels with one or two TDRs implanted.

CONCLUSIONS

The present results demonstrated that one or two Maverick discs implanted subjacent to a long thoracolumbar fusion preserved considerable and intact-like ranges of motion and maintained motion patterns similar to the intact specimen, in this ex vivo study with applied pure moments and compressive follower preload. The hybrid analysis demonstrated that, after fusion, the TDR-implanted levels are required to undergo large rotations, relative to those necessary before fusion, in order to achieve the same motion between T8 and S1. Additional clinical and biomechanical research is necessary to determine if such a kinematic demand would be made on these levels clinically and the biomechanical performance of these implants if it were.

Biomechanical evaluation of the Total Facet Arthroplasty System® (TFAS®): loading as compared to a rigid posterior instrumentation system

PURPOSE

To gain insight into a new technology, a novel facet arthroplasty device (TFAS) was compared to a rigid posterior fixation system (UCR). The axial and bending loads through the implants and at the bone-implant interfaces were evaluated using an ex vivo biomechanical study and matched finite element analysis. Kinematic behaviour has been reported for TFAS, but implant loads have not. Implant loads are important indicators of an implant's performance and safety. The rigid posterior fixation system is used for comparison due to the extensive information available about these systems.

METHODS

Unconstrained pure moments were applied to 13 L3-S1 cadaveric spine segments. Specimens were tested intact, following decompression, UCR fixation and TFAS implantation at L4-L5. UCR fixation was via standard pedicle screws and TFAS implantation was via PMMA-cemented transpedicular stems. Three-dimensional 10 Nm moments and a 600 N follower load were applied; L4-L5 disc pressures and implant loads were measured using a pressure sensor and strain gauges, respectively. A finite element model was used to calculate TFAS bone-implant interface loads.

RESULTS

UCR experienced greater implant loads in extension (p < 0.004) and lateral bending (p < 0.02). Under flexion, TFAS was subject to greater implant moments (p < 0.04). At the bone-implant interface, flexion resulted in the smallest TFAS (average = 0.20 Nm) but greatest UCR (1.18 Nm) moment and axial rotation resulted in the greatest TFAS (3.10 Nm) and smallest UCR (0.40 Nm) moments. Disc pressures were similar to intact for TFAS but not for UCR (p < 0.04).

CONCLUSIONS

These results are most applicable to the immediate post-operative period prior to remodelling of the bone-implant interface since the UCR and TFAS implants are intended for different service lives (UCR--until fusion, TFAS--indefinitely). TFAS reproduced intact-like anterior column load-sharing--as measured by disc pressure. The highest bone-implant moment of 3.1 Nm was measured in TFAS and for the same loading condition the UCR interface moment was considerably lower (0.4 Nm). For other loading conditions, the differences between TFAS and UCR were smaller, with the UCR sometimes having larger values and for others the TFAS was larger. The long-term physiological meaning of these findings is unknown and demonstrates the need for a better understanding of the relationship between spinal arthroplasty devices and the host tissue as development of next generation motion-preserving posterior devices that hope to more accurately replicate the natural functions of the native tissue continues.

Spinal facet joint biomechanics and mechanotransduction in normal, injury and degenerative conditions

The facet joint is a crucial anatomic region of the spine owing to its biomechanical role in facilitating articulation of the vertebrae of the spinal column. It is a diarthrodial joint with opposing articular cartilage surfaces that provide a low friction environment and a ligamentous capsule that encloses the joint space. Together with the disc, the bilateral facet joints transfer loads and guide and constrain motions in the spine due to their geometry and mechanical function. Although a great deal of research has focused on defining the biomechanics of the spine and the form and function of the disc, the facet joint has only recently become the focus of experimental, computational and clinical studies. This mechanical behavior ensures the normal health and function of the spine during physiologic loading but can also lead to its dysfunction when the tissues of the facet joint are altered either by injury, degeneration or as a result of surgical modification of the spine. The anatomical, biomechanical and physiological characteristics of the facet joints in the cervical and lumbar spines have become the focus of increased attention recently with the advent of surgical procedures of the spine, such as disc repair and replacement, which may impact facet responses. Accordingly, this review summarizes the relevant anatomy and biomechanics of the facet joint and the individual tissues that comprise it. In order to better understand the physiological implications of tissue loading in all conditions, a review of mechanotransduction pathways in the cartilage, ligament and bone is also presented ranging from the tissue-level scale to cellular modifications. With this context, experimental studies are summarized as they relate to the most common modifications that alter the biomechanics and health of the spine-injury and degeneration. In addition, many computational and finite element models have been developed that enable more-detailed and specific investigations of the facet joint and its tissues than are provided by experimental approaches and also that expand their utility for the field of biomechanics. These are also reviewed to provide a more complete summary of the current knowledge of facet joint mechanics. Overall, the goal of this review is to present a comprehensive review of the breadth and depth of knowledge regarding the mechanical and adaptive responses of the facet joint and its tissues across a variety of relevant size scales.

Development of a versatile intra-articular pressure sensing array

A new sensor array intended to accurately and directly measure spatial and time-dependent pressures within a highly curved biological intra-articular joint was developed and tested. To evaluate performance of the new sensor array for application within intra-articular joints generally, and specifically to fit within the relatively restrictive space of the lumbar spine facet joint, geometric constraints of length, width, thickness and sensor spatial resolution were evaluated. Additionally, the effects of sensor array curvature, frequency response, linearity, drift, hysteresis, repeatability, and total system cost were assessed. The new sensor array was approximately 0.6mm in thickness, scalable to below the nominal 12 mm wide by 15 high lumbar spine facet joint size, offered no inherent limitations on the number or spacing of the sensors with less than 1.7% cross talk with sensor immediately adjacent to one another. No difference was observed in sensor performance down to a radius of curvature of 7 mm and a 0.66±0.97% change in sensor sensitivity was observed at a radius of 5.5mm. The sensor array had less than 0.07 dB signal loss up to 5.5 Hz, linearity was 0.58±0.13% full scale (FS), drift was less than 0.2% FS at 250 s and less than 0.6% FS at 700 s, hysteresis was 0.78±0.18%. Repeatability was excellent with a coefficient of variation less than 2% at pressures between 0 and 1.000 MPa. Total system cost was relatively small as standard commercially available data acquisition systems could be utilized, with no specialized software, and individual sensors within an array can be replaced as needed. The new sensor array had small and scalable geometry and very acceptable intrinsic performance including minimal to no alteration in performance at physiologically relevant ranges of joint curvature.

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.