Force Closure Mechanism Modeling for Musculoskeletal Multibody Simulation

Patella height influences patellofemoral contact and kinematics following cruciate-retaining total knee replacement

The role of patella height is discussed controversially in total knee arthroplasty (TKA). Therefore, this computational study aims to systematically analyze the biomechanical effect of different patella heights on patellofemoral (PF) forces and kinematics after cruciate-retaining (CR) TKA. We implemented a CR bicondylar TKA with a dome patellar button in a validated dynamic musculoskeletal multibody model of a male human knee joint. Retropatellar dynamics (contact force [N], shear force [N], patellar shift [mm], tilt [°], and rotation [°]) were evaluated during dual-limb squat motion (flexion from 0° to 90°) with simulated active muscle forces and the effects of different patella heights (Blackburne-Peel [BP] ratio of 0.39, 0.49, 0.65, 0.85, 1.01, and 1.1 were systematically examined). As active knee flexion increased, PF contact force also increased. Patella alta (BP = 1.1) resulted in higher PF contact forces compared to normal patella height (BP = 0.65) by up to 16%. Contrarily, patella baja was associated with decreased PF forces by 7%. Compared to patella baja (BP = 0.39), patella alta (BP = 1.1) considerably increased the contact force by up to 25%. Different patellar heights mainly affected PF shear forces during early knee flexion. Concerning PF kinematics, patella alta (BP = 1.1) yielded a greater lateral tilt of more than 4° and higher patellar rotation by up to 3° during deep knee flexion, compared to normal patella height (BP = 0.65). Our computational study indicates that patella alta is associated with the highest PF contact and shear force after the implantation of a CR bicondylar TKA. This should be considered in PF disorders following TKA.

Computer-based analysis of different component positions and insert thicknesses on tibio-femoral and patello-femoral joint dynamics after cruciate-retaining total knee replacement

BACKGROUND

Positioning of the implant components and tibial insert thickness constitute critical aspects of total knee replacement (TKR) that influence the postoperative knee joint dynamics. This study aimed to investigate the impact of implant component positioning (anterior-posterior and medio-lateral shift) and varying tibial insert thickness on the tibio-femoral (TF) and patello-femoral (PF) joint kinematics and contact forces after cruciate-retaining (CR)-TKR.

METHOD

A validated musculoskeletal multibody simulation (MMBS) model with a fixed-bearing CR-TKR during a squat motion up to 90° knee flexion was deployed to calculate PF and TF joint dynamics for varied implant component positions and tibial insert thicknesses. Evaluation was performed consecutively by comparing the respective knee joint parameters (e.g. contact force, quadriceps muscle force, joint kinematics) to a reference implant position.

RESULTS

The PF contact forces were mostly affected by the anterior-posterior as well as medio-lateral positioning of the femoral component (by 3 mm anterior up to 31 % and by 6 mm lateral up to 14 %). TF contact forces were considerably altered by tibial insert thickness (24 % in case of + 4 mm increase) and by the anterior-posterior position of the femoral component (by 3 mm posterior up to 16 %). Concerning PF kinematics, a medialised femoral component by 6 mm increased the lateral patellar tilt by more than 5°.

CONCLUSIONS

Our results indicate that regarding PF kinematics and contact forces the positioning of the femoral component was more critical than the tibial component. The positioning of the femoral component in anterior-posterior direction on and PF contact force was evident. Orthopaedic surgeons should strictly monitor the anterior-posterior as well as the medio-lateral position of the femoral component and the insert thickness.

The Influence of Mathematical Definitions on Patellar Kinematics Representations

A correlation between patellar kinematics and anterior knee pain is widely accepted. However, there is no consensus on how they are connected or what profile of patellar kinematics would minimize anterior knee pain. Nevertheless, answering this question by merging existing studies is further complicated by the variety of ways to describe patellar kinematics. Therefore, this study describes the most frequently used conventions for defining patellar kinematics, focusing on the rotations. The similarities and differences between the Cardan sequences and angles calculated by projecting axes are analyzed. Additionally, a tool is provided to enable the conversion of kinematic data between definitions in different studies. The choice of convention has a considerable impact on the absolute values and the clinical characteristics of the patello-femoral angles. In fact, the angles that result from using different mathematical conventions to describe a given patello-femoral rotation from our analyses differ up to a Root Mean Squared Error of 111.49° for patellar flexion, 55.72° for patellar spin and 35.39° for patellar tilt. To compare clinical kinematic patello-femoral results, every dataset must follow the same convention. Furthermore, researchers should be aware of the used convention’s implications to ensure reproducibility when interpreting and comparing such data.

Musculoskeletal Multibody Simulation Analysis on the Impact of Patellar Component Design and Positioning on Joint Dynamics after Unconstrained Total Knee Arthroplasty

Patellofemoral (PF) disorders are considered a major clinical complication after total knee replacement (TKR). Malpositioning and design of the patellar component impacts knee joint dynamics, implant fixation and wear propagation. However, only a limited number of studies have addressed the biomechanical impact of the patellar component on PF dynamics and their results have been discussed controversially. To address these issues, we implemented a musculoskeletal multibody simulation (MMBS) study for the systematical analysis of the patellar component’s thickness and positioning on PF contact forces and kinematics during dynamic squat motion with virtually implanted unconstrained cruciate-retaining (CR)-TKR. The patellar button thickness clearly increased the contact forces in the PF joint (up to 27%). Similarly, the PF contact forces were affected by superior–inferior positioning (up to 16%) and mediolateral positioning (up to 8%) of the patellar button. PF kinematics was mostly affected by the mediolateral positioning and the thickness of the patellar component. A medialization of 3 mm caused a lateral patellar shift by up to 2.7 mm and lateral patellar tilt by up to 1.6°. However, deviations in the rotational positioning of the patellar button had minor effects on PF dynamics. Aiming at an optimal intraoperative patellar component alignment, the orthopedic surgeon should pay close attention to the patellar component thickness in combination with its mediolateral and superior–inferior positioning on the retropatellar surface. Our generated MMBS model provides systematic and reproducible insight into the effects of patellar component positioning and design on PF dynamics and has the potential to serve as a preoperative analysis tool.

Estimating the Instantaneous Screw Axis and the Screw Axis Invariant Descriptor of Motion by Means of Inertial Sensors: An Experimental Study with a Mechanical Hinge Joint and Comparison to the Optoelectronic System

The motion of a rigid body can be represented by the instantaneous screw axis (ISA, also known as the helical axis). Recently, an invariant representation of motion based on the ISA, namely, the screw axis invariant descriptor (SAID), was proposed in the literature. The SAID consists of six scalar features that are independent from the coordinate system chosen to represent the motion. This method proved its usefulness in robotics; however, a high sensitivity to noise was observed. This paper aims to explore the performance of inertial sensors for the estimation of the ISA and the SAID for a simple experimental setup based on a hinge joint. The free swing motion of the mechanical hinge was concurrently recorded by a marker-based optoelectronic system (OS) and two magnetic inertial measurement units (MIMUs). The ISA estimated by the MIMU was more precise, while the OS was more accurate. The mean angular error was ≈2.2° for the OS and was ≈4.4° for the MIMU, while the mean standard deviation was ≈2.3° for the OS and was ≈0.2° for the MIMU. The SAID features based on angular velocity were better estimated by the MIMU, while the features based on translational velocity were better estimated by the OS. Therefore, a combination of both measurements systems is recommended to accurately estimate the complete SAID.

Effect of surgical parameters on the biomechanical behaviour of bicondylar total knee endoprostheses - A robot-assisted test method based on a musculoskeletal model

The complicated interplay of total knee replacement (TKR) positioning and patient-specific soft tissue conditions still causes a considerable number of unsatisfactory outcomes. Therefore, we deployed a robot-assisted test method, in which a six-axis robot moved and loaded a bicondylar cruciate-retaining (CR)-TKR in a virtual lower extremity emulated by a musculoskeletal multibody model. This enabled us to systematically analyse the impact of the posterior cruciate ligament (PCL), tibial slope, and tibial component rotation on TKR function while considering the physical implant components and physiological-like conditions during dynamic motions. The PCL resection yielded a decrease of femoral rollback by 4.5 mm and a reduction of tibiofemoral contact force by 50 N. A reduced tibial slope led to an increase of tibiofemoral contact force by about 170 N and a decrease of femoral rollback up to 1.7 mm. Although a higher tibial slope reduced the contact force, excessive tibial slopes should be avoided to prevent joint instability. Contrary to an external rotation of the tibial component, an internal rotation clearly increased the contact force and lateral femoral rollback. Our data contribute to improved understanding the biomechanics of TKRs and show the capabilities of the robot-assisted test method based on a musculoskeletal multibody model as a preoperative planning tool.

Neuro-musculoskeletal flexible multibody simulation yields a framework for efficient bone failure risk assessment

Fragility fractures are a major socioeconomic problem. A non-invasive, computationally-efficient method for the identification of fracture risk scenarios under the representation of neuro-musculoskeletal dynamics does not exist. We introduce a computational workflow that integrates modally-reduced, quantitative CT-based finite-element models into neuro-musculoskeletal flexible multibody simulation (NfMBS) for early bone fracture risk assessment. Our workflow quantifies the bone strength via the osteogenic stresses and strains that arise due to the physiological-like loading of the bone under the representation of patient-specific neuro-musculoskeletal dynamics. This allows for non-invasive, computationally-efficient dynamic analysis over the enormous parameter space of fracture risk scenarios, while requiring only sparse clinical data. Experimental validation on a fresh human femur specimen together with femur strength computations that were consistent with literature findings provide confidence in the workflow: The simulation of an entire squat took only 38 s CPU-time. Owing to the loss (16% cortical, 33% trabecular) of bone mineral density (BMD), the strain measure that is associated with bone fracture increased by 31.4%; and yielded an elevated risk of a femoral hip fracture. Our novel workflow could offer clinicians with decision-making guidance by enabling the first combined in-silico analysis tool using NfMBS and BMD measurements for optimized bone fracture risk assessment.