Computational wear simulation of patellofemoral articular cartilage during in vitro testing

Though changes in normal joint motions and loads (e.g., following anterior cruciate ligament injury) contribute to the development of knee osteoarthritis, the precise mechanism by which these changes induce osteoarthritis remains unknown. As a first step toward identifying this mechanism, this study evaluates computational wear simulations of a patellofemoral joint specimen wear tested on a knee simulator machine. A multibody dynamic model of the specimen mounted in the simulator machine was constructed in commercial computer-aided engineering software. A custom elastic foundation contact model was used to calculate contact pressures and wear on the femoral and patellar articular surfaces using geometry created from laser scan and MR data. Two different wear simulation approaches were investigated--one that wore the surface geometries gradually over a sequence of 10 one-cycle dynamic simulations (termed the "progressive" approach), and one that wore the surface geometries abruptly using results from a single one-cycle dynamic simulation (termed the "non-progressive" approach). The progressive approach with laser scan geometry reproduced the experimentally measured wear depths and areas for both the femur and patella. The less costly non-progressive approach predicted deeper wear depths, especially on the patella, but had little influence on predicted wear areas. Use of MR data for creating the articular and subchondral bone geometry altered wear depth and area predictions by at most 13%. These results suggest that MR-derived geometry may be sufficient for simulating articular cartilage wear in vivo and that a progressive simulation approach may be needed for the patella and tibia since both remain in continuous contact with the femur.     

Potential for thermal damage to articular cartilage by PMMA reconstruction of a bone cavity following tumor excision: A finite element study

Benign, giant cell tumors are often treated by intralesional excision and reconstruction with polymethylmethacrylate (PMMA) bone cement. The exothermic reaction of the in-situ polymerizing PMMA is believed to beneficially kill remaining tumor cells. However, at issue is the extent of this necrotic effect into the surrounding normal bone and the adjacent articular cartilage. Finite element analysis (ABAQUS 6.4-1) was used to determine the extent of possible thermal necrosis around prismatically shaped, PMMA implants (8–24 cc in volume), placed into a peripheral, sagittally symmetric, metaphyseal defect in the proximal tibia. Temperature/exposure time conditions indicating necrotic potential during the exotherm of the polymerizing bone cement were found in regions of the cancellous bone within 3 mm of the superior surface of the PMMA implant. If less than 3 mm of cancellous bone existed between the PMMA implant and the subchondral bone layer, regions of the subchondral bone were also exposed to thermally necrotic conditions. However, as long as there were at least 2 mm of uniform subchondral bone above the PMMA implant, the necrotic regions did not extend into the overlying articular cartilage. This was the case even when the PMMA was in direct contact with the subchondral bone. If the subchondral bone is not of sufficient thickness, or is not continuous, then care should be taken to protect the articular cartilage from thermal damage as a result of the reconstruction of the tumor cavity with PMMA bone cement.     

Double-step registration of in vivo stereophotogrammetry with both in vitro 6-DOFs electrogoniometry and CT medical imaging

Standard registration techniques of bone morphology to motion analysis data often lead to unsatisfactory motion simulation because of discrepancies during the location of anatomical landmarks in the datasets. This paper describes an iterative registration method of a three-dimensional (3D) skeletal model with both 6 degrees-of-freedom joint kinematics and standard motion analysis data. The method is demonstrated in this paper on the lower limb. The method includes two steps. A primary registration allowed synchronization of in vitro kinematics of the knee and ankle joints using flexion/extension angles from in vivo gait analysis. Results from primary registration were then improved by a so-called advanced registration, which integrated external constraints obtained from experimental gait pre-knowledge. One cadaver specimen was analyzed to obtain both joint kinematics of knee and ankle joints using 3D electrogoniometry, and 3D bone morphology from medical imaging data. These data were registered with motion analysis data from a volunteer during the execution of locomotor tasks. Computer graphics output was implemented to visualize the results for a motion of sitting on a chair. Final registration results allowed the observation of both in vivo motion data and joint kinematics from the synchronized specimen data. The method improved interpretation of gait analysis data, thanks to the combination of realistic 3D bone models and joint mechanism. This method should be of interest both for research in gait analysis and medical education. Validation of the overall method was performed using RMS of the differences between bone poses estimated after registration and original data from motion analysis.     

Three-dimensional dynamic hip contact area and pressure distribution during activities of daily living

Estimation of the hip joint contact area and pressure distribution during activities of daily living is important in predicting joint degeneration mechanism, prosthetic implant wear, providing biomechanical rationales for preoperative planning and postoperative rehabilitation. These biomechanical data were estimated utilizing a generic hip model, the Discrete Element Analysis technique, and the in vivo hip joint contact force data. The three-dimensional joint potential contact area was obtained from the anteroposterior radiograph of a subject and the actual joint contact area and pressure distribution in eight activities of daily living were calculated. During fast, normal, and slow walking, the peak pressure of moderate magnitude was located at the lateral roof of the acetabulum during mid-stance. In standing up and sitting down, and during knee bending, the peak pressures were located at the edge of the posterior horn and the magnitude of the peak pressure during sitting down was 2.8 times that of normal walking. The peak pressure was found at the lateral roof in climbing up stairs which was higher than that in going down stairs. These results can be used to rationalize rehabilitation protocols, functional restrictions after complex acetabular reconstructions, and prosthetic component wear and fatigue test set up. The same model and analysis can provide further insight to soft tissue loading and pathology such as labral injury. When the pressure distribution on the acetabulum is inverted onto the femoral head, prediction of subchondral bone collapse associated with avascular necrosis can be achieved with improved accuracy.     

Analysis of surface-to-surface distance mapping during three-dimensional motion at the ankle and subtalar joints

Joint surface interaction and ligament constraints determine the kinematic characteristics of the ankle and subtalar joints. Joint surface interaction is characterized by joint contact mechanics and by relative joint surface position potentially characterized by distance mapping. While ankle contact mechanics was investigated, limited information is available on joint distance mapping and its changes during motion. The purpose of this study was to use image-based distance mapping to quantify this interaction at the ankle and subtalar joints during tri-planar rotations of the ankle complex. Five cadaveric legs were scanned using Computed Tomography and the images were processed to produce 3D bone models of the tibia, fibula, talus and calcaneus. Each leg was tested on a special linkage through which the ankle complex was loaded in dorsiflexion/plantarflexion, inversion/eversion, and internal/external rotation and the resulting bone movements were recorded. Fiduciary bone markers data and 3D bone models were combined to generate color-coded distance maps for the ankle and subtalar joints. The results were processed focusing on the changes in surface-to-surface distance maps between the extremes of the range of motion and neutral. The results provided detailed insight into the three-dimensional highly coupled nature of these joints showing significant and unique changes in distance mapping from neutral to extremes of the range of motion. The non-invasive nature of the image-based distance mapping technique could result, after proper modifications, in an effective diagnostic and clinical evaluation technique for application such as ligament injuries and quantifying the effect of arthrodesis or total ankle replacement surgery.     

Trabecular bone strain changes associated with subchondral stiffening of the proximal tibia

Subchondral stiffening is a hallmark pathologic feature of osteoarthritis but its mechanical and temporal relationship to the initiation or the progression of osteoarthritis is not established. The mechanical effect of subchondral stiffening on the surrounding trabecular bone is poorly understood. This study employs a relatively new application of digital image correlation to measure strain in the trabecular region of the proximal medial tibia in normal specimens and in specimens with simulated subchondral bone stiffening. Coronal sections from eight normal human cadaveric proximal tibiae were loaded in static compression and high resolution contact radiographs were made. Repeat contact radiographs were collected after the subchondral bone near the jointline was stiffened using polymethylmethacrylate. Digital images, made from loaded and unloaded contact radiographs, were compared using in-house software to measure trabecular displacement and calculate trabecular bone strain. Overall strain was higher in the stiffened specimens suggesting experimental artifiact significantly affected our results. Consistent increases in median maximum shear strain, median maximum principal strain, median minimum principal strain, and peak shear strain were measured near the inner and outer edges of the stiffened segment. Our experiment provides direct experimental measurement of increases in trabecular bone strain caused by subchondral stiffening, however, the clinical and biologic importance of strain increases is unknown.     

Automated measurement of neural foramen cross-sectional area during in vivo functional movement

An automated technique to measure neural foramen cross-sectional area during in vivo, multi-planar movements is presented. This method combines three-dimensional (3D) models of each vertebra obtained from CT scans with in vivo movement data collected using high-speed biplane radiography. A novel computer algorithm that automatically traces a path around the bony boundary that defines the neural foramen at every frame of X-ray data is described. After identifying the neural foramen boundary, the cross-sectional area is calculated. The technique is demonstrated using data collected from a patient with cervical radiculopathy who is tested before and after conservative treatment. The technique presented here can be applied when 3D, dynamic, functional movements are performed. Neural foramen cross-sectional area can be quantified at specific angles of intervertebral rotation, allowing for matched comparisons between two trials or two test sessions. The present technique is ideal for longitudinal studies involving subjects who receive conservative or surgical treatments that may affect spine motion.     

Static and dynamic error of a biplanar videoradiography system using marker-based and markerless tracking techniques

The use of biplanar videoradiography technology has become increasingly popular for evaluating joint function in vivo. Two fundamentally different methods are currently employed to reconstruct 3D bone motions captured using this technology. Marker-based tracking requires at least three radio-opaque markers to be implanted in the bone of interest. Markerless tracking makes use of algorithms designed to match 3D bone shapes to biplanar videoradiography data. In order to reliably quantify in vivo bone motion, the systematic error of these tracking techniques should be evaluated. Herein, we present new markerless tracking software that makes use of modern GPU technology, describe a versatile method for quantifying the systematic error of a biplanar videoradiography motion capture system using independent gold standard instrumentation, and evaluate the systematic error of the W.M. Keck XROMM Facility's biplanar videoradiography system using both marker-based and markerless tracking algorithms under static and dynamic motion conditions. A polycarbonate flag embedded with 12 radio-opaque markers was used to evaluate the systematic error of the marker-based tracking algorithm. Three human cadaveric bones (distal femur, distal radius, and distal ulna) were used to evaluate the systematic error of the markerless tracking algorithm. The systematic error was evaluated by comparing motions to independent gold standard instrumentation. Static motions were compared to high accuracy linear and rotary stages while dynamic motions were compared to a high accuracy angular displacement transducer. Marker-based tracking was shown to effectively track motion to within 0.1?mm and 0.1 deg under static and dynamic conditions. Furthermore, the presented results indicate that markerless tracking can be used to effectively track rapid bone motions to within 0.15 deg for the distal aspects of the femur, radius, and ulna. Both marker-based and markerless tracking techniques were in excellent agreement with the gold standard instrumentation for both static and dynamic testing protocols. Future research will employ these techniques to quantify in vivo joint motion for high-speed upper and lower extremity impacts such as jumping, landing, and hammering.     

Computed Tomographic Imaging of Subchondral Fatigue Cracks in the Distal End of the Third Metacarpal Bone in the Thoroughbred Racehorse Can Predict Crack Micromotion in an Ex-Vivo Model

Articular stress fracture arising from the distal end of the third metacarpal bone (MC3) is a common serious injury in Thoroughbred racehorses. Currently, there is no method for predicting fracture risk clinically. We describe an ex-vivo biomechanical model in which we measured subchondral crack micromotion under compressive loading that modeled high speed running. Using this model, we determined the relationship between subchondral crack dimensions measured using computed tomography (CT) and crack micromotion. Thoracic limbs from 40 Thoroughbred racehorses that had sustained a catastrophic injury were studied. Limbs were radiographed and examined using CT. Parasagittal subchondral fatigue crack dimensions were measured on CT images using image analysis software. MC3 bones with fatigue cracks were tested using five cycles of compressive loading at -7,500N (38 condyles, 18 horses). Crack motion was recorded using an extensometer. Mechanical testing was validated using bones with 3 mm and 5 mm deep parasagittal subchondral slots that modeled naturally occurring fatigue cracks. After testing, subchondral crack density was determined histologically. Creation of parasagittal subchondral slots induced significant micromotion during loading (p<0.001). In our biomechanical model, we found a significant positive correlation between extensometer micromotion and parasagittal crack area derived from reconstructed CT images (S_R = 0.32, p<0.05). Correlations with transverse and frontal plane crack lengths were not significant. Histologic fatigue damage was not significantly correlated with crack dimensions determined by CT or extensometer micromotion. Bones with parasagittal crack area measurements above 30 mm² may have a high risk of crack propagation and condylar fracture in vivo because of crack micromotion. In conclusion, our results suggest that CT could be used to quantify subchondral fatigue crack dimensions in racing Thoroughbred horses in-vivo to assess risk of condylar fracture. Horses with parasagittal crack arrays that exceed 30 mm² may have a high risk for development of condylar fracture.     

Quantitative Assessment of Murine Articular Cartilage and Bone Using X-Ray Phase-Contrast Imaging

Murine models for rheumatoid arthritis (RA) research can provide important insights for understanding RA pathogenesis and evaluating the efficacy of novel treatments. However, simultaneously imaging both murine articular cartilage and subchondral bone using conventional techniques is challenging because of low spatial resolution and poor soft tissue contrast. X-ray phase-contrast imaging (XPCI) is a new technique that offers high spatial resolution for the visualisation of cartilage and skeletal tissues. The purpose of this study was to utilise XPCI to observe articular cartilage and subchondral bone in a collagen-induced arthritis (CIA) murine model and quantitatively assess changes in the joint microstructure. XPCI was performed on the two treatment groups (the control group and CIA group, n = 9 per group) to monitor the progression of damage to the femur from the knee joint in a longitudinal study (at 0, 4 and 8 weeks after primary injection). For quantitative assessment, morphologic parameters were measured in three-dimensional (3D) images using appropriate image analysis software. Our results showed that the average femoral cartilage volume, surface area and thickness were significantly decreased (P<0.05) in the CIA group compared to the control group. Meanwhile, these decreases were accompanied by obvious destruction of the surface of subchondral bone and a loss of trabecular bone in the CIA group. This study confirms that XPCI technology has the ability to qualitatively and quantitatively evaluate microstructural changes in mouse joints. This technique has the potential to become a routine analysis method for accurately monitoring joint damage and comprehensively assessing treatment efficacy.     

Apparent density of the primate calcaneo‐cuboid joint and its association with locomotor mode,foot posture,and the “midtarsal break”

Primates use a range of locomotor modes during which they incorporate various foot postures. Humans are unique compared with other primates in that humans lack a mobile fore‐ and midfoot. Rigidity in the human foot is often attributed to increased propulsive and stability requirements during bipedalism. Conversely, fore‐ and midfoot mobility in nonhuman primates facilitates locomotion in arboreal settings. Here, we evaluated apparent density (AD) in the subchondral bone of human, ape, and monkey calcanei exhibiting different types of foot loading. We used computed tomography osteoabsorptiometry and maximum intensity projection (MIP) maps to visualize AD in subchondral bone at the cuboid articular surface of calcanei. MIPs represent 3D volumes (of subchondral bone) condensed into 2D images by extracting AD maxima from columns of voxels comprising the volumes. False‐color maps are assigned to MIPs by binning pixels in the 2D images according to brightness values. We compared quantities and distributions of AD pixels in the highest bin to test predictions relating AD patterns to habitual locomotor modes and foot posture categories of humans and several nonhuman primates. Nonhuman primates exhibit dorsally positioned high AD concentrations, where maximum compressive loading between the calcaneus and cuboid likely occurs during “midtarsal break” of support. Humans exhibit less widespread areas of high AD, which could reflect reduced fore‐ and midfoot mobility. Analysis of the internal morphology of the tarsus, such as subchondral bone AD, potentially offers new insights for evaluating primate foot function during locomotion. Am J Phys Anthropol, 2010. © 2009 Wiley‐Liss, Inc.     

A methodology to accurately quantify patellofemoral cartilage contact kinematics by combining 3D image shape registration and cine-PC MRI velocity data

Patellofemoral osteoarthritis and its potential precursor patellofemoral pain syndrome (PFPS) are common, costly, and debilitating diseases. PFPS has been shown to be associated with altered patellofemoral joint mechanics; however, an actual variation in joint contact stresses has not been established due to challenges in accurately quantifying in vivo contact kinematics (area and location). This study developed and validated a method for tracking dynamic, in vivo cartilage contact kinematics by combining three magnetic resonance imaging (MRI) techniques, cine-phase contrast (CPC), multi-plane cine (MPC), and 3D high-resolution static imaging. CPC and MPC data were acquired from 12 healthy volunteers while they actively extended/flexed their knee within the MRI scanner. Since no gold standard exists for the quantification of in vivo dynamic cartilage contact kinematics, the accuracy of tracking a single point (patellar origin relative to the femur) represented the accuracy of tracking the kinematics of an entire surface. The accuracy was determined by the average absolute error between the PF kinematics derived through registration of MPC images to a static model and those derived through integration of the CPC velocity data. The accuracy ranged from 0.47 mm to 0.77 mm for the patella and femur and from 0.68 mm to 0.86 mm for the patellofemoral joint. For purely quantifying joint kinematics, CPC remains an analytically simpler and more accurate (accuracy <0.33 mm) technique. However, for application requiring the tracking of an entire surface, such as quantifying cartilage contact kinematics, this combined imaging approach produces accurate results with minimal operator intervention.     

Three-dimensional knee model: constrained by isometric ligament bundles and experimentally obtained tibio-femoral contacts

The purpose of this study is to investigate the effect of anterior portion of anterior cruciate ligament, posterior cruciate ligament, anterior and deep portions of medial collateral ligament and the tibio-femoral articular contacts on passive knee motion. A well-accepted reference model for a normal tibio-femoral joint is reconstructed from the literature. The proposed three-dimensional dynamic tibio-femoral model includes the isometric fascicles, ligament bundles and irregularly shaped medial-lateral contact surfaces. With the approach we aim to analyze bone shape and ligament related abnormalities of knee kinematics. The rotations, translations and the contact forces during passive knee flexion were compared against a reference model and the results were found in close accordance. This study demonstrated that isometric ligament bundles play an important role in understanding the femur shape from contact points on tibia. Femoral condyles are not necessarily spherical. The surgical treatments should consider both ligament bundle lengths and contact surface geometries to achieve a problem free knee kinematics after a knee surgery.     

Analysis of skin movement with respect to flexional bone motion using MR images of a hand

This study describes the methodology for analyzing skin movement with respect to the predefined flexional motion of the hand. The 3D in vivo geometric data were acquired through magnetic resonance (MR) scanning of multiple hands in various postures. Custom software was developed to segment MR images and to generate polygonal isosurfaces. A new bone coordinate system (BCS) was defined to describe the skin movement measured with the skin-surface markers. The BCS was defined by the surface registration technique, which could minimize the distance among the 3D polygonal surfaces of different postures. MR images of six subjects in the 3 postures were captured with skin-surface markers for the experiment. Skin movement on the second metacarpal bone, which was fixed during the flexional motion, ranged from 1.169 to 10.885 mm. While the distance of the skin movement varied from one subject to another, a common skin movement direction was observed among the six subjects. For the repeatability evaluation of the proposed method, the MR images of the same subject were independently processed 10 times by a trained operator. The standard deviation of the skin movement on the second metacarpal of the index finger with respect to the predefined flexional motion ranged from 0.45 to 0.74 mm.     

Three-dimensional modeling and animation of two carpal bones: a technique

The objectives of this study were to (a). create 3D reconstructions of two carpal bones from single CT data sets and animate these bones with experimental in vitro motion data collected during dynamic loading of the wrist joint, (b). develop a technique to calculate the minimum interbone distance between the two carpal bones, and (c). validate the interbone distance calculation process. This method utilized commercial software to create the animations and an in-house program to interface with three-dimensional CAD software to calculate the minimum distance between the irregular geometries of the bones. This interbone minimum distance provides quantitative information regarding the motion of the bones studied and may help to understand and quantify the effects of ligamentous injury.     

In vivo contact kinematics and contact forces of the knee after total knee arthroplasty during dynamic weight-bearing activities

Analysis of polyethylene component wear and implant loosening in total knee arthroplasty (TKA) requires precise knowledge of in vivo articular motion and loading conditions. This study presents a simultaneous in vivo measurement of tibiofemoral articular contact forces and contact kinematics in three TKA patients. These measurements were accomplished via a dual fluoroscopic imaging system and instrumented tibial implants, during dynamic single leg lunge and chair rising-sitting. The measured forces and contact locations were also used to determine mediolateral distribution of axial contact forces. Contact kinematics data showed a medial pivot during flexion of the knee, for all patients in the study. Average axial forces were higher for lunge compared to chair rising-sitting (224% vs. 187% body weight). In this study, we measured peak anteroposterior and mediolateral forces averaging 13.3% BW during lunge and 18.5% BW during chair rising-sitting. Mediolateral distributions of axial contact force were both patient and activity specific. All patients showed equitable medial-lateral loading during lunge but greater loads at the lateral compartment during chair rising-sitting. The results of this study may enable more accurate reproduction of in vivo loads and articular motion patterns in wear simulators and finite element models. This in turn may help advance our understanding of factors limiting longevity of TKA implants, such as aseptic loosening and polyethylene component wear, and enable improved TKA designs.     

An improved design of electrodes for measurement of streaming potentials on wet bone in vitro and in vivo

Streaming potentials are generated by mechanical stress in wet bone and may constitute a control mechanism for bone remodeling. Measurement of streaming potentials in bone has attracted considerable effort in past years but quantitative studies have been hampered by relatively poor repeatability when using Ag.AgCl electrodes which contact bone via a wick moistened with electrolyte. Improvement now has been achieved with an electrode design that limits the specific area of contact of an agar/salt bridge by means of a silastic seal, thus permitting the same equipotential surface to be contacted for each set of measurements. This reduces variations caused by bone structure and impedance, and facilitates quantitative comparisons of the response of bone samples to selected variables. The new design also permits considerable qualitative improvement in recordings made from bone during locomotor function in experimental animals in vivo.     

Measurement of in vivo anterior cruciate ligament strain during dynamic jump landing

Despite recent attention in the literature, anterior cruciate ligament (ACL) injury mechanisms are controversial and incidence rates remain high. One explanation is limited data on in vivo ACL strain during high-risk, dynamic movements. The objective of this study was to quantify ACL strain during jump landing. Marker-based motion analysis techniques were integrated with fluoroscopic and magnetic resonance (MR) imaging techniques to measure dynamic ACL strain non-invasively. First, eight subjects' knees were imaged using MR. From these images, the cortical bone and ACL attachment sites of the tibia and femur were outlined to create 3D models. Subjects underwent motion analysis while jump landing using reflective markers placed directly on the skin around the knee. Next, biplanar fluoroscopic images were taken with the markers in place so that the relative positions of each marker to the underlying bone could be quantified. Numerical optimization allowed jumping kinematics to be superimposed on the knee model, thus reproducing the dynamic in vivo joint motion. ACL length, knee flexion, and ground reaction force were measured. During jump landing, average ACL strain peaked 55±14 ms (mean and 95% confidence interval) prior to ground impact, when knee flexion angles were lowest. The peak ACL strain, measured relative to its length during MR imaging, was 12±7%. The observed trends were consistent with previously described neuromuscular patterns. Unrestricted by field of view or low sampling rate, this novel approach provides a means to measure kinematic patterns that elevate ACL strains and that provide new insights into ACL injury mechanisms.     

Overload arthrosis: strain patterns in the equine metacarpal condyle

An overload arthrosis occurs consistently in the palmar region of the metacarpal condyle of the equine fetlock (metacarpophalangeal) joint characterized by subchondral bone sclerosis, devitalization and mechanical failure leading to collapse of the overlying articular cartilage. Samples were selected of joints with mild, moderate, and severe subchondral sclerosis, in which cartilage collapse had not yet occurred. An additional group that had severe sclerosis with focal rarefaction suggesting impending collapse was also studied (n=5/group). Parasagittal slices were milled to 2.0 mm thickness and subjected to palmar forces 50 to 200% of those applied by the sesamoid bone at angles corresponding to early, mid and late stance support phases of the gait cycle. From contact radiographs in the loaded and unloaded samples, strains were determined by recognizing displacements in the trabecular patterns using texture correlation analysis. Failure did not occur in any of the samples. Strains were generally proportional to the forces applied and greatest at midstance. Strain patterns varied between samples and with the different loading positions. With increased subchondral bone sclerosis there was greater shear strain in overlying trabeculae. Strain patterns were not consistently different within the sclerotic bone at the site of failure. Focally higher strains at the surface were sometimes related to the edge of the platen which was molded to mimic the sesamoid bone in vivo. These results indicate that sclerotic thickening of subchondral bone transmits stresses to overlying trabeculae. No consistent strain pattern was recognized where devitalization and mechanical failure occurs. Focally higher strains related to the edge of the opposing sesamoid bone may play a role.