Scan-Free and Fully Automatic Tracking of Native Knee Anatomy from Dynamic Stereo-Radiography with Statistical Shape and Intensity Models.

Kinematic tracking of native anatomy from stereo-radiography provides a quantitative basis for evaluating human movement. Conventional tracking procedures require significant manual effort and call for acquisition and annotation of subject-specific volumetric medical images. The current work introduces a framework for fully automatic tracking of native knee anatomy from dynamic stereo-radiography which forgoes reliance on volumetric scans. The method consists of three computational steps. First, captured radiographs are annotated with segmentation maps and anatomic landmarks using a convolutional neural network. Next, a non-convex polynomial optimization problem formulated from annotated landmarks is solved to acquire preliminary anatomy and pose estimates. Finally, a global optimization routine is performed for concurrent refinement of anatomy and pose. An objective function is maximized which quantifies similarities between masked radiographs and digitally reconstructed radiographs produced from statistical shape and intensity models. The proposed framework was evaluated against manually tracked trials comprising dynamic activities, and additional frames capturing a static knee phantom. Experiments revealed anatomic surface errors routinely below 1.0 mm in both evaluation cohorts. Median absolute errors of individual bone pose estimates were below 1.0 ∘ or mm for 15 out of 18 degrees of freedom in both evaluation cohorts. Results indicate that accurate pose estimation of native anatomy from stereo-radiography may be performed with significantly reduced manual effort, and without reliance on volumetric scans.

Fully automatic tracking of native knee kinematics from stereo-radiography with digitally reconstructed radiographs.

Precise measurement of joint-level motion from stereo-radiography facilitates understanding of human movement. Conventional procedures for kinematic tracking require significant manual effort and are time intensive. The current work introduces a method for fully automatic tracking of native knee kinematics from stereo-radiography sequences. The framework consists of three computational steps. First, biplanar radiograph frames are annotated with segmentation maps and key points using a convolutional neural network. Next, initial bone pose estimates are acquired by solving a polynomial optimization problem constructed from annotated key points and anatomic landmarks from digitized models. A semidefinite relaxation is formulated to realize the global minimum of the non-convex problem. Pose estimates are then refined by registering computed tomography-based digitally reconstructed radiographs to masked radiographs. A novel rendering method is also introduced which enables generating digitally reconstructed radiographs from computed tomography scans with inconsistent slice widths. The automatic tracking framework was evaluated with stereo-radiography trials manually tracked with model-image registration, and with frames which capture a synthetic leg phantom. The tracking method produced pose estimates which were consistently similar to manually tracked values; and demonstrated pose errors below 1.0 degree or millimeter for all femur and tibia degrees of freedom in phantom trials. Results indicate the described framework may benefit orthopaedics and biomechanics applications through acceleration of kinematic tracking.

Kinematic Performance of Medial Pivot Total Knee Arthroplasty.

BACKGROUND: Total knee arthroplasty (TKA) implants have continued to evolve to accommodate new understandings of knee mechanics. The medial-pivot implant is a newer design, which is intended to limit anterior-posterior translation in the medial compartment while allowing lateral compartment translation. However, evidence for a generalized medial-pivot characteristic across all activities is limited. The purpose of the study was to quantify and compare in vivo knee joint kinematics using high-speed stereo radiography during activities of daily living in patients who have undergone a TKA with a cruciate sacrificing medial-pivot implant to age-matched and sex-matched native controls. METHODS: Fifteen participants (7 patients, 4 women, mean age 70 years and 8 nonsymptomatic controls, 4 women, mean age 64 years) performed 6 functional tasks in high-speed stereo radiography: deep-knee lunge, chair rise, step down, gait, gait with 90° turn, and seated knee extension. Translational differences between groups (surgical versus control) were assessed for the medial and lateral condyle, while pivot location was normalized to subject-specific tibial plateau geometry. RESULTS: The surgical cohort displayed a more constrained medial condyle that provided greater stability of the medial compartment and did not result in the paradoxical anterior translation at mid-flexion angles during weight-bearing activities, but was associated with less condylar translation than native knees. Additionally, the transverse tibial pivot location occurs most commonly in the middle third of the tibial plateau and secondarily on the medial third. CONCLUSIONS: Some variability in pivot location occurs between activities and is more in nonsymptomatic, native knee controls.

Fully automatic tracking of native glenohumeral kinematics from stereo-radiography.

The current work introduces a system for fully automatic tracking of native glenohumeral kinematics in stereo-radiography sequences. The proposed method first applies convolutional neural networks to obtain segmentation and semantic key point predictions in biplanar radiograph frames. Preliminary bone pose estimates are computed by solving a non-convex optimization problem with semidefinite relaxations to register digitized bone landmarks to semantic key points. Initial poses are then refined by registering computed tomography-based digitally reconstructed radiographs to captured scenes, which are masked by segmentation maps to isolate the shoulder joint. A particular neural net architecture which exploits subject-specific geometry is also introduced to improve segmentation predictions and increase robustness of subsequent pose estimates. The method is evaluated by comparing predicted glenohumeral kinematics to manually tracked values from 17 trials capturing 4 dynamic activities. Median orientation differences between predicted and ground truth poses were 1.7∘ and 8.6∘ for the scapula and humerus, respectively. Joint-level kinematics differences were less than 2∘ in 65%, 13%, and 63% of frames for XYZ orientation DoFs based on Euler angle decompositions. Automation of kinematic tracking can increase scalability of tracking workflows in research, clinical, or surgical applications.

Modeling of active skeletal muscles: a 3D continuum approach incorporating multiple muscle interactions.

Skeletal muscles have a highly organized hierarchical structure, whose main function is to generate forces for movement and stability. To understand the complex heterogeneous behaviors of muscles, computational modeling has advanced as a non-invasive approach to evaluate relevant mechanical quantities. Aiming to improve musculoskeletal predictions, this paper presents a framework for modeling 3D deformable muscles that includes continuum constitutive representation, parametric determination, model validation, fiber distribution estimation, and integration of multiple muscles into a system level for joint motion simulation. The passive and active muscle properties were modeled based on the strain energy approach with Hill-type hyperelastic constitutive laws. A parametric study was conducted to validate the model using experimental datasets of passive and active rabbit leg muscles. The active muscle model with calibrated material parameters was then implemented to simulate knee bending during a squat with multiple quadriceps muscles. A computational fluid dynamics (CFD) fiber simulation approach was utilized to estimate the fiber arrangements for each muscle, and a cohesive contact approach was applied to simulate the interactions among muscles. The single muscle simulation results showed that both passive and active muscle elongation responses matched the range of the testing data. The dynamic simulation of knee flexion and extension showed the predictive capability of the model for estimating the active quadriceps responses, which indicates that the presented modeling pipeline is effective and stable for simulating multiple muscle configurations. This work provided an effective framework of a 3D continuum muscle model for complex muscle behavior simulation, which will facilitate additional computational and experimental studies of skeletal muscle mechanics. This study will offer valuable insight into the future development of multiscale neuromuscular models and applications of these models to a wide variety of relevant areas such as biomechanics and clinical research.

Joint Track Machine Learning: An Autonomous Method of Measuring Total Knee Arthroplasty Kinematics From Single-Plane X-Ray Images.

BACKGROUND: Dynamic radiographic measurements of 3-dimensional (3-D) total knee arthroplasty (TKA) kinematics have provided important information for implant design and surgical technique for over 30 years. However, current methods of measuring TKA kinematics are too cumbersome, inaccurate, or time-consuming for practical clinical application. Even state-of-the-art techniques require human-supervision to obtain clinically reliable kinematics. Eliminating human supervision could potentially make this technology practical for clinical use. METHODS: We demonstrate a fully autonomous pipeline for quantifying 3D-TKA kinematics from single-plane radiographic imaging. First, a convolutional neural network (CNN) segmented the femoral and tibial implants from the image. Second, those segmented images were compared to precomputed shape libraries for initial pose estimates. Lastly, a numerical optimization routine aligned 3D implant contours and fluoroscopic images to obtain the final implant poses. RESULTS: The autonomous technique reliably produces kinematic measurements comparable to human-supervised measures, with root-mean-squared differences of less than 0.7 mm and 4° for our test data, and 0.8 mm and 1.7° for external validation studies. CONCLUSION: A fully autonomous method to measure 3D-TKA kinematics from single-plane radiographic images produces results equivalent to a human-supervised method, and may soon make it practical to perform these measurements in a clinical setting.

Evaluation of single-stage vision models for pose estimation of surgical instruments.

PURPOSE: Multiple applications in open surgical environments may benefit from adoption of markerless computer vision depending on associated speed and accuracy requirements. The current work evaluates vision models for 6-degree of freedom pose estimation of surgical instruments in RGB scenes. Potential use cases are discussed based on observed performance. METHODS: Convolutional neural nets were developed with simulated training data for 6-degree of freedom pose estimation of a representative surgical instrument in RGB scenes. Trained models were evaluated with simulated and real-world scenes. Real-world scenes were produced by using a robotic manipulator to procedurally generate a wide range of object poses. RESULTS: CNNs trained in simulation transferred to real-world evaluation scenes with a mild decrease in pose accuracy. Model performance was sensitive to input image resolution and orientation prediction format. The model with highest accuracy demonstrated mean in-plane translation error of 13 mm and mean long axis orientation error of 5[Formula: see text] in simulated evaluation scenes. Similar errors of 29 mm and 8[Formula: see text] were observed in real-world scenes. CONCLUSION: 6-DoF pose estimators can predict object pose in RGB scenes with real-time inference speed. Observed pose accuracy suggests that applications such as coarse-grained guidance, surgical skill evaluation, or instrument tracking for tray optimization may benefit from markerless pose estimation.

Knee pivot location in asymptomatic older adults.

Representative data of asymptomatic, native-knee kinematics is important when studying changes in knee function across the lifespan. High-speed stereo radiography (HSSR) provides a reliable measure of knee kinematics to <1 mm of translation and 1° of rotation, but studies often have limited statistical power to make comparisons between groups or measure the contribution of individual variability. The purpose of this study is to examine in vivo condylar kinematics to quantify the transverse center-of-rotation, or pivot, location across the flexion range and challenge the medial-pivot paradigm in asymptomatic knee kinematics. We quantified the pivot location during supine leg press, knee extension, standing lunge, and gait for 53 middle-aged and older adults (27 men; 26 women: 50.8 ± 7.0 yrs, 1.75 ± 0.1 m, 79.1 ± 15.4 kg). A central- to medial-pivot location was identified for all activities with increased knee flexion associated with posterior translation of the center-of-rotation. The association between knee angle and anterior-posterior center-of-rotation location was not as strong as the relation between medial-lateral and anterior-posterior location, excluding gait. The Pearson's correlation for gait was stronger between knee angle and anterior-posterior center-of-rotation location (P < 0.001) than medial-lateral and anterior-posterior location (P = 0.0122). Individual variability accounted for a measurable proportion in variance explained of center-of-rotation location. Unique to gait, the lateral translation of center-of-rotation location resulted in the anterior translation of center-of-rotation at <10° knee flexion. Furthermore, no association between vertical ground-reaction force and center-of-rotation was identified.

Computational comparison of medializing tibial tubercle osteotomy and trochleoplasty in patients with trochlear dysplasia.

Medial patellofemoral ligament reconstruction (MPFLR) has emerged as the procedure of choice for recurrent patellar dislocation. This addresses soft tissue injury but does not address underlying anatomic factors, including trochlear dysplasia, that are commonly present and increase risk of dislocation. Quantification of the stability offered by other surgical interventions, namely, medializing tibial tubercle osteotomy (mTTO) and trochleoplasty, with and without MPFLR, may provide insight for surgical choices in patients with trochlear dysplasia. We developed subject-specific finite element models based on magnetic resonance scans from a cohort of 20 patients with trochlear dysplasia and recurrent patellar dislocation. The objectives of this study were (1) to compare patella stability after mTTO and trochleoplasty procedures; (2) to evaluate whether it is necessary to perform an MPFLR in combination with the mTTO or trocheoplasty procedure; and (3) to quantify the robustness of patellar stability to variability in knee kinematics. Trochleoplasty performed better than mTTO at stabilizing the patella between 5° and 30° flexion. For both mTTO and trochleoplasty procedures, it was beneficial to also perform MPFLR-inclusion of MPFLR halved the magnitude of patellar laxity predicted in the simulations. Simulations that did not include any medial patellofemoral ligament restraint were also more sensitive to variation in tibiofemoral internal-external kinematics.

Impact of patient, surgical, and implant design factors on predicted tray-bone interface micromotions in cementless total knee arthroplasty.

Micromotion magnitudes exceeding 150 µm may prevent bone formation and limit fixation after cementless total knee arthroplasty (TKA). Many factors influence the tray-bone interface micromotion but the critical parameters and sensitivities are less clear. In this study, we assessed the impacts of surgical (tray alignment, tibial coverage, and resection surface preparation), patient (bone properties and tibiofemoral kinematics), and implant design (tray feature and surface friction) factors on tray-bone interface micromotions during a series of activities of daily living. Micromotion was estimated via three previously validated implant-bone finite element models and tested under gait, deep knee bending, and stair descent loads. Overall, the average micromotion across the tray-bone cementless contact interface ranged from 9.3 to 111.4 µm, and peak micromotion was consistently found along the anterior tray edge. Maximizing tibial coverage above a properly sized tibial tray (an average of 12.3% additional area) had minimal impact on micromotion. A 1 mm anterior tray alignment change reduced the average micromotion by an average of 16.1%. Two-degree tibial angular resection errors reduced the area for bone ingrowth up to 48.1%. Differences on average micromotion from ±25% changes in bone moduli were up to 75.5%. A more posterior tibiofemoral contact due to additional 100 N posterior force resulted in an average of 79.3% increase on average micromotion. Overall, careful surgical technique, patient selection, and controlling kinematics through articular design all contribute meaningfully to minimizing micromotion in cementless TKA, with centralizing the load transfer to minimize the resulting moment at the anterior tray perimeter a consistent theme.

Impact of bone health on the mechanics of plate fixation for Vancouver B1 periprosthetic femoral fractures.

BACKGROUND: Condyle-spanning plate-screw constructs have shown potential to lower the risks of femoral refractures after the healing of a primary Vancouver type B1 periprosthetic femoral fracture. Limited information exists to show how osteoporosis (a risk factor for periprosthetic femoral fractures) may affect the plate fixation during activities of daily living. METHODS: Using total hip arthroplasty and plate-implanted finite element models of three osteoporotic femurs, this study simulated physiological loads of three activities of daily living, as well as osteoporosis associated muscle weakening, and compared the calculated stress/strain, load transfer and local stiffness with experimentally validated models of three healthy femurs. Two plating systems and two construct lengths (a diaphyseal construct and a condyle-spanning construct) were modeled. FINDINGS: Osteoporotic femurs showed higher bone strain (21.9%) and higher peak plate stress (144.3%) as compared with healthy femurs. Compared with shorter diaphyseal constructs, condyle-spanning constructs of two plating systems reduced bone strains in both healthy and osteoporotic femurs (both applying 'the normal' and 'the weakened muscle forces') around the most distal diaphyseal screw and in the distal metaphysis, both locations where secondary fractures are typically reported. The lowered resultant compressive force and the increased local compressive stiffness in the distal diaphysis and metaphysis may be associated with strain reductions via condyle-spanning constructs. INTERPRETATION: Strain reductions in condyle-spanning constructs agreed with the clinically reported lowered risks of femoral refractures in the distal diaphysis and metaphysis. Multiple condylar screws may mitigate the concentrated strains in the lateral condyle, especially in osteoporotic femurs.

Drivers of initial stability in cementless TKA: Isolating effects of tibiofemoral conformity and fixation features.

The initial fixation of cementless tibial trays after total knee arthroplasty is critical to ensure bony ingrowth and long-term fixation. Various fixed-bearing implant designs that utilize different fixation features, surface coatings, and bony preparations to facilitate this initial stability are currently used clinically. However, the role of tibiofemoral conformity and the effect of different tray fixation features on initial stability are still unclear. This study assessed the implant stability of two TKA designs during a series of simulated daily activities including experimental testing and corresponding computational models. Tray-bone interface micromotions and the porous area ideal for bone ingrowth were investigated computationally and compared between the two designs. The isolated effect of femoral-insert conformity and fixation features on the micromotion was examined separately by virtually exchanging design features. The peak interface micromotions predicted were at least 47% different for the two designs, which was a combined result of different femoral-insert conformity (contributed 79% of the micromotion difference) and fixation features (21%). A more posterior femoral-insert contact due to lower tibiofemoral conformity in a force-controlled simulation significantly increased the micromotion and reduced the surface area ideal for bone ingrowth. The maximum difference in peak micromotions caused by only changing the fixation features was up to 33%. Overall, the moment arm from the insert articular contact point to the anterolateral tray perimeter was the primary factor correlated to peak and average micromotion. Our results indicated that tray-bone micromotion could be minimized by centralizing the load transfer and optimizing the fixation features.

Supine leg press as an alternative to standing lunge in high-speed stereo radiography.

The standing lunge is an activity commonly used to quantify in-vivo knee kinematics with fluoroscopy. The ability to perform the standing lunge varies between subjects and can necessitate movement accommodations to successfully complete the desired range of motion. We proposed a supine leg press as an alternative to the standing lunge that aimed to provide a similar evaluation of knee motion while increasing the measured range of motion. Tibiofemoral kinematics of 53 non-symptomatic adults (27 men, 26 women, 50.8 ± 7.0 yrs.) were calculated from the tracked high-speed stereo radiography (HSSR) images for supine leg press and standing lunge using CT-segmented bony geometries of the right lower limb. The supine leg press proved to be a useful alternative to the standing lunge while providing 46.2° greater range of motion in knee flexion. The difference in angle-matched kinematics across a 100° flexion range between the leg press and lunge was 0.70° in varus-valgus rotation, 1.5° in internal-external rotation, 1.0 mm in medial-lateral translation, 2.3 mm in anterior-posterior translation, and 0.46 mm in superior-inferior translation for men. The angle-matched difference for women across 100° was 0.58° in varus-valgus rotation, 2.4° internal-external rotation, 0.70 mm medial-lateral translation, 2.1 mm anterior-posterior translation, and 0.78 mm superior-inferior translation. The similar kinematics, while having a greater range of motion, and control of the applied load makes the supine leg press an alternative for quantifying in-vivo knee kinematics.

Simplified Mechanical Tests Can Simulate Physiological Mechanics of a Fixation Construct for Periprosthetic Femoral Fractures.

Plate fractures after fixation of a Vancouver Type B1 periprosthetic femoral fracture (PFF) are difficult to treat and could lead to severe disability. However, due to the lack of direct measurement of in vivo performance of the PFF fixation construct, it is unknown whether current standard mechanical tests or previous experimental and computational studies have appropriately reproduced the in vivo mechanics of the plate. To provide a basis for the evaluation and development of appropriate mechanical tests for assessment of plate fracture risk, this study applied loads of common activities of daily living (ADLs) to implanted femur finite element (FE) models with PFF fixation constructs with an existing or a healed PFF. Based on FE simulated plate mechanics, the standard four-point-bend test adequately matched the stress state and the resultant bending moment in the plate as compared with femur models with an existing PFF. In addition, the newly developed constrained three-point-bend tests were able to reproduce plate stresses in models with a healed PFF. Furthermore, a combined bending and compression cadaveric test was appropriate for risk assessment including both plate fracture and screw loosening after the complete healing of PFF. The result of this study provides the means for combined experimental and computational preclinical evaluation of PFF fixation constructs.

Effect of intraoperative treatment options on hip joint stability following total hip arthroplasty.

Dislocation remains the leading indication for revision of total hip arthroplasty (THA). The objective of this study was to use a computational model to compare the overall resistance to both anterior and posterior dislocation for the available THA constructs commonly considered by surgeons attempting to produce a stable joint. Patient-specific musculoskeletal models of THA patients performing activities consistent with anterior and posterior dislocation were developed to calculate joint contact forces and joint positions used for simulations of dislocation in a finite element model of the implanted hip that included an experimentally calibrated hip capsule representation. Dislocations were then performed with consideration of offset using +5 and +9 offset, iteratively with three lipped liner variations in jump distance (10°, 15°, and 20° lips), a size 40 head, and a dual-mobility construct. Dislocation resistance was quantified as the moment required to dislocate the hip and the integral of the moment-flexion angle (dislocation energy). Increasing head diameter increased resistive moment on average for anterior and posterior dislocation by 22% relative to a neutral configuration. A lipped liner resulted in increases in the resistive moment to posterior dislocation of 9%, 19%, and 47% for 10°, 15°, and 20° lips, a sensitivity of approximately 2.8 Nm/mm of additional jump distance. A dual-mobility acetabular design resulted in an average 38% increase in resistive moment and 92% increase in dislocation energy for anterior and posterior dislocation. A quantitative understanding of tradeoffs in the dislocation risk inherent to THA construct options is valuable in supporting surgical decision making.

Impact of periprosthetic femoral fracture fixation plating constructs on local stiffness, load transfer, and bone strains.

Secondary femoral fractures after the successful plate-screw fixation of a primary Vancouver type B1 periprosthetic femoral fracture (PFF) have been associated with the altered state of stress/strain in the femur as the result of plating. The laterally implanted condyle-spanning plate-screw constructs have shown promises clinically in avoiding secondary bone and implant failures as compared with shorter diaphyseal plates. Though the condyle-spanning plating has been hypothesized to avoid stress concentration in the femoral diaphysis through increasing the working length of the plate, biomechanical evidence is lacking on how plate length may impact the stress/strain state of the implanted femur. Through developing and experimentally validating finite element (FE) models of 3 cadaveric femurs, this study investigated the impact of plating on bone strains, load transfer and local stiffness, which were compared between FE models of 2 different plating systems that each had a diaphyseal configuration and a condyle-spanning configuration. Under simulated gait-loading, the condyle-spanning constructs of both plating systems were shown to lower the bone strains around the distal fixation screws (up to 24.8% reduction in maximum principal strain and 26.6% reduction in minimum principal strain) and in the distal metaphyseal shaft of the femur (up to 15.9% and 25.7% reductions in maximum and minimum principal strains, respectively), where secondary bone fractures have been typically reported. In the distal diaphyseal and metaphyseal shaft of femur, FE models of the condyle-spanning constructs were shown to increase the local compressive stiffness (up to 152.9% increases under simulated gait-loading) and decrease the transfer of compressive load (37.1% decreases under simulated gait-loading), which may be indicative of the lowered risks of bone damage.

Assessment of the primary stability of revision tibial trays augmented with a cementless sleeve in AORI Type III defects.

BACKGROUND: Previous studies have evaluated the initial stability of uncemented tibial components in revision total knee replacement (rTKR) in the presence of an Anderson Orthopaedic Research Institute (AORI) Type II tibial defect. This study sought to evaluate similar metrics in the severe Type III (AORI TIII) defects with combined uncemented stem and sleeve fixation, specifically, the effect of varying the stem's length and tibial canal engagement upon stability and bone strain. METHOD: Finite element models generated from the CT scans of 4 tibias with Type III defects were used to investigate the primary stability, in terms of the bone-implant composite peak micromotion (CPM) and microstrains (CPS), achieved after virtual implantations with and without stems. RESULTS: A stemless rTKR had increased metaphyseal CPM and CPS compared to all stemmed implants. Significant area of the bone supporting the stemless rTKR had CPS greater than bone yield (7000 µÎµ). Short engaging stems (≤150 mm construct length), could not achieve reliable engagement in the diaphysis (canal fill ≤ 50%), leading to insufficient reduction of CPS (≥5000 µÎµ). Longer engaging stems (170-220 mm construct length), were able to reliably engage the diaphysis (fill ratio ≥ 75%) resulting in CPS ≤ 5000 µÎµ. Although, non-engaging stems resulted in increased CPM and CPS compared to engaging stems, long non-engaging stems (170-220 mm construct) appeared to provide additional stability to the rTKR compared to stemless rTKR. CONCLUSION: The results indicate a likely correlation between uncemented stem engagement and metaphyseal CPS in Type III defects. Excessive strain within the supporting metaphyseal bone is likely to lead to rTKR migration and loosening.

Automatic tracking of healthy joint kinematics from stereo-radiography sequences.

Kinematic tracking of healthy joints in radiography sequences is frequently performed by maximizing similarities between computed perspective projections of 3D computer models and corresponding objects' appearances in radiographic images. Significant human effort associated with manual tracking presents a major bottleneck in biomechanics research methods and limits the scale of target applications. The current work introduces a method for fully-automatic tracking of tibiofemoral and patellofemoral kinematics in stereo-radiography sequences for subjects performing dynamic activities. The proposed method involves the application of convolutional neural networks for annotating radiographs and a multi-stage optimization pipeline for estimating bone pose based on information provided by neural net predictions. Predicted kinematics are evaluated by comparing against manually-tracked trends across 20 distinct trials. Median absolute differences below 1.5 millimeters or degrees for 6 tibiofemoral and 3 patellofemoral degrees of freedom demonstrate the utility of our approach, which improves upon previous semi-automatic methods by enabling end-to-end automation. Implementation of a fully-automatic pipeline for kinematic tracking will benefit evaluation of human movement by enabling large-scale studies of healthy knee kinematics.

Machine learning for rapid estimation of lower extremity muscle and joint loading during activities of daily living.

Joint contact and muscle forces estimated with musculoskeletal modeling techniques offer useful metrics describing movement quality that benefit multiple research and clinical applications. The expensive processing of laboratory data associated with generating these outputs presents challenges to researchers and clinicians, including significant time and expertise requirements that limit the number of subjects typically evaluated. The objective of the current study was to develop and compare machine learning techniques for rapid, data-driven estimation of musculoskeletal metrics from derived gait lab data. OpenSim estimates of patient joint and muscle forces during activities of daily living were simulated using laboratory data from 70 total knee replacement patients and used to develop 4 different machine learning algorithms. Trained machine learning models predicted both trend and magnitude of estimated joint contact (mean correlation coefficients ranging from 0.93 to 0.94 during gait) and muscle forces (mean correlation coefficients ranging from 0.83 to 0.91 during gait) based on anthropometrics, ground reaction forces, and joint angle data. Patient mechanics were accurately predicted by recurrent neural networks, even after removing dependence on key subsets of predictor features. The ability to quickly estimate patient mechanics from derived measurements of movement has the potential to broaden the impact of musculoskeletal modeling by enabling faster assessment in both clinical and research settings.

Development of axial compression and combined axial compression and torque loading configurations to reproduce strain in the implanted femur during activities of daily living.

Femoral strain is indicative of the potential for bone remodeling (strain energy density, SED) and periprosthetic femoral fracture (magnitude of principal strains) after total hip arthroplasty (THA). Previous modeling studies have evaluated femoral strains in THA-implanted femurs under gait loads including both physiological hip contact force and femoral muscle forces. However, experimental replication of the complex muscle forces during activities of daily living (ADLs) is difficult for in vitro assessment of femoral implant or fixation hardware. Alternatively, cadaveric tests using simplified loading configurations have been developed to assess post-THA bone mechanics, although no current studies have demonstrated simplified loading configurations used in mechanical tests may simulate the physiological femoral strains under ADL loads. Using an optimization approach integrated with finite element analysis, this study developed axial compression and combined axial compression and torque testing configurations for three common ADLs (gait, stair-descent and sit-to-stand) via matching the SED profile of the femur in THA-implanted models of three specimens. The optimized simplified-loading models showed good agreement in predicting bone remodeling stimuli (post-THA change in SED per unit mass) and fatigue regions as compared with the ADL-loading models, as well as other modeling and clinical studies. The optimized simplified test configurations can provide a physiological-loading based pre-clinical platform for the evaluation of implant/fixation devices of the femur.