Abstract
Introduction
Geometric variations of the hip joint can give rise to repeated abnormal contact between the femur and acetabular rim, resulting in cartilage and labrum damage. Population-based geometric parameterisation can facilitate the flexible and automated in silico generation of a range of clinically relevant hip geometries, allowing the position and size of cams to be defined precisely in three dimensions. This is advantageous compared to alpha angles, which are unreliable for stratifying populations by cam type. Alpha angles provide an indication of cam size in a single two-dimensional view, and high alpha angles have been observed in asymptomatic individuals.
Parametric geometries can be developed into finite element models to assess the potential effects of morphological variations in bone on soft tissue strains. The aim of this study was to demonstrate the capabilities of our parameterisation research tool by assessing impingement severity resulting from a range of parametrically varied femoral and acetabular geometries.
Methods
Custom made MATLAB (MathWorks) and Python codes[1] were used to generate bone surfaces, which were developed into finite element models in Abaqus (SIMULIA). Parametric femoral surfaces were defined by a spherical proximal head and ellipse sections through the neck/cam region. This method produced surfaces that were well fitted to bone geometry segmented from CT scans of cam patients and capable of producing trends in results similar to those found using segmented models. A simplified spherical geometry, including the labrum and acetabular cartilage, represented the acetabulum.
Femoral parameters were adjusted to define relevant variations in cam size and position. Two radii (small and large cams) and two positions (anterior and superior cams) were defined resulting in four models. Alpha angles of these parametric femurs were measured in an anterior-posterior view and a cross-table lateral view using ImageJ (NIH). A further model was developed using a femur with a medium cam size and position, and the level of acetabular coverage and labrum length were varied.
Bones were modelled as rigid bodies and soft tissues were modelled as transversely isotropic linearly elastic materials. With the acetabulum fully constrained in all cases, the femurs were constrained in translation and rotated to simulate flexion followed by internal rotation to cause impingement against the labrum.
Results and Discussion
Models generated using the parametric approach showed that potential for tissue damage, indicated through local strain, was not predicted by measured alpha angle, but resulted from cam extent and position as defined by the ellipses. When variations were made to the acetabular rim, an increase in bone coverage had the greatest effect on impingement severity, indicated by strain in the cartilage labral-junction. An increase in labral length increased labral displacement, but had less effect on cartilage-labral strain. Patient specific models currently require full image segmentation, but there is potential to further develop these parametric methods to assess likely impingement severity based on a series of measures of the neck and acetabulum when three-dimensional imaging of patients is available.
