Validation of radiocarpal joint contact models based on images from a clinical MRI scanner

This study was undertaken to assess magnetic resonance imaging (MRI)-based radiocarpal surface contact models of functional loading in a clinical MRI scanner for future in vivo studies, by comparison with experimental measures from three cadaver forearm specimens. Experimental data were acquired using a Tekscan sensor during simulated light grasp. Magnetic resonance (MR) images were used to obtain model geometry and kinematics (image registration). Peak contact pressures (PPs) and average contact pressures (APs), contact forces and contact areas were determined in the radiolunate and radioscaphoid joints. Contact area was also measured directly from MR images acquired with load and compared with model data. Based on the validation criteria (within 25% of experimental data), out of the six articulations (three specimens with two articulations each), two met the criterion for AP (0%, 14%); one for peak pressure (20%); one for contact force (5%); four for contact area with respect to experiment (8%, 13%, 19% and 23%), and three contact areas met the criterion with respect to direct measurements (14%, 21% and 21%). Absolute differences between model and experimental PPs were reasonably low (within 2.5 MPa). Overall, the results indicate that MRI-based models generated from 3T clinical MR scanner appear sufficient to obtain clinically relevant data.     

MRI-based modeling for evaluation of in vivo contact mechanics in the human wrist during active light grasp

Investigations of in vivo joint mechanics are important for understanding the joint function under functional loading and the mechanisms of pathology. In this study we used magnetic resonance imaging (MRI) based joint contact modeling to evaluate in vivo joint contact mechanics in the human wrist. MRI scans were performed on the wrists of four subjects while they maintained light grasp of a cylinder, and with the same wrist relaxed. 3D models of the radius, scaphoid and lunate, including cartilage surface data, were constructed from the relaxed image data. These models were transformed into the loaded configuration, as determined from the grasp image data, and contact mechanics were evaluated. The resulting contact pressures, areas and forces were then analyzed for each articulation and for each subject. Contact areas were measured directly from grasp MRI images for comparison to the model predictions. The first-ever estimates for in vivo radioscaphoid and radiolunate contact pressure agreed reasonably well with previous cadaveric studies. This investigation also produced novel in vivo scapholunate contact results that were similar to radiolunate data. The specimen-specific contact area comparison generally showed substantial variability between the models and the direct measurements from MRI. On average, the models were within about 10% of the direct MRI measurements for radioscaphoid and scapholunate contact areas, but radiolunate contact areas from the model were only within 55% of the direct measurements. Overall, the results of the study suggest that MRI-based modeling has substantial potential for evaluation of in vivo joint contact mechanics, especially as technology and methodology improve.     

Physical validation of a patient-specific contact finite element model of the ankle

A validation study was conducted to determine the extent to which computational ankle contact finite element (FE) results agreed with experimentally measured tibio-talar contact stress. Two cadaver ankles were loaded in separate test sessions, during which ankle contact stresses were measured with a high-resolution (Tekscan) pressure sensor. Corresponding contact FE analyses were subsequently performed for comparison. The agreement was good between FE-computed and experimentally measured mean (3.2% discrepancy for one ankle, 19.3% for the other) and maximum (1.5% and 6.2%) contact stress, as well as for contact area (1.7% and 14.9%). There was also excellent agreement between histograms of fractional areas of cartilage experiencing specific ranges of contact stress. Finally, point-by-point comparisons between the computed and measured contact stress distributions over the articular surface showed substantial agreement, with correlation coefficients of 90% for one ankle and 86% for the other. In the past, general qualitative, but little direct quantitative agreement has been demonstrated with articular joint contact FE models. The methods used for this validation enable formal comparison of computational and experimental results, and open the way for objective statistical measures of regional correlation between FE-computed contact stress distributions from comparison articular joint surfaces (e.g., those from an intact versus those with residual intra-articular fracture incongruity).     

A new technique to measure the dynamic contact pressures on the Tibial Plateau

Studies of the load transfer role of the meniscus have been limited to static experimental and analytical approaches. The objective of this study was to develop an experimental technique to allow the contact pressures on the tibial plateau of cadaveric knees to be measured under dynamic physiological loads. Accordingly, we adapted a load-controlled knee joint simulator to accept a cadaveric sheep knee, programmed the simulator with sheep gait kinematics data, and utilized a pressure sensor array to measure the contact pressure distribution on the lateral tibial plateau during gait. The technique was applied to six sheep knees that were tested intact and after meniscectomy. Meniscectomy resulted in a 267% increase in average contact pressure, a 117% increase in peak contact pressure, and an 80% decrease in contact area, all measured at the point of maximum peak contact stress in the gait cycle. It is envisaged that the experimental model herein developed will allow for the screening of candidate materials prior to more expensive and time-consuming animal models.     

Explicit finite element modeling of total knee replacement mechanics

Joint kinematics and contact mechanics dictate the success of current total knee replacement (TKR) devices. Efficient computer models present an effective way of evaluating these characteristics. Predicted contact stress and area due to articulations at the tibio-femoral and patello-femoral interfaces indicate potential clinical performance. Previous finite element (FE) knee models have generally been used to predict contact stresses and/or areas during static or quasi-static loading conditions. Explicit dynamic FE analyses have recently been used to efficiently predict TKR kinematics and contact mechanics during dynamic loading conditions. The objective of this study was to develop and experimentally validate an explicit FE TKR model that incorporates tibio-femoral and patello-femoral articulations. For computational efficiency, we developed rigid body analyses that can reasonably reproduce the kinematics, contact pressure distribution, and contact area of a fully deformable system. Results from the deformable model showed that the patello-femoral and tibio-femoral kinematics were in good agreement with experimental knee simulator measurements. Kinematic results from the rigid body analyses were nearly identical to those from the fully deformable model, and the contact pressure and contact area correlation was acceptable given the great reduction in analysis time. Component mesh density studied had little effect on the predicted kinematics, particularly for the patellar component, and small effects on the predicted contact pressure and area. These analyses have shown that, at low computational cost, a force-control dynamic simulation of a gait cycle can yield useful and predictable results.     

Effect of implementing magnetic resonance imaging for patient-specific OpenSim models on lower-body kinematics and knee ligament lengths

BackgroundOpenSim models are typically based on cadaver findings that are generalized to represent a wide range of populations, which curbs their validity. Patient-specific modelling through incorporating magnetic resonance imaging (MRI) improves the model’s biofidelity with respect to joint alignment and articulations, muscle wrapping, and ligament insertions. The purpose of this study was to determine if the inclusion of an MRI-based knee model would elicit differences in lower limb kinematics and resulting knee ligament lengths during a side cut task.MethodsEleven participants were analyzed with the popular Rajagopal OpenSim model, two variations of the same model to include three and six degrees of freedom knee (DOF), and a fourth version featuring a four DOF MRI-based knee model. These four models were used in an inverse kinematics analysis of a side cut task and the resulting lower limb kinematics and knee ligament lengths were analyzed.ResultsThe MRI-based model was more responsive to the movement task than the original Rajagopal model while less susceptible to soft tissue artifact than the unconstrained six DOF model. Ligament isometry was greatest in the original Rajagopal model and smallest in the six DOF model.ConclusionsWhen using musculoskeletal modelling software, one must acutely consider the model choice as the resulting kinematics and ligament lengths are dependent on this decision. The MRI-based knee model is responsive to the kinematics and ligament lengths of highly dynamic tasks and may prove to be the most valid option for continuing with late-stage modelling operations such as static optimization.     

Use of distraction loading to estimate subject-specific knee ligament slack lengths

Knee ligaments guide and restrain joint motion, and their properties influence joint mechanics. Inverse modeling schemes have been used to estimate specimen-specific ligament properties, where external joint forces are assumed to balance with internal ligament and contact forces. This study simplifies this assumption by adjusting experimental loads to remove internal contact forces. The purpose of this study was to use novel experimental loading in an inverse modeling scheme to estimate ligament slack lengths, perform validation using additional loading scenarios, and evaluate sensitivity to the applied loading. Joint kinematics and kinetics were experimentally measured for a set of load cases. An optimization scheme used a specimen-specific forward kinematics model to estimate ligament slack lengths by minimizing the residual between model and experimentally measured kinetics. The calibrated model was used for a form of validation by evaluating non-optimized load cases. Additionally, uncertainty analysis related kinetic errors to previously reported kinematic errors. The six DOF tibial reactions realized RMS errors less than 23 N and 0.75 Nm for optimized load cases, and 33 N and 2.25 Nm for the non-optimized load cases. The uncertainty analysis, which was performed using the optimized load cases, showed average kinetic RMS errors less than 26 N and 0.45 Nm. The model’s recruitment patterns were similar to those found in clinical and cadaveric studies. This study demonstrated that experimental distraction loading can be used in an inverse modeling scheme to estimate ligament slack lengths with a forward kinematics model.     

A knee-specific finite element analysis of the human anterior cruciate ligament impingement against the femoral intercondylar notch

This work presents a finite element analysis of anterior cruciate ligament (ACL) impingement against the intercondylar notch during tibial external rotation and abduction, as a mechanism of noncontact ACL injuries. Experimentally, ACL impingement was measured in a cadaveric knee in terms of impingement contact pressure and six degrees-of-freedom tibiofemoral kinematics. Three-dimensional geometries of the ACL, femur and tibia were incorporated into the finite element model of the individual knee specimen. A fiber-reinforced model was adopted, which accounts for the anisotropy, large deformation, nonlinearity and incompressibility of the ACL. With boundary conditions specified based on the experimental tibiofemoral kinematics, the finite element analysis showed that impingement between the ligament and the lateral wall of intercondylar notch could occur when qthe knee at 45° was externally rotated at 29.1° and abducted at 10.0°. Strong contact pressure and tensile stress occurred at the impinging and nonimpinging sides of the ligament, respectively. The impingement force and contact area estimated from the model matched their counterparts from the corresponding cadaver experiment. The modeling and experimental approach provides a useful tool to characterize potential ACL impingement on a knee-specific basis, which may help elucidate the ACL injury mechanism and develop more effective treatments.     

Specimen-specific predictions of contact stress under physiological loading in the human hip: validation and sensitivity studies

Hip osteoarthritis may be initiated and advanced by abnormal cartilage contact mechanics, and finite element (FE) modeling provides an approach with the potential to allow the study of this process. Previous FE models of the human hip have been limited by single specimen validation and the use of quasi-linear or linear elastic constitutive models of articular cartilage. The effects of the latter assumptions on model predictions are unknown, partially because data for the instantaneous behavior of healthy human hip cartilage are unavailable. The aims of this study were to develop and validate a series of specimen-specific FE models, to characterize the regional instantaneous response of healthy human hip cartilage in compression, and to assess the effects of material nonlinearity, inhomogeneity and specimen-specific material coefficients on FE predictions of cartilage contact stress and contact area. Five cadaveric specimens underwent experimental loading, cartilage material characterization and specimen-specific FE modeling. Cartilage in the FE models was represented by average neo-Hookean, average Veronda Westmann and specimen- and region-specific Veronda Westmann hyperelastic constitutive models. Experimental measurements and FE predictions compared well for all three cartilage representations, which was reflected in average RMS errors in contact stress of less than 25 %. The instantaneous material behavior of healthy human hip cartilage varied spatially, with stiffer acetabular cartilage than femoral cartilage and stiffer cartilage in lateral regions than in medial regions. The Veronda Westmann constitutive model with average material coefficients accurately predicted peak contact stress, average contact stress, contact area and contact patterns. The use of subject- and region-specific material coefficients did not increase the accuracy of FE model predictions. The neo-Hookean constitutive model underpredicted peak contact stress in areas of high stress. The results of this study support the use of average cartilage material coefficients in predictions of cartilage contact stress and contact area in the normal hip. The regional characterization of cartilage material behavior provides the necessary inputs for future computational studies, to investigate other mechanical parameters that may be correlated with OA and cartilage damage in the human hip. In the future, the results of this study can be applied to subject-specific models to better understand how abnormal hip contact stress and contact area contribute to OA.     

Changes in talocrural joint contact stress characteristics after simulated rotationplasty

In order to assess the changes in talocrural joint contact stress after rotationplasty, 10 lower-leg cadaver specimens were axially loaded with 600 N and investigated in two loading situations: (1) Normal loading with a plantigrade foot; (2) in an equinus position of a simulated rotationplasty. Joint contact stress in the talar facet of the talocrural joint was determined with Fuji Prescale film cut to size and analyzed with digital image analysis for joint contact area, mean and peak pressure, contact force, and location of the load application on the trochlea tali. The results demonstrate a significant transfer on the loading zone to the posterior part of the talus (p = 0.005), a significant reduction of the contact area (p = 0.005) and force (p = 0.005), and a significant increase of the mean (p = 0.022) and maximum pressures (p = 0.013). These results indicate that the rotationplasty causes pronounced changes in joint loading characteristics.     

Effect of fluid boundary conditions on joint contact mechanics and applications to the modeling of osteoarthritic joints

The long-term goal of our research is to understand the mechanism of osteoarthritis (OA) initiation and progress through experimental and theoretical approaches. In previous theoretical models, joint contact mechanics was implemented without consideration of the fluid boundary conditions and with constant permeability. The primary purpose of this study was to investigate the effect of fluid boundary conditions at the articular surfaces on the contact mechanics, in terms of load sharing and fluid flow properties using variable permeability. The tested conditions included totally sealed surfaces, open surfaces, and open surfaces with variable permeability. While the sealed surface model failed to predict relaxation times and load sharing properly, the class of open surface models (open surfaces with constant permeability, and surfaces with variable permeability) gave good agreement with experiments, in terms of relaxation time and load sharing between the solid and the fluid phase. In particular, the variable permeability model was judged to be the most realistic of the three models, from a biological and physical point of view. This model was then used to simulate joint contact in the early and late stages of OA. In the early stages of OA, the model predicted a decrease in peak contact pressure and an increase in contact area, while in the late stages of OA, peak pressures were increased and contact areas were decreased compared to normal. These findings agree well with experimental observations.     

Modeling of laser coagulation of tissue with MRI temperature monitoring

Light energy from a laser source that is delivered into body tissue via a fiber-optic probe with minimal invasiveness has been used to ablate solid tumors. This thermal coagulation process can be guided and monitored accurately by continuous magnetic resonance imaging (MRI) since the laser energy delivery system does not interfere with MRI. This report deals with mathematical modeling and analysis of laser coagulation of tissue. This model is intended for "real-time" analysis of magnetic resonance images obtained during the coagulation process to guide clinical treatment. A mathematical model is developed to simulate the thermal response of tissue to a laser light heating source. For fast simulation, an approximate solution of the thermal model is used to predict the dynamics of temperature distribution and tissue damage induced by a laser energy line source. The validity of these simulations is tested by comparison with MRI-based temperature data acquired from in vivo experiments in rabbits. The model-simulated temperature distribution and predicted lesion dynamics correspond closely with MRI-based data. These results demonstrate the potential for using this combination of fast modeling and MRI technologies during laser heating of tissue for online prediction of tumor lesion size during laser heating.     

Experimental evaluation of an elastic foundation model to predict contact pressures in knee replacements

Computational wear prediction is an attractive concept for evaluating new total knee replacement designs prior to physical testing and implementation. An important hurdle to such technology is the lack of in vivo contact pressure predictions. To address this issue, this study evaluates a computationally efficient simulation approach that combines the advantages of rigid and deformable body modeling. The hybrid method uses rigid body dynamics to predict body positions and orientations and elastic foundation theory to predict contact pressures between general three-dimensional surfaces. To evaluate the method, we performed static pressure experiments with a commercial knee implant in neutral alignment using flexion angles of 0, 30, 60, and 90 degrees and loads of 750, 1500, 2250, and 3000N. Using manufacturer CAD geometry for the same implant, an elastic foundation model with linear or nonlinear polyethylene material properties was implemented within a commercial multibody dynamics software program. The model's ability to predict experimental peak and average contact pressures simultaneously was evaluated by performing dynamic simulations to find the static configuration. Both the linear and nonlinear material models predicted the average contact pressure data well, while only the linear material model could simultaneously predict the trends in the peak contact pressure data. This novel modeling approach is sufficiently fast and accurate to be used in design sensitivity and optimization studies of knee implant mechanics and ultimately wear.     

QCT/FEA predictions of femoral stiffness are strongly affected by boundary condition modeling

Quantitative computed tomography-based finite element models of proximal femora must be validated with cadaveric experiments before using them to assess fracture risk in osteoporotic patients. During validation, it is essential to carefully assess whether the boundary condition (BC) modeling matches the experimental conditions. This study evaluated proximal femur stiffness results predicted by six different BC methods on a sample of 30 cadaveric femora and compared the predictions with experimental data. The average stiffness varied by 280% among the six BCs. Compared with experimental data, the predictions ranged from overestimating the average stiffness by 65% to underestimating it by 41%. In addition, we found that the BC that distributed the load to the contact surfaces similar to the expected contact mechanics predictions had the best agreement with experimental stiffness. We concluded that BC modeling introduced large variations in proximal femora stiffness predictions.     

An In Vivo Experimental Validation of a Computational Model of Human Foot

Reliable computational foot models offer an alternative means to enhance knowledge on the biomechanics of human foot.Model validation is one of the most critical aspects of the entire foot modeling and analysis process.This paper presents an invivo experiment combining motion capture system and plantar pressure measure platform to validate a three-dimensional finiteelement model of human foot.The Magnetic Resonance Imaging(MRI)slices for the foot modeling and the experimental datafor validation were both collected from the same volunteer subject.The validated components included the comparison of staticmodel predictions of plantar force,plantar pressure and foot surface deformation during six loading conditions,to equivalentmeasured data.During the whole experiment,foot surface deformation,plantar force and plantar pressure were recorded simultaneouslyduring six different loaded standing conditions.The predictions of the current FE model were in good agreementwith these experimental results.     

MRI-based finite element modeling of facial mimics: a case study on the paired zygomaticus major muscles

Finite element simulation of facial mimics provides objective indicators about soft tissue functions for improving diagnosis, treatment and follow-up of facial disorders. There is a lack of in vivo experimental data for model development and validation. In this study, the contribution of the paired Zygomaticus Major (ZM) muscle contraction on the facial mimics was investigated using in vivo experimental data derived from MRI. Maximal relative differences of 7.7% and 37% were noted between MRI-based measurements and numerical outcomes for ZM and skin deformation behaviors respectively. This study opens a new direction to simulate facial mimics with in vivo data.     

Comparison of MRI-based estimates of articular cartilage contact area in the tibiofemoral joint

Knee osteoarthritis (OA) detrimentally impacts the lives of millions of older Americans through pain and decreased functional ability. Unfortunately, the pathomechanics and associated deviations from joint homeostasis that OA patients experience are not well understood. Alterations in mechanical stress in the knee joint may play an essential role in OA; however, existing literature in this area is limited. The purpose of this study was to evaluate the ability of an existing magnetic resonance imaging (MRI)-based modeling method to estimate articular cartilage contact area in vivo. Imaging data of both knees were collected on a single subject with no history of knee pathology at three knee flexion angles. Intra-observer reliability and sensitivity studies were also performed to determine the role of operator-influenced elements of the data processing on the results. The method's articular cartilage contact area estimates were compared with existing contact area estimates in the literature. The method demonstrated an intra-observer reliability of 0.95 when assessed using Pearson's correlation coefficient and was found to be most sensitive to changes in the cartilage tracings on the peripheries of the compartment. The articular cartilage contact area estimates at full extension were similar to those reported in the literature. The relationships between tibiofemoral articular cartilage contact area and knee flexion were also qualitatively and quantitatively similar to those previously reported. The MRI-based knee modeling method was found to have high intra-observer reliability, sensitivity to peripheral articular cartilage tracings, and agreeability with previous investigations when using data from a single healthy adult. Future studies will implement this modeling method to investigate the role that mechanical stress may play in progression of knee OA through estimation of articular cartilage contact area.     

A finite element beam-model for efficient simulation of large-scale porous structures

This paper presents a new method for the generation of a beam finite element (FE) model from a three-dimensional (3D) data set acquired by micro-computed tomography (micro-CT). This method differs from classical modeling of trabecular bone because it models a specific sample only and differs from conventional solid hexahedron element-based FE approaches in its computational efficiency. The stress-strain curve, characterizing global mechanical properties of a porous structure, could be well predicted (R(2)=0.92). Furthermore, validation of the method was achieved by comparing local displacements of element nodes with the displacements directly measured by time-lapsed imaging methods of failure, and these measures were in good agreement. The presented model is a first step in modeling specific samples for efficient strength analysis by FE modeling. We believe that with upcoming high-resolution in-vivo imaging methods, this approach could lead to a novel and accurate tool in the risk assessment for osteoporotic fractures.     

Quantitative detection of cartilage surfaces and ligament geometry of the wrist using an imaging cryomicrotome system

Biomechanical models may aid in improving diagnosis and treatment of wrist joint disorders. As input, geometrical information is required for model development. Previous studies acquired some elements of the average wrist joint geometry. However, there is a close geometric functional match between articulating surfaces and ligament geometry. Therefore, biomechanical models need to be fed with the geometric data of individual joints. This study is aimed at acquiring geometric data of cartilage surfaces and ligaments from individual wrist joints by using a cryomicrotome imaging system and the evaluation of inter- and intra-observer variability of the data.The 3D geometry of 30 cartilage surfaces and 15 ligaments in three cadaver wrists was manually detected and quantitatively reconstructed. The inter- and intra-observer variability of the cartilage surface detection was 0.14 and 0.19 mm, respectively. For the position of the radius attachment of the dorsal radiocarpal ligament (DRC), the observer variations were 0.12 and 0.65 mm, for intra-/inter-observer, respectively. For the DRC attachment on the triquetrum, the observer variations were 0.22 and 1.19 mm.Anatomic reconstruction from 3D cryomicrotome images offer a method to obtain unique geometry data of the entire wrist joint for modeling purposes.     

Effect of subtalar arthroereisis on the tibiotalar contact characteristics in a cadaveric flatfoot model

Previous studies reported the effect of flatfoot deformity on tibiotalar joint contact characteristics. The lateral shift of the load which occurred in flatfeet may be responsible for degenerative changes in the ankle joint. The purpose was to assess the pattern of joint contact stress of the tibiotalar joint in intact, flat, and corrected specimens with subtalar arthroereisis. Seven fresh-frozen cadaver specimens were studied in the intact and flat-footed condition after transection of ligaments which support the medial arch. Ankle joint contact stress and plantar pressure patterns were determined from a capacitive pressure sensor inserted in the tibiotalar joint and a pressure distribution platform when the specimens were axially loaded in simulated mid-stance. Contact pressure was also assessed after subtalar arthroereisis with a 12 mm Kalix implant for correction of the flatfoot deformity. A maximum contact pressure of 1414.2±319.2 kPa was recorded in the middle-medial region in intact specimens. Flatfoot caused a lateral shift in the pressure distribution (p<0.05). In the flat specimens, the maximum contact pressure of 1394.7 8±470.5 kPa was in the anterior-central region. After subtalar arthroereisis with the Kalix implant for correction of the flatfoot deformity a maximum contact pressure of 1323.3±497.5 kPa was observed in the middle-lateral region. In a cadaver model, subtalar arthroereisis with Kalix implant failed to restore a normal intraarticular ankle joint pressure pattern. Further interventions should be considered to restore a normal pressure pattern.