Predictive structural models resulting in a trabecular bone topology closely resembling real bone would be a step toward 3D printing of sympathetic prosthetics. This study modifies an established trabecular bone structural adaptation approach, with the objective of achieving an improved adapted topology, specifically connectivity, compared to CT imaging studies; whilst retaining continuum level mechanical properties consistent with those reported in experimental studies. Strain driven structural adaptation models successfully identify trabecular trajectories, although tend to overpredict connectivity and skew trabecular radii distribution towards the smallest radius included in the adaptation. Radius adaptation of each trabecula is driven by a mechanostat approach with a target strain (1250 µɛ) below which radius is decreased (resorption), and above which radius is increased (apposition). Simulations include a lazy zone, in which neither resorption nor apposition takes place (1000 to 1500 µɛ); and a dead zone (<250 µɛ) in which complete resorption of trabeculae with the smallest included radius takes place. This study assesses the impact of increasing the dead zone threshold from <250 µɛ to <1000 µɛ, the lower limit of the lazy zone. In-silico structural models with an initial connectivity (number of trabeculae connecting at each joint) of 14 were generated using a nearest neighbour approach applied to a random cloud of points. Trabeculae were modelled using circular beams whose radii were adapted in response to normal strains caused by the axial force and bending moments due to a vertical pressure of 1 MPa applied to the top of the lattice, with the bottom of the lattice fixed in the vertical direction. Lattices in which nodes are either able (rigid jointed) or unable (pin jointed) to transmit bending moments were considered. Five virtual samples of each lattice type were used, and each simulation repeated twice: with a dead zone of either <250 µɛ or <1000 µɛ.Abstract
1.0 Objectives
2.0 Methods
A fibril reinforced multiphasic cartilage model was developed to improve the understanding of the depth-dependent cartilage internal structure and its through thickness biomechanical response. The heterogeneous model of cartilage was validated against full-field strain measurement obtained via Digital Image Correlation (DIC) during free swelling experiments. Hemi-cylindrical cartilage cores of 5mm diameter were obtained from porcine femoral condyles and humeral heads. The full field behaviour of these samples was monitored using DIC during an osmotic free swelling experiment performed following a standardised protocol [1]. Computational models were created in FEBio (version 2.8, Abstract
Objectives
Methods
Prediction tools are instruments which are commonly used to estimate the prognosis in oncology and facilitate clinical decision-making in a more personalized manner. Their popularity is shown by the increasing numbers of prediction tools, which have been described in the medical literature. Many of these tools have been shown to be useful in the field of soft-tissue sarcoma of the extremities (eSTS). In this annotation, we aim to provide an overview of the available prediction tools for eSTS, provide an approach for clinicians to evaluate the performance and usefulness of the available tools for their own patients, and discuss their possible applications in the management of patients with an eSTS. Cite this article:
The screw fastening torque applied during bone fracture fixation has a decisive influence on subsequent bone healing. Insufficient screw tightness can result in device/construct instability; conversely, excessive torques risk damaging the bone causing premature fixation failure. This effect is even more prominent in osteoporotic bone, a condition associated annually with almost 9 million fractures worldwide. During fracture fixation, screw tightening torque is applied using subjective feel. This approach may not be optimal for patient”s recovery, increasing risk of fixation failure, particularly in osteoporotic bone, and potentially require revision surgical interventions. Besides bone density, various factors influence the performance of screw fixation. These factors include bone geometry, cortical thickness and time-dependant relaxation behaviour of the bone. If the influence of screw fastening torque on the bone and relationships between these factors was better understood, the surgical technique could be optimised to reduce the risk of complications. Within this study, we developed an axisymmetric finite element (FE) model of bone screw tightening incorporating viscoelastic behaviour of the cortical bone such as creep and stress relaxation. The model anticipated time-dependent behaviour of the bone for different bone thickness and density after a typical bone fixation screw had been inserted. The idealised model has been developed based on CT scans of bones with varying densities and inserted screws. The model was validated through a series of experiments involving bovine tibiae (4–5 months) to evaluate the evolution of surface strains with time (Ncorr v1.2). Stress distribution was assessed in photoelastic experiments using acrylic analogues. Relaxation tests have been performed in aqueous environment for up to 48 hours to ensure the relaxation would be complete. The creep behaviour (maximum principal strain) was compared against computational predictions. Our early simulations predicted relaxation strains on the surface of the bone to be 1.1% within 24 hours comparing favourably to 1.3% measured experimentally. Stress distribution patterns were in agreement with photoelastic results. Using experimentally derived viscoelastic properties, the model has the potential to predict creep and stress relaxation patterns after screw insertion with different fastening torques for bones with varying density and geometry. We aim to develop this into a planning tool providing guidance to surgeons for optimal tightening when using screw fixation, particularly in reduced quality bone.
The results of the original mobile bearing Oxford unicompartmental knee replacement (UKR) in the lateral compartment have been disappointing because of high dislocation rates (11%). This original implant used a flat bearing articulation on the tibial tray. To address the issue of dislocation a new implant (domed tibia with biconcave bearing to increase entrapment) was introduced with a modified surgical technique. The aim of this study was to compare the risk of dislocation between a domed and flat lateral UKR. Separate geometric computer models of an Oxford mobile bearing lateral UKR were generated for the two types of articulation between the tibial component and the meniscal bearing: Flat-on-flat (flat) and Concave-on-convex (domed). Each type of mobile bearing was used to investigate three distinct dislocation modes observed clinically: lateral to medial dislocation, with the bearing resting on the tray wall (L-M-Wall); medial to lateral dislocation, out of the joint space (M-L); anterior to posterior dislocation, out of the joint space (A-P). A size C tray and a medium femoral component and bearing were used in all models. The femoral component, tibial tray and bearing were first aligned in a neutral position. For each dislocation the tibial tray was restrained in all degrees of freedom. The femoral component was restrained from moving in the anterior-posterior directions and in the medial-lateral directions. The femoral component was also restrained from rotating about the anterior-posterior, medial-lateral and superior-inferior directions. This meant that the femoral component was only able to move in the superior-inferior direction. Different bearing sizes were inserted into the model and the effect that moving the femoral component medially and laterally had on the amount of distraction required to cause bearing dislocation was investigated.Introduction
Methods
Patients with cancer and bone metastases can have an increased risk of fracturing their femur. Treatment is based on the impending fracture risk: patients with a high fracture risk are considered for prophylactic surgery, whereas low fracture risk patients are treated conservatively with radiotherapy to decrease pain. Current clinical guidelines suggest to determine fracture risk based on axial cortical involvement of the lesion on conventional radiographs, but that appears to be difficult. Therefore, we developed a patient-specific finite element (FE) computer model that has shown to be able to predict fracture risk in an experimental setting and in patients. The goal of this study was to determine whether patient-specific finite element (FE) computer models are better at predicting fracture risk for femoral bone metastases compared to clinical assessments based on axial cortical involvement on conventional radiographs, as described in current clinical guidelines. 45 patients (50 affected femurs) affected with predominantly lytic bone metastases who were treated with palliative radiotherapy for pain were included. CT scans were made and patients were followed for six months to determine whether or not they fractured their femur. Non-linear isotropic FE models were created with the patient-specific geometry and bone density obtained from the CT scans. Subsequently, an axial load was simulated on the models mimicking stance. Failure loads normalized for bodyweight (BW) were calculated for each femur. High and low fracture risks were determined using a failure load of 7.5 × BW as a threshold. Experienced assessors measured axial cortical involvement on conventional radiographs. Following clinical guidelines, patients with lesions larger than 30 mm were identified as having a high fracture risk. FE predictions were compared to clinical assessments by means of diagnostic accuracy values (sensitivity, specificity and positive (PPV) and negative predictive values (NPV)). Seven femurs (14%) fractured during follow-up. Median time to fracture was 8 weeks. FE models were better at predicting fracture risk in comparison to clinical assessments based on axial cortical involvement (sensitivity 100% vs. 86%, specificity 74% vs. 42%, PPV 39% vs. 19%, and NPV 100% vs. 95%, for the FE computer model vs. axial cortical involvement, respectively). We concluded that patient-specific FE computer models improve fracture risk predictions of femoral bone metastases in advanced cancer patients compared to clinical assessments based on axial cortical involvement, which is currently used in clinical guidelines. Therefore, we are initiating a pilot for clinical implementation of the FE model.
Total shoulder arthroplasty is a well-tested procedure that offers pain relief and restores the joint function. However, failure rate is still high, and glenoid loosening is pointed as the main reason in orthopedic registers. In order to understand the principles of failure, the principal strain distributions after implantation with Comprehensive® Total Shoulder System of Biomet® were experimental and numerically studied to predict bone behavior. Fourth generation composite left humerus and scapula from Sawbones® were used. These were implanted with Comprehensive® Total Shoulder System (Biomet®) with a modular Hybrid® glenoid base and Regenerex® glenoid and placed in situ by an experienced surgeon. The structures were placed in order to simulate 90º abduction, including principal muscular actions. Muscle forces used were as follows: Deltoideus 300N, Infraspinatus 120N, Supraspinatus 90N, Subscapularis 225N. All bone structures were modeled considering cortical and the trabecular bone of the scapula. The components of prosthesis were placed in the same positions than those in the in vitro models. Geometries were meshed with tetrahedral linear elements, with material properties as follows: Elastic modulus of cortical bone equal to 16 GPa, elastic modulus of trabecular bone equal to 0.155 GPa, polyethylene equal to 1GPa and titanium equal to 110 GPa. The assumed Poisson's ratio was 0.3 in all except for polyethylene where we assumed a value of 0.4. The prosthesis was considered as glued to the adjacent bone. The finite element model was composed of 336 024 elements. At the glenoid cavity, the major influence of the strain distributions was observed at the posterior-superior region, in both cortical and trabecular bone structures. The system presents critical region around holes of fixation in glenoid component. At the trabecular bone, the maximum principal strains at the posterior-superior region ranged from 2250 µε to 3000 µε. While at the cortical bone, the maximum principal strains were 300 µε to 400 µε. The results observed evidence some critical regions of concern and the effect of implant in the bone strains mainly at the posterior-superior region of the glenoid cavity is pronounced. This indicates that this region is more affected by the implant if bone remodeling is a concern and it is due to the strain-shielding effect, which has been connected with loosening of the glenoid component.
Polyethylene wear of joint replacements can cause severe clinical complications, including; osteolysis, implant loosening, inflammation and pain. Wear simulator testing is often used to assess new designs, but it is expensive and time consuming. It is possible to predict the volume of polyethylene implant wear from finite element models using a modification of Archard's classic wear law [1–2]. Typically, linear elastic isotropic, or elasto-plastic material models are used to represent the polyethylene. The purpose of this study was to investigate whether use of a viscoelastic material model would significantly alter the predicted volumetric wear of a mobile-bearing unicompartmental knee replacement. Tensile creep-recovery experiments were performed to characterise the creep and relaxation behaviour of the polyethylene (moulded GUR 4150 samples machined to 180×20×1 mm). Samples were loaded to 3 MPa stress in 4 minutes, and then held for 6 hours, the tensile stress was removed and samples were left to relax for 6 hours. The mechanical test data was used fit to a validated three–dimensional fractional Maxwell viscoelastic constitutive material model [3]. An explicit finite element model of a mobile–bearing unicompartmental knee replacement was created, which has been described previously [4]. The medial knee replacement was loaded to 1200 N over a period of 0.2 s. The bearing was meshed using quadratic tetrahedral elements (1.5 mm seeding size based on results of a mesh convergence study), and the femoral component was represented as an analytical rigid body. Wear predictions were made from the contact stress and sliding distance using Archard's law, as has been described in the literature [1–2]. A wear factor of 5.24×10−11 was used based upon the work by Netter et al. [2]. All models were created and solved using ABAQUS finite element software (version 6.14, Simulia, Dassault Systemes). The fractional viscoelastic material model predicted almost twice as much wear (0.119 mm3/million cycles) compared to the elasto-plastic model (0.069 mm3/million cycles). The higher wear prediction was due to both an increased sliding distance and higher contact pressures in the viscoelastic model. These preliminary findings indicate the simplified elasto-plastic polyethylene material representation can underestimate wear predictions from numerical simulations. Polyethylene is known to be a viscoelastic material which undergoes creep clinically, and it is not surprising that it is necessary to represent that viscoelastic behaviour to accurately predict implant wear. However, it does increase the complexity and run time of such computational studies, which may be prohibitive.
The Charnley Elite femoral component was first introduced in 1992 as a new design variant of the original Charnley femoral component (De Puy, Leeds, UK) with modified neck and stem geometry. The original component had undergone few changes in nearly forty years and has excellent long-term results. Early migration of the new stem design was determined by Roentgen Stereophotogrammetric Analysis (RSA)1. Rapid early migration of a component relative to the bone, measured by RSA, is predictive of subsequent aseptic loosening for a number of femoral stems. As there was rapid early migration and rotation of the Charnley Elite stem, we predicted that the long-term results would be poor. An outcome assessment is indicated as stems of this type are still being implanted. One hundred Charnley Elite stems, implanted in our centre between 1994 and 1997 were included in a prospective, cross-sectional follow-up study. Outcome measures include validated clinical scores (Charnley hip score, Harris hip score and Oxford hip score) and radiological scores (Gruen classification) as well as revision rates over the past 10 years. The clinical follow-up supports the RSA predictions of early failure of the Charnley Elite femoral stem.
To review the outcomes of patients undergoing manipulation under anaesthetic (MUA) after primary total knee arthroplasty (TKA) and predict those that may require such a procedure. We prospectively analysed all patients who required MUA post TKA performed by 2 surgeons using the same prosthesis from 2003 to 2008 and compared them to a control group of primary TKA matched for age, gender and surgeon. All patients in both groups had pre- and post-operative measurements of range of movement. In addition risk factors were identified including warfarin and statin use, diabetes and body mass index.Purpose
Methods
Early migration of the new stem design was determined by Roentgen Stereophotogrammetric Analysis (RSA). Rapid early migration of a component relative to the bone, measured by RSA, is predictive of subsequent aseptic loosening for a number of femoral stems. As there was rapid early migration and rotation of the Charnley Elite stem, we predicted that the long-term results would be poor. An outcome assessment is required as stems of this type are still being implanted.
Preliminary clinical scores in the patients who had not undergone any subsequent surgery were adequate (Oxford Hip Score mean average of 23.9). Thirteen percent of radiographs analysed had evidence of loosening, giving an overall loosening rate of 14% at 8 years.
The Stanford Upper Extremity Model (SUEM) (Holzbauer, Murray, Delp 2005, Ann Biomed Eng) includes the major muscles of the upper limb and has recently been described in scientific literature for various biomechanical purposes including modeling the muscle behavior after shoulder arthroplasty (Hoenecke, Flores-Hernandez, D'Lima 2014, J Shoulder Elbow Surg; Walker, Struk, Banks 2013, ISTA Proceedings). The initial publication of the SUEM compared the muscle moment arm predictions of the SUEM against various moment arm studies and all with the scapula fixed. A more recent study (Ackland, Pak, and Pandy 2008, J Anat) is now available that can be used to compare SUEM moment arm predictions to cadaver data for similar muscle sub-regions, during abduction and flexion motions, and with simulated scapular motion. SUEM muscle moment arm component vectors were calculated using the OpenSim Analyze Tool for an idealized abduction and an idealized flexion motion from 10° to 90° that corresponded to the motions described in Ackland for the cadaver arms. The normalized, averaged muscle moment arm data for the cadavers was manually digitized from the published figures and then resampled into uniform angles matching the SUEM data. Standard deviations of the muscle moment arms from the cadaver study were calculated from source data provided by the study authors. Python code was then used to calculate the differences, percent differences, and root-mean-square (RMS) values between the data sets. Of the 14 muscle groups in the SUEM, the smallest difference in predicted and measured moment arm was for the supraspinatus during the abduction task, with an RMS of the percent difference of 11.4%. In contrast, the middle latissimus dorsi had an RMS percent difference over 400% during the flexion task. The table presents the RMS difference and the RMS of the percent difference for the muscles with the largest abduction and adduction moment arms (during abduction) and the largest flexion and extension moment arms (during flexion). The moment arm data for the SUEM model and the cadaver data (with 1 standard deviation band) during the motion of the same muscles are provided in Figure 1 for the Abduction motion task and in Figure 2 for the Flexion motion task. It is challenging to simulate the three dimensional, time variant geometries of shoulder muscles while maintaining model fidelity and optimizing computational cost. Dividing muscles in to sub regions and using wrapping line segment approximations appears a reasonable strategy though more work could improve model accuracy especially during complex three dimensional motions.
FDA approval of metal-on-metal (MOM: 28, 32mm) bearings has provided 10 years of clinical experiences in USA. However there has been no detailed mapping of wear phenomena in retrieval cases. We present an analysis of 28 cases, MOM retrievals with 1 to 10 years follow-up, radiographic reviews and metal ion studies. Ball diameters ranged from 28mm to 42mm. Two balls were the early design with skirts. Main indicators for revision were the progressive radiographic changes indicative of osteolysis, with associated hip pain. Approximately 54% of patients were males and ages ranged from 36 to 76 years of age. Only 7 femoral stems were recovered but all had impingement marks. Only three cases lacked any evidence of stripe wear and these were in very elderly patients. Approximately 85% of these cases showed some evidence of stripe wear and multiple stripes were clearly visible on 50% of the femoral balls. The medial ball stripes were twice as common as lateral. Stripe wear was identified in 25% of CoCr liners. In the hip simulator studies generally show ‘run-in’ wear rates of 1–7mm3 per million cycles (Mc). We noted that above the 5mm3/Mc threshold, the serum generally appeared black. In contrast, the ‘steady-state’ wear rates of 0.1–1.6 mm3/Mc showed the true potential of MOM bearings. However there were often examples of higher wear (7–20 mm3/Mc), which gave confounding trends in published studies. Our studies of metal ions in the simulator lubricant provided a very accurate representation of MOM wear. There are many limitations in comparing in-vitro to in-vivo wear performance. Our retrieval data are biased to cases that failed due to hip pain, had radiographic signs of progressive osteolysis and some showed high levels of metal ions. There was also the bias of having predominantly a CoCr sandwich design (polyethylene adaptor). Use of the small ball added the well-known risks of impingement, subluxation and dislocation with rigid cups. Using the ‘damage modes’ from McKellop, we found only normal Mode-1 wear to be rare in these cases, whereas Modes# 2–4 had an incidence approaching 30% each. Signs of impingement were evident in 85% of our cases. Thus summarizing these MOM wear phenomena in retrieved 28mm sandwich cups, the evidence implicated impingement and 3rd-body wear modes (#2–4) as the clinical risk for adverse wear effects at 10 years follow-up. The in-vitro wear studies have not yet simulated such adverse clinical effects.
Many finite element (FE) studies have been performed in the past to assess the biomechanical performance of TKA and THA components. The boundary conditions have often been simplified to a few peak loads. With the availability of personalized musculoskeletal (MS) models we becomes possible to estimate dynamic muscle and prosthetic forces in a patient specific manner. By combining this knowledge with FE models, truly patient specific failure analyses can be performed. In this study we applied this combined technique to the femoral part of a cementless THR and calculated the cyclic micro-motions of the stem relative to the bone in order to assess the potential for bone ingrowth. An FE model of a complete femur with a CLS Spotorno stem inserted was generated. An ideal fit between the implant and the bone was modeled proximally, whereas distally an interface gap of 100μm was created to simulate a more realistic interface condition obtained during surgery. Furthermore, a gait analysis was performed on a young subject and fed into the Anybody™ MS modeling system. The anatomical data set (muscle attachment points) used by the Anybody™ system was morphed to the shape of the femoral reconstruction. In this way a set of muscle attachment points was obtained which was consistent with the FE model. The predicted muscle and hip contact forces by the Anybody™ modeling system were dynamic and divided into 37 increments including two stance phases and a swing phase of the right leg.Introduction
Methods
Summary Statement. A coupled finite element - analytical model is presented to predict and to elucidate a clinical healing scenario where bone regenerates in a critical-sized femoral defect, bounded by periosteum or a periosteum substitute implant and stabilised via an intramedullary nail. Introduction. Bone regeneration and maintenance processes are intrinsically linked to mechanical environment. However, the cellular and subcellular mechanisms of mechanically-modulated bone (re-) generation are not fully understood. Recent studies with periosteum osteoprogenitor cells exhibit their mechanosensitivity in vitro and in situ. In addtion, while a variety of growth factors are implicated in bone healing processes, bone morphogenetic protein-2 (BMP-2) is recognised to be involved in all stages of bone regeneration. Furthermore, periosteal injuries heal predominantly via endochondral ossification mechanisms. With this background in mind, the current study aims to understand the role of mechanical environment on BMP-2 production and periosteally-mediated bone regeneration. The one-stage bone transport model [1] provides a clinically relevant experimental platform on which to model the mechanobiological process of periosteum-mediated bone regeneration in a critical-sized defect. Here we develop a model framework to study the cellular-, extracellular- and mechanically-modulated process of defect infilling, governed by the mechanically-modulated production of BMP-2 by osteoprogenitor cells located in the periosteum. Methods. Material properties of the healing callus and periosteum contribute to the strain stimulus sensed by osteoprogenitor cells therein. Using a mechanical finite element model, periosteal surface strains are first predicted as a function of callus properties. Strains are then input to a mechanistic mathematical model, where mechanical regulation of BMP-2 production mediates rates of cellular proliferation, differentiation and extracellular matrix (ECM) production, to predict healing outcomes. A parametric approach enables the spatial and temporal prediction of tissue regeneration via endochondral ossification.
Introduction. Gait laboratory measurement of whole-body kinematics and ground reaction forces during a wide range of activities is frequently performed in joint replacement patient diagnosis, monitoring, and rehabilitation programs. These data are commonly processed in musculoskeletal modeling platforms such as OpenSim and Anybody to estimate muscle and joint reaction forces during activity. However, the processing required to obtain musculoskeletal estimates can be time consuming, requires significant expertise, and thus seriously limits the patient populations studied. Accordingly, the purpose of this study was to evaluate the potential of deep learning methods for estimating muscle and joint reaction forces over time given kinematic data, height, weight, and ground reaction forces for total knee replacement (TKR) patients performing activities of daily living (ADLs). Methods. 70 TKR patients were fitted with 32 reflective markers used to define anatomical landmarks for 3D motion capture. Patients were instructed to perform a range of tasks including gait, step-down and sit-to-stand. Gait was performed at a self-selected pace, step down from an 8” step height, and sit-to-stand using a chair height of 17”. Tasks were performed over a force platform while force data was collected at 2000 Hz and a 14 camera motion capture system collected at 100 Hz. The resulting data was processed in OpenSim to estimate joint reaction and muscle forces in the hip and knee using static optimization. The full set of data consisted of 135 instances from 70 patients with 63 sit-to-stands, 15 right-sided step downs, 14 left-sided step downs, and 43 gait sequences. Two classes of neural networks (NNs), a recurrent neural network (RNN) and temporal convolutional neural network (TCN), were trained to predict activity classification from joint angle, ground reaction force, and anthropometrics. The NNs were trained to predict muscle and joint reaction forces over time from the same input metrics. The 135 instances were split into 100 instances for training, 15 for validation, and 20 for testing. Results. The RNN and TCN yielded classification accuracies of 90% and 100% on the test set. Correlation coefficients between ground truth and predictions from the test set ranged from 0.81–0.95 for the RNN, depending on the activity.
Objectives. Opening wedge high tibial osteotomy (HTO) is an established surgical procedure for the treatment of early-stage knee arthritis. Other than infection, the majority of complications are related to mechanical factors – in particular, stimulation of healing at the osteotomy site. This study used finite element (FE) analysis to investigate the effect of plate design and bridging span on interfragmentary movement (IFM) and the influence of fracture healing on plate stress and potential failure. Materials and Methods. A 10° opening wedge HTO was created in a composite tibia. Imaging and strain gauge data were used to create and validate FE models. Models of an intact tibia and a tibia implanted with a custom HTO plate using two different bridging spans were validated against experimental data. Physiological muscle forces and different stages of osteotomy gap healing simulating up to six weeks postoperatively were then incorporated.
Disorders of bone integrity carry a high global disease burden, frequently requiring intervention, but there is a paucity of methods capable of noninvasive real-time assessment. Here we show that miniaturized handheld near-infrared spectroscopy (NIRS) scans, operated via a smartphone, can assess structural human bone properties in under three seconds. A hand-held NIR spectrometer was used to scan bone samples from 20 patients and predict: bone volume fraction (BV/TV); and trabecular (Tb) and cortical (Ct) thickness (Th), porosity (Po), and spacing (Sp).Aims
Methods
Introduction: Vertebral body osteophytes are common in elderly spines, but their mechanical function is unclear. Do they act primarily to reduce compressive stress on the vertebral body, or to stabilise the spine in bending?. Methods: Spines were obtained from cadavers aged 51–92yrs (mean 77yrs) with radiographic evidence of vertebral osteophytes (mostly antero-lateral). Twenty motion segments, from T5-T6 to L3-L4, were dissected and loaded a) in compression to 1.5kN, and b) in bending to 10–25Nm. Vertebral movements were tracked at 50Hz using an optical MacReflex system. Bending tests were performed in random order, in flexion, extension, and lateral bending. Resistance to bending and compression was measured before and after surgical excision of all osteophytes. Bone mineral content (BMC) of osteophytes was measured using DXA. ANOVA was used to detect changes after osteophyte excision, and regression was used to examine the influence of osteophyte size. Results: Compressive stiffness was reduced by an average 17% following osteophyte removal (p<
0.05). In flexion and extension, bending stiffness was reduced by 60% and 79% respectively (p<
0.01). In left and right lateral bending, stiffness was reduced by 42% and 49% respectively. Osteophyte removal increased the neutral zone and range of motion in each mode of bending, and changes were proportional to osteophyte mass and BMC (p<
0.01). Conclusion: Vertebral body osteophytes primarily stabilise the spine in bending, and do little to resist compression, despite their considerable BMC.