Ward Bartels

Subject-specific hip geometry and hip joint centre location affects calculated contact forces at the hip during gait

Hip loading affects the development of hip osteoarthritis, bone remodelling and osseointegration of implants. In this study, we analyzed the effect of subject-specific modelling of hip geometry and hip joint centre (HJC) location on the quantification of hip joint moments, muscle moments and hip contact forces during gait, using musculoskeletal modelling, inverse dynamic analysis and static optimization. For 10 subjects, hip joint moments, muscle moments and hip loading in terms of magnitude and orientation were quantified using three different model types, each including a different amount of subject-specific detail: (1) a generic scaled musculoskeletal model, (2) a generic scaled musculoskeletal model with subject-specific hip geometry (femoral anteversion, neck-length and neck-shaft angle) and (3) a generic scaled musculoskeletal model with subject-specific hip geometry including HJC location. Subject-specific geometry and HJC location were derived from CT. Significant differences were found between the three model types in HJC location, hip flexion-extension moment and inclination angle of the total contact force in the frontal plane. No model agreement was found between the three model types for the calculation of contact forces in terms of magnitude and orientations, and muscle moments. Therefore, we suggest that personalized models with individualized hip joint geometry and HJC location should be used for the quantification of hip loading. For biomechanical analyses aiming to understand modified hip joint loading, and planning hip surgery in patients with osteoarthritis, the amount of subject-specific detail, related to bone geometry and joint centre location in the musculoskeletal models used, needs to be considered.

Muscle optimization techniques impact the magnitude of calculated hip joint contact forces

In musculoskeletal modelling, several optimization techniques are used to calculate muscle forces, which strongly influence resultant hip contact forces (HCF). The goal of this study was to calculate muscle forces using four different optimization techniques, i.e., two different static optimization techniques, computed muscle control (CMC) and the physiological inverse approach (PIA). We investigated their subsequent effects on HCFs during gait and sit to stand and found that at the first peak in gait at 15–20% of the gait cycle, CMC calculated the highest HCFs (median 3.9 times peak GRF (pGRF)). When comparing calculated HCFs to experimental HCFs reported in literature, the former were up to 238% larger. Both static optimization techniques produced lower HCFs (median 3.0 and 3.1 pGRF), while PIA included muscle dynamics without an excessive increase in HCF (median 3.2 pGRF). The increased HCFs in CMC were potentially caused by higher muscle forces resulting from co‐contraction of agonists and antagonists around the hip. Alternatively, these higher HCFs may be caused by the slightly poorer tracking of the net joint moment by the muscle moments calculated by CMC. We conclude that the use of different optimization techniques affects calculated HCFs, and static optimization approached experimental values best.

Computed tomography-based joint locations affect calculation of joint moments during gait when compared to scaling approaches

Hip joint moments are an important parameter in the biomechanical evaluation of orthopaedic surgery. Joint moments are generally calculated using scaled generic musculoskeletal models. However, due to anatomical variability or pathology, such models may differ from the patients anatomy, calling into question the accuracy of the resulting joint moments. This study aimed to quantify the potential joint moment errors caused by geometrical inaccuracies in scaled models, during gait, for eight test subjects. For comparison, a semi-automatic computed tomography (CT)-based workflow was introduced to create models with subject-specific joint locations and inertial parameters. 3D surface models of the femora and hemipelves were created by segmentation and the hip joint centres and knee axes were located in these models. The scaled models systematically located the hip joint centre (HJC) up to 33.6 mm too inferiorly. As a consequence, significant and substantial peak hip extension and abduction moment differences were recorded, with, respectively, up to 23.1% and 15.8% higher values in the image-based models. These findings reaffirm the importance of accurate HJC estimation, which may be achieved using CT- or radiography-based subject-specific modelling. However, obesity-related gait analysis marker placement errors may have influenced these results and more research is needed to overcome these artefacts.

Subject-specific geometrical detail rather than cost function formulation affects hip loading calculation*

Abstract This study assessed the relative importance of introducing an increasing level of medical image-based subject-specific detail in bone and muscle geometry in the musculoskeletal model, on calculated hip contact forces during gait. These forces were compared to introducing minimization of hip contact forces in the optimization criterion. With an increasing level of subject-specific detail, specifically MRI-based geometry and wrapping surfaces representing the hip capsule, hip contact forces decreased and were more comparable to contact forces measured using instrumented prostheses (average difference of 0.69 BW at the first peak compared to 1.04 BW for the generic model). Inclusion of subject-specific wrapping surfaces in the model had a greater effect than altering the cost function definition.

Towards automatic computer‐aided knee surgery by innovative methods for processing the femur surface model

The femoral shaft (FDA) and transepicondylar (TA), anterior–posterior (WL) and posterior condylar (PCL) axes are fundamental quantities in planning knee arthroplasty surgery. As an alternative to the TA, we introduce the anatomical flexion axis (AFA). Obtaining such axes from image data without any manual supervision remains a practical objective. We propose a novel method that automatically computes the axes of the distal femur by processing the femur mesh surface.

Changes in contact area characteristics of the ankle after a cartilage biopsy at the postero-medial rim of the talar dome

OBJECTIVE Study the changes in local and generalized biomechanical characteristics of the ankle joint, associated with a well defined cartilage biopsy at the postero-medial rim of the talar dome, to evaluate its safety. METHODS Ten cadaver ankles were (sub-) physiologically loaded pre- and post-biopsy; in neutral position, 10 degrees of plantar-flexion (PF) and 10 degrees of dorsi-flexion (DF). Fuji film was used as transducer. Qualitatively, the coverage of the biopsy by the tibial plafond, and changes in the shape of the footprint were analyzed. Quantitatively, the pressure profile plot, normalized-tibio-talar contact area and the centroid position of pressure were examined. Results were reported as a mean for all specimens, and as individual values for every single specimen as well. RESULTS Mean results did not show significant changes, but those of some single specimens did. The majority of those changes were in PF. Some occurred in N, and besides two exceptions none occurred in DF. Two specimens did not show any change. One specimen showed an isolated quantitative change. Seven specimens showed both qualitative and quantitative changes. However, all changes were of low-magnitude and contact stresses did not show any rebound effect. CONCLUSIONS Although biopsies at the postero-medial rim of the talar dome did not induce on average significant changes in quantitative contact characteristics, few specimens did show some alterations. Currently, the investigated biopsy site seems safe, but long term follow-up studies in patients are needed for confirmation.

Sensitivity analysis of hip joint centre estimation based on three-dimensional CT scans

In morphological analysis of the femur, the hip joint centre (HJC) is generally determined using a 3D model of the femoral head based on medical images. However, the portion of the image selected to represent the femoral head may influence the HJC. We determined if this influence invalidates the results of three HJC calculation methods, one of which we introduce here. To isolate femoral heads in cadaver CT images, thresholds were applied to the distance between femur and acetabulum models. The sensitivity of the HJC to these thresholds and the differences between methods were quantified. For thresholds between 6 and 9 mm and healthy hips, differences between methods were below 1 mm and all methods were insensitive to threshold changes. For higher thresholds, the fovea capitis femoris disturbed the HJC. In two deformed hips, the new method performed superiorly. We conclude that for normal hips all methods produce valid results.

Patient-Specific Reconstruction of Large Bone Defects: Clinical Success Due to an Integrated Bioengineering Workflow

Surgical reconstruction of major acetabular defects (Paprosky type IIIb) in view of total hip replacement is a challenge. Existing systems have a documented failure rate between 20% and 50% at five years post-operatively. Our first hypothesis is that this limited success is due to the fact that standard components have to be adapted and assembled intraoperatively to provide a patient solution. Hence our second hypothesis is that a fully personalized solution for acetabular reconstruction will provide an added value for the patient’s outcome.

Changes in the length of the sacrospinous and sacrotuberous ligaments induced by Salter osteotomy: a computer simulation

Severe developmental dysplasia of the hip may be treated by a Salter osteotomy. This osteotomy separates the origins of the sacrotuberous and sacrospinous ligaments from their insertions, thus changing the length of these ligaments. This study uses a multibody computer model in order to quantify this effect in relation to three different surgical parameters: lateral wedge height, wedge insertion angle and osteotomy height. The model, created in MSC.Adams, was based on medical images of an 8-year-old girl with severe dysplasia of both hips. Further data analysis was performed in the MATLAB matrix calculation environment. Every combination of surgical parameters that was tested, resulted in both ligaments getting shortened with the lateral wedge height having the strongest influence on these length changes. Thus, a Salter osteotomy may remove the stabilizing effects of these ligaments, but ligament strain does not appear to occur.