Integrating modelling, motion capture and x-ray fluoroscopy to investigate patellofemoral function during dynamic activity

The non-invasive evaluation technique of patellofemoral joint stress: a systematic literature review

Introduction: Patellofemoral joint stress (PFJS) is an important parameter for understanding the mechanism of patellofemoral joint pain, preventing patellofemoral joint injury, and evaluating the therapeutic efficacy of PFP rehabilitation programs. The purpose of this systematic review was to identify and categorize the non-invasive technique to evaluate the PFJS. Methods: Literature searches were conducted from January 2000 to October 2022 in electronic databases, namely, PubMed, Web of Science, and EBSCO (Medline, SPORTDiscus). This review includes studies that evaluated the patellofemoral joint reaction force (PJRF) or PFJS, with participants including both healthy individuals and those with patellofemoral joint pain, as well as cadavers with no organic changes. The study design includes cross-sectional studies, case-control studies, and randomized controlled trials. The JBI quality appraisal criteria tool was used to assess the risk of bias in the included studies. Results: In total, 5016 articles were identified in the database research and the citation network, and 69 studies were included in the review. Discussion: Researchers are still working to improve the accuracy of evaluation for PFJS by using a personalized model and optimizing quadriceps muscle strength calculations. In theory, the evaluation method of combining advanced computational and biplane fluoroscopy techniques has high accuracy in evaluating PFJS. The method should be further developed to establish the "gold standard" for PFJS evaluation. In practical applications, selecting appropriate methods and approaches based on theoretical considerations and ecological validity is essential.

Deciphering the "Art" in Modeling and Simulation of the Knee Joint: Variations in Model Development

The use of computational modeling to investigate knee joint biomechanics has increased exponentially over the last few decades. Developing computational models is a creative process where decisions have to be made, subject to the modelers' knowledge and previous experiences, resulting in the "art" of modeling. The long-term goal of the KneeHub project is to understand the influence of subjective decisions on the final outcomes and the reproducibility of computational knee joint models. In this paper, we report on the model development phase of this project, investigating model development decisions and deviations from initial modeling plans. Five teams developed computational knee joint models from the same dataset, and we compared each teams' initial uncalibrated models and their model development workflows. Variations in the software tools and modeling approaches were found, resulting in differences such as the representation of the anatomical knee joint structures in the model. The teams consistently defined the boundary conditions and used the same anatomical coordinate system convention. However, deviations in the anatomical landmarks used to define the coordinate systems were present, resulting in a large spread in the kinematic outputs of the uncalibrated models. The reported differences and similarities in model development and simulation presented here illustrate the importance of the "art" of modeling and how subjective decision-making can lead to variation in model outputs. All teams deviated from their initial modeling plans, indicating that model development is a flexible process and difficult to plan in advance, even for experienced teams.

Validated Computational Framework for Evaluation of In Vivo Knee Mechanics

Dynamic, in vivo evaluations of knee mechanics are important for understanding knee injury and repair, and developing successful treatments. Computational models have been used with in vivo experiments to quantify joint mechanics, but they are typically not predictive. The current study presents a novel integrated approach with high-speed stereo radiography, musculoskeletal modeling, and finite element (FE) modeling for evaluation of subject-specific, in vivo knee mechanics in a healthy subject performing a seated knee extension and weight-bearing lunge. Whole-body motion capture, ground reaction forces, and radiography-based kinematics were used to drive musculoskeletal and predictive FE models for load-controlled simulation of in vivo knee mechanics. A predictive simulation of knee mechanics was developed in four stages: (1) in vivo measurements of one subject performing a lunge and a seated knee extension, (2) rigid-body musculoskeletal modeling to determine muscle forces, (3) FE simulation of knee extension for knee-ligament calibration, and (4) predictive FE simulation of a lunge. FE models predicted knee contact and ligament mechanics and evaluated the impact of cruciate ligament properties on joint kinematics and loading. Calibrated model kinematics demonstrated good agreement to the experimental motion with root-mean-square differences of tibiofemoral flexion-extension <3 deg, internal-external <4 deg, and anterior-posterior <2 mm. Ligament reference strain and attachment locations were the most critical properties in the calibration process. The current work advances previous in vivo knee modeling through simulation of dynamic activities, modeling of subject-specific knee behavior, and development of a load-controlled knee model.

Prediction of patellofemoral joint kinematics and contact through co-simulation of rigid body dynamics and nonlinear finite element analysis

Joint-level rigid body dynamics simulations, when coupled with tissue-level finite element analyses, can simultaneously provide movement and tissue deformation metrics to understand mechanical interactions within the joint on a multi-scale level. In this study, a co-simulation workflow of a joint-level rigid body model that predicts the relative motion as a function of the non-linear cartilage response predicted by a non-linear implicit finite element solver is presented. Predictions are compared to in-vitro measurements (The Open Knee(s) project) in terms of the mean error and level-of-agreement: pressure_error = 0.46 MPa (level-of-agreement, -0.23 - 1.1 MPa); area_error = -89 mm² (level-of-agreement, -280 - 98 mm²) and contact force_error = 93 N (level-of-agreement, 7.8 - 180 N). The automated co-simulation control algorithm enables multiscale coupling between joint and tissue-level models with real-time two-way communication as opposed to the traditional feed-forward approach of multi-scale models.

Analytical and computational sliding wear prediction in a novel knee implant: a case study

Osteoarthritis (OA) is a commonly occurring cartilage degenerative disease. The end stage treatment is Total Knee Arthroplasty (TKA), which can be costly in terms of initial surgery, but also in terms of revision knee arthroplasty, which is quite often required. A novel conceptual knee implant has been proposed to function as a reducer of stress across the joint surface, to extend the period of time before TKA becomes necessary. The objective of this paper is to develop a computational model which can be used to assess the wear arising at the implant articulating surfaces. Experimental wear coefficients were determined from physical testing, the results of which were verified using a semi-analytical model. Experimental results were incorporated into an anatomically correct computational model of the knee and implant. The wear-rate predicted for the implant was 27.74 mm³ per million cycles (MC) and the wear depth predicted was 1.085 mm/MC. Whereas the wear-rate is comparable to that seen in conventional knee implants, the wear depth is significantly higher than for conventional knee prostheses, and indicates that, in order to be viable, wear-rates should be reduced in some way, perhaps by using low-wear polymers.

Dynamic Evaluation of Patellofemoral Instability: A Clinical Reality or Just a Research Field? A Literature review

Patellofemoral instability (PFI) is one of the most disabling conditions in the knee, often affecting young individuals. Despite its not uncommon presentation, the underlying biomechanical features leading to this entity are not entirely understood. The suitability of classic physical examination manoeuvres and imaging tests is a matter of discussion among treating surgeons, and so are the findings provided by these means. A potential cause for this lack of consensus is the fact that, classically, the diagnostic approach for PFI has relied on statically obtained data. Many authors advocate for the study of this entity in a dynamic scenario, closer to the actual situation in which the instability episodes occur. In this literature review, we have compiled the available data from the last decades regarding dynamic evaluation methods for PFI and related conditions. Several categories are presented, grouping the related techniques and devices: physical examination, imaging modalities (ultrasound (US), magnetic resonance imaging (MRI), computed tomography (CT) and combined methods), arthroscopic evaluation, and others. In conclusion, although a vast number of quality studies are presented, in which comprehensive data about the biomechanics of the patellofemoral joint (PFJ) are described, this evidence has not yet reached clinical practice universally. Most of the data still stays in the research field and is seldom employed to assist a better understanding of the PFI cases and their ideal treatment targets.

Analysis of the moment arms and kinematics of ostrich (Struthio camelus) double patellar sesamoids

The patella ("kneecap") is a biomechanically important feature of the tendinous insertion of the knee extensor muscles, able to alter the moment arm lengths between its input and output tendons, and so modify the mechanical advantage of the knee extensor muscle. However, patellar gearing function is little-explored outside of humans, and the patella is often simplified or ignored in biomechanical models. Here, we investigate patellar gearing and kinematics in the ostrich-frequently used as an animal analogue to human bipedal locomotion and unusual in its possession of two patellae at the knee joint. We use x-ray reconstruction of moving morphology (XROMM) techniques to capture the kinematics of the patellae in an adult ostrich cadaver, passively manipulated in flexion-extension. Moment arm ratios between the input and output tendons of each patella are calculated from kinematically determined centers of patellofemoral joint rotation. Both patellae are found to decrease the mechanical advantage of the extensor muscle-tendon complex, decreasing the tendon output force for a given muscle input force, but potentially increasing the relative speed of knee extension. Mechanically and kinematically, the proximal patella behaves similarly to the single patella of most other species, whereas the distal patella has properties of both a fixed retroarticular process and a moving sesamoid. It is still not clear why ostriches possess two patellae, but we suggest that the configuration in ostriches benefits their rapid locomotion and provides tendon protection.

Combined measurement and modeling of specimen-specific knee mechanics for healthy and ACL-deficient conditions

Quantifying the mechanical environment at the knee is crucial for developing successful rehabilitation and surgical protocols. Computational models have been developed to complement in vitro studies, but are typically created to represent healthy conditions, and may not be useful in modeling pathology and repair. Thus, the objective of this study was to create finite element (FE) models of the natural knee, including specimen-specific tibiofemoral (TF) and patellofemoral (PF) soft tissue structures, and to evaluate joint mechanics in intact and ACL-deficient conditions. Simulated gait in a whole joint knee simulator was performed on two cadaveric specimens in an intact state and subsequently repeated following ACL resection. Simulated gait was performed using motor-actuated quadriceps, and loads at the hip and ankle. Specimen-specific FE models of these experiments were developed in both intact and ACL-deficient states. Model simulations compared kinematics and loading of the experimental TF and PF joints, with average RMS differences [max] of 3.0° [8.2°] and 2.1° [8.4°] in rotations, and 1.7 [3.0] and 2.5 [5.1] mm in translations, for intact and ACL-deficient states, respectively. The timing of peak quadriceps force during stance and swing phase of gait was accurately replicated within 2° of knee flexion and with an average error of 16.7% across specimens and pathology. Ligament recruitment patterns were unique in each specimen; recruitment variability was likely influenced by variations in ligament attachment locations. ACL resections demonstrated contrasting joint mechanics in the two specimens with altered knee motion shown in one specimen (up to 5mm anterior tibial translation) while increased TF joint loading was shown in the other (up to 400N).

Mobile Biplane X-Ray Imaging System for Measuring 3D Dynamic Joint Motion During Overground Gait

Most X-ray fluoroscopy systems are stationary and impose restrictions on the measurement of dynamic joint motion; for example, knee-joint kinematics during gait is usually measured with the subject ambulating on a treadmill. We developed a computer-controlled, mobile, biplane, X-ray fluoroscopy system to track human body movement for high-speed imaging of 3D joint motion during overground gait. A robotic gantry mechanism translates the two X-ray units alongside the subject, tracking and imaging the joint of interest as the subject moves. The main aim of the present study was to determine the accuracy with which the mobile imaging system measures 3D knee-joint kinematics during walking. In vitro experiments were performed to measure the relative positions of the tibia and femur in an intact human cadaver knee and of the tibial and femoral components of a total knee arthroplasty (TKA) implant during simulated overground gait. Accuracy was determined by calculating mean, standard deviation and root-mean-squared errors from differences between kinematic measurements obtained using volumetric models of the bones and TKA components and reference measurements obtained from metal beads embedded in the bones. Measurement accuracy was enhanced by the ability to track and image the joint concurrently. Maximum root-mean-squared errors were 0.33 mm and 0.65° for translations and rotations of the TKA knee and 0.78 mm and 0.77° for translations and rotations of the intact knee, which are comparable to results reported for treadmill walking using stationary biplane systems. System capability for in vivo joint motion measurement was also demonstrated for overground gait.

Dynamic CT technique for assessment of wrist joint instabilities

PURPOSE

To develop a 4D [three-dimensional (3D) + time] CT technique to capture high spatial and temporal resolution images of wrist joint motion so that dynamic joint instabilities can be detected before the development of static joint instability and onset of osteoarthritis (OA).

METHODS

A cadaveric wrist was mounted onto a custom motion simulator and scanned with a dual source CT scanner during radial-ulnar deviation. A dynamic 4D CT technique was utilized to reconstruct images at 20 equidistant time points from one motion cycle. 3D images of carpal bones were generated using volume rendering techniques (VRT) at each of the 20 time points and then 4D movies were generated to depict the dynamic joint motion. The same cadaveric wrist was also scanned after cutting all portions of the scapholunate interosseus ligament to simulate scapholunate joint instability. Image quality were assessed on an ordinal scale (1-4, 4 being excellent) by three experienced orthopedic surgeons (specialized in hand surgery) by scoring 2D axial images. Dynamic instability was evaluated by the same surgeons by comparing the two 4D movies of joint motion. Finally, dose reduction was investigated using the cadaveric wrist by scanning at different dose levels to determine the lowest radiation dose that did not substantially alter diagnostic image quality.

RESULTS

The mean image quality scores for dynamic and static CT images were 3.7 and 4.0, respectively. The carpal bones, distal radius and ulna, and joint spaces were clearly delineated in the 3D VRT images, without motion blurring or banding artifacts, at all time points during the motion cycle. Appropriate viewing angles could be interactively selected to view any articulating structure using different 3D processing techniques. The motion of each carpal bone and the relative motion among the carpal bones were easily observed in the 4D movies. Joint instability was correctly and easily detected in the scan performed after the ligament was cut by observing the relative motion between the scaphoid and lunate bones. Diagnostic capability was not sacrificed with a volume CT dose index (CTDI(vol)) as low as 18 mGy for the whole scan, with estimated skin dose of approximately 33 mGy, which is much lower than the threshold for transient skin erythema (2000 mGy).

CONCLUSIONS

The proposed dynamic 4D CT imaging technique generated high spatial and high temporal resolution images without requiring periodic joint motion. Preliminary results from this cadaveric study demonstrate the feasibility of detecting joint instability using this technique.

Three-dimensional geometry of the human biceps femoris long head measured in vivo using magnetic resonance imaging

BACKGROUND

The human biceps femoris long head is susceptible to injury, especially when sprinting. The potential mechanical action of this muscle at a critical stage in the stride cycle was evaluated by calculating three-dimensional lines-of-action and moment arms about the hip and knee joints in vivo.

METHODS

Axial magnetic resonance images of the right lower-limb (pelvis to proximal tibia) were recorded from four participants under two conditions: a reference pose, with the lower-limb in the anatomical position and the hamstrings relaxed; and a terminal swing pose, with the hip and knee joints flexed to mimic the lower-limb orientation during the terminal swing phase of sprinting and the hamstrings isometrically activated. Images were used to segment biceps femoris long head and the relevant bones. The musculotendon path and joint coordinate systems were defined from which lines-of-action and moment arms were computed.

FINDINGS

Biceps femoris long head displayed hip extensor and adductor moment arms as well as knee flexor, abductor and external-rotator moment arms. Sagittal-plane moment arms were largest, whereas transverse-plane moment arms were smallest. Moment arms remained consistent in polarity across all participants and testing conditions, except in the transverse-plane about the hip. For the terminal swing pose compared to the reference pose, sagittal-plane moment arms for biceps femoris long head increased by 19.9% to 48.9% about the hip and 42.3% to 93.9% about the knee.

INTERPRETATION

Biceps femoris long head has the potential to cause hip extension and adduction as well as knee flexion during the terminal swing phase of sprinting.

Radiological method for measuring patellofemoral tracking and tibiofemoral kinematics before and after total knee replacement

OBJECTIVES

Numerous complications following total knee replacement (TKR) relate to the patellofemoral (PF) joint, including pain and patellar maltracking, yet the options for in vivo imaging of the PF joint are limited, especially after TKR. We propose a novel sequential biplane radiological method that permits accurate tracking of the PF and tibiofemoral (TF) joints throughout the range of movement under weightbearing, and test it in knees pre- and post-arthroplasty.

METHODS

A total of three knees with end-stage osteoarthritis and three knees that had undergone TKR at more than one year's follow-up were investigated. In each knee, sequential biplane radiological images were acquired from the sagittal direction (i.e. horizontal X-ray source and 10° below horizontal) for a sequence of eight flexion angles. Three-dimensional implant or bone models were matched to the biplane images to compute the six degrees of freedom of PF tracking and TF kinematics, and other clinical measures.

RESULTS

The mean and standard deviation for the six degrees of freedom of PF tracking and TF kinematics were computed. TF and PF kinematics were highly accurate (< 0.9 mm, < 0.6°) and repeatable.

CONCLUSIONS

The developed method permitted measuring of in vivo PF tracking and TF kinematics before and after TKR throughout the range of movement. This method could be a useful tool for investigating differences between cohorts of patients (e.g., with and without pain) impacting clinical decision-making regarding surgical technique, revision surgery or implant design.

Mechanics of the foot Part 1: a continuum framework for evaluating soft tissue stiffening in the pathologic foot

Soft tissue stiffening is a common mechanical observation reported in foot pathologies including diabetes mellitus and gout. These material changes influence the spatial distribution of stress and affect blood flow, which is essential to nutrient entry and waste removal. An anatomically-based subject-specific foot model was developed to explore the influence of tissue stiffening on plantar pressure and internal von Mises stress at heel-strike, midstance and toe-off. This work draws on the model database developed for the Physiome project consisting of muscles, bones, soft tissue and other structures such as sensory nerves. The anisotropic structure of soft tissue was embedded in a single continuum as an efficient model for finite soft tissue deformation, and customisation methods were used to capture the unique foot profile. The model was informed by kinetics from an instrumented treadmill and kinematics from motion capture, synchronised together. Foot sole pressure predictions were evaluated against a commercial pressure platform. Key outcomes showed that internal stress can be up to 1.6 times the surface pressure with implications for internal soft tissue damage not observed at the surface. The main nerve branch stimulated during gait was the lateral plantar nerve. This subject-specific modelling framework can play an integral part in therapeutic treatments by informing assistive strategies such as mechanical noise stimulation and orthotics.

Integrated assessment techniques for linking kinematics, kinetics and muscle activation to early migration: a pilot study

The goal of this pilot study was to develop and test an integrated method to assess kinematics, kinetics and muscle activation of total knee prostheses during dynamic activities, by integrating fluoroscopic measurements with force plate, electromyography and external motion registration measurements. Subsequently, this multi-instrumental analysis was then used to assess the relationship between kinematics, kinetics and muscle activation and early migration of the tibial component of total knee prostheses. This pilot study showed that it is feasible to integrate fluoroscopic, kinematic and kinetic measurements and relate findings to early migration data. Results showed that there might be an association between deviant kinematics and early migration in patients with a highly congruent mobile-bearing total knee prosthesis. Patients that showed high levels of coactivation, diverging axial rotations of the insert and a deviant pivot point showed increased migration and might be at higher risk for tibial component loosening. In the future, to confirm our findings, the same integrated measurements have to be performed in larger patient groups and different prosthesis designs.

PATELLA TRACKING WITH PERIPATELLAR SOFT TISSUE STABILIZERS AS A FUNCTION OF DYNAMIC SUBJECT-SPECIFIC KNEE FLEXION

Musculoskeletal modeling has found wide application in joint biomechanics investigations. This technique has been improved by incorporating subject-specific skeletal elements and passive patellofemoral stabilizers in a dynamic analysis. After trochlear engagement, the volunteers' patellae displaced laterally, whereas tilt was subject specific. Comparison of the tilt and mediolateral position values to in vivo MRI values at 30° knee flexion showed a mean accuracy of 84.4% and 96.9%, respectively. Medial patellofemoral ligament tension decreased with knee flexion, while the patellar tendon–quadriceps tendon ratio ranged from 0.4 to 1.2. The patellofemoral contact load–quadriceps tendon load ratio ranged from 0.7 to 1.3, whereas the mediolateral load component–resultant load ratio ranged from 0 to 0.4. Three validated subject-specific musculoskeletal models facilitated the analysis of patellofemoral biomechanics: Subject-specific patella tracking and passive stabilizer response was analyzed as a function of dynamic knee flexion.

The quality of bone surfaces may govern the use of model based fluoroscopy in the determination of joint laxity

The assessment of knee joint laxity is clinically important but its quantification remains elusive. Calibrated, low dosage fluoroscopy, combined with registered surfaces and controlled external loading may offer possible solutions for quantifying relative tibio-femoral motion without soft tissue artefact, even in native joints. The aim of this study was to determine the accuracy of registration using CT and MRI derived 3D bone models, as well as metallic implants, to 2D single-plane fluoroscopic datasets, to assess their suitability for examining knee joint laxity. Four cadaveric knees and one knee implant were positioned using a micromanipulator. After fluoroscopy, the accuracy of registering each surface to the 2D fluoroscopic images was determined by comparison against known translations from the micromanipulator measurements. Dynamic measurements were also performed to assess the relative tibio-femoral error. For CT and MRI derived 3D femur and tibia models during static testing, the in-plane error was 0.4 mm and 0.9 mm, and out-of-plane error 2.6 mm and 9.3 mm respectively. For metallic implants, the in-plane error was 0.2 mm and out-of-plane error 1.5 mm. The relative tibio-femoral error during dynamic measurements was 0.9 mm, 1.2 mm and 0.7 mm in-plane, and 3.9 mm, 10.4 mm and 2.5 mm out-of-plane for CT and MRI based models and metallic implants respectively. The rotational errors ranged from 0.5° to 1.9° for CT, 0.5-4.3° for MRI and 0.1-0.8° for metallic implants. The results of this study indicate that single-plane fluoroscopic analysis can provide accurate information in the investigation of knee joint laxity, but should be limited to static or quasi-static evaluations when assessing native bones, where possible. With this knowledge of registration accuracy, targeted approaches for the determination of tibio-femoral laxity could now determine objective in vivo measures for the identification of ligament reconstruction candidates as well as improve our understanding of the consequences of knee joint instability in TKA.

The development of lower limb musculoskeletal models with clinical relevance is dependent upon the fidelity of the mathematical description of the lower limb. Part 2: patient-specific geometry

Musculoskeletal models have the potential to evolve into sensitive clinical tools that provide relevant therapeutic guidance. A key impediment to this is the lack of understanding as to the function of such models. In order to improve this it is useful to recognise that musculoskeletal modelling is the mathematical description of musculoskeletal movement – a process that involves the construction and solution of equations of motion. These equations are derived from standard mechanical considerations and the mathematical representation of anatomy. The fidelity of musculoskeletal models is highly dependent on the assumption that such representations also describe the function of the musculoskeletal geometry. In addition, it is important to understand the sensitivity of such representations to patient-specific variations in anatomy. The exploration of these twin considerations will be fundamental to the creation of musculoskeletal modelling tools with clinical relevance and a systematic enquiry of these key parameters is recommended.

Future trends in the use of X-ray fluoroscopy for the measurement and modelling of joint motion

Knowledge of three-dimensional skeletal kinematics during functional activities such as walking, is required for accurate modelling of joint motion and loading, and is important in identifying the effects of injury and disease. For example, accurate measurement of joint kinematics is essential in understanding the pathogenesis of osteoarthritis and its symptoms and for developing strategies to alleviate joint pain. Bi-plane X-ray fluoroscopy has the capacity to accurately and non-invasively measure human joint motion in vivo. Joint kinematics obtained using bi-plane X-ray fluoroscopy will aid in the development of more complex musculoskeletal models, which may be used to assess joint function and disease and plan surgical interventions and post-operative rehabilitation strategies. At present, however, commercial C-arm systems constrain the motion of the subject within the imaging field of view, thus precluding recording of motions such as overground gait. These fluoroscopy systems also operate at low frame rates and therefore cannot accurately capture high-speed joint motion during tasks such as running and throwing. In the future, bi-plane fluoroscopy systems may include computer-controlled tracking for the measurement of joint kinematics over entire cycles of overground gait without constraining motion of the subject. High-speed cameras will facilitate measurement of high-impulse joint motions, and computationally efficient pose-estimation software may provide a fast and fully automated process for quantification of natural joint motion.

Model predictions of increased knee joint loading in regions of thinner articular cartilage after patellar tendon adhesion

Patellar tendon adhesion is a complication from anterior cruciate ligament (ACL) reconstruction that may affect patellofemoral and tibiofemoral biomechanics. A computational model was used to investigate the changes in knee joint mechanics due to patellar tendon adhesion under normal physiological loading during gait. The calculations showed that patellar tendon adhesion up to the level of the anterior tibial plateau led to patellar infera, increased patellar flexion, and increased anterior tibial translation. These kinematic changes were associated with increased patellar contact force, a distal shift in peak patellar contact pressure, a posterior shift in peak tibial contact pressure, and increased peak tangential contact sliding distance over one gait cycle (i.e., contact slip). Postadhesion, patellar and tibial contact locations corresponded to regions of thinner cartilage. The predicted distal shift in patellar contact was in contrast to other patellar infera studies. Average patellar and tibial cartilage pressure did not change significantly following patellar tendon adhesion; however, peak medial tibial pressure increased. These results suggest that changes in peak tibial cartilage pressure, contact slip, and the migration of contact to regions of thinner cartilage are associated with patellar tendon adhesion and may be responsible for initiating patellofemoral pain and knee joint structural damage observed following ACL reconstruction.

Assessing the accuracy and precision of musculoskeletal motion tracking using cine-PC MRI on a 3.0T platform

The rising cost of musculoskeletal pathology, disease, and injury creates a pressing need for accurate and reliable methods to quantify 3D musculoskeletal motion, fostering a renewed interest in this area over the past few years. To date, cine-phase contrast (PC) MRI remains the only technique capable of non-invasively tracking in vivo 3D musculoskeletal motion during volitional activity, but current scan times are long on the 1.5T MR platform (∼ 2.5 min or 75 movement cycles). With the clinical availability of higher field strength magnets (3.0T) that have increased signal-to-noise ratios, it is likely that scan times can be reduced while improving accuracy. Therefore, the purpose of this study is to validate cine-PC MRI on a 3.0T platform, in terms of accuracy, precision, and subject-repeatability, and to determine if scan time could be minimized. On the 3.0T platform it is possible to limit scan time to 2 min, with sub-millimeter accuracy (<0.33 mm/0.97°), excellent technique precision (<0.18°), and strong subject-repeatability (<0.73 mm/1.10°). This represents reduction in imaging time by 25% (42 s), a 50% improvement in accuracy, and a 72% improvement in technique precision over the original 1.5T platform. Scan time can be reduced to 1 min (30 movement cycles), but the improvements in accuracy are not as large.

Automatic determination of anatomical coordinate systems for three-dimensional bone models of the isolated human knee

The combination of three-dimensional (3-D) models with dual fluoroscopy is increasingly popular for evaluating joint function in vivo. Applying these modalities to study knee motion with high accuracy requires reliable anatomical coordinate systems (ACSs) for the femur and tibia. Therefore, a robust method for creating ACSs from 3-D models of the femur and tibia is required. We present and evaluate an automated method for constructing ACSs for the distal femur and proximal tibia based solely on 3-D bone models. The algorithm requires no observer interactions and uses model cross-sectional area, center of mass, principal axes of inertia, and cylindrical surface fitting to construct the ACSs. The algorithm was applied to the femur and tibia of 10 (unpaired) human cadaveric knees. Due to the automated nature of the algorithm, the within specimen variability is zero for a given bone model. The algorithm's repeatability was evaluated by calculating variability in ACS location and orientation across specimens. Differences in ACS location and orientation between specimens were low (<1.5mm and <2.5 degrees). Variability arose primarily from natural anatomical and morphological differences between specimens. The presented algorithm provides an alternative method for automatically determining subject-specific ACSs from the distal femur and proximal tibia.

Predicting the effects of knee focal articular surface injury with a patient-specific finite element model

Successful focal articular surface injury (FAI) repair depends on appropriate matching of the geometrical/material properties of the repaired site, and on the overall dynamic response of the knee to in-vivo loading. There is evidence linking the pathogenesis of lesion progression (e.g. osteoarthritis) to weightbearing site and defect size. The paper investigates further this link by studying the effects of osteochondral defect size on the load distribution at the human knee. Experimental data from cadaver knees (n=8) loaded at 30 degrees of flexion was used as input to a validated finite element (FE) model. Contact pressure was assessed for the intact knees and over a range of circular osteochondral defects (5 mm to 20 mm) at 30 degrees of flexion with 700 N axial load. Patient specific FE models and the specific boundary conditions of the experimental set-up were used to analyze the osteochondral defects. Stress concentration around the rims of defects 8 mm and smaller was not significant and pressure distribution was dominated by the menisci. Experimental data was confirmed by the model. For defects 10 mm and greater, distribution of peak pressures followed the rim of the defect with a mean distance from the rim of 2.64 mm on the medial condyle and 2.90 mm on the lateral condyle (model predictions were 2.63 and 2.87 mm respectively). Statistical significance was reported when comparing defects that differed by 4 mm or greater (except for the 5 mm case). Peak rim pressure did not significantly increase as defects were enlarged from 10 mm to 20 mm. Peak values were always significantly higher over the medial femoral condyle. Although the decision to treat osteochondral lesions is multifactorial, the results of this finite element analysis indicate that a size threshold of 10 mm, may be a useful early adjunct to guide clinical decision-making. This modified FE method can be employed for in-vivo studies.

Evaluation of predicted knee-joint muscle forces during gait using an instrumented knee implant

Musculoskeletal modeling and optimization theory are often used to determine muscle forces in vivo. However, convincing quantitative evaluation of these predictions has been limited to date. The present study evaluated model predictions of knee muscle forces during walking using in vivo measurements of joint contact loading acquired from an instrumented implant. Joint motion, ground reaction force, and tibial contact force data were recorded simultaneously from a single subject walking at slow, normal, and fast speeds. The body was modeled as an 8-segment, 21-degree-of-freedom articulated linkage, actuated by 58 muscles. Joint moments obtained from inverse dynamics were decomposed into leg-muscle forces by solving an optimization problem that minimized the sum of the squares of the muscle activations. The predicted knee muscle forces were input into a 3D knee implant contact model to calculate tibial contact forces. Calculated and measured tibial contact forces were in good agreement for all three walking speeds. The average RMS errors for the medial, lateral, and total contact forces over the entire gait cycle and across all trials were 140 +/- 40 N, 115 +/- 32 N, and 183 +/- 45 N, respectively. Muscle coordination predicted by the model was also consistent with EMG measurements reported for normal walking. The combined experimental and modeling approach used in this study provides a quantitative framework for evaluating model predictions of muscle forces in human movement.