Corroboration of coupled musculoskeletal model and finite element predictions with in vivo RSA migration of an uncemented acetabular component.

While finite element (FE) models have been used extensively in orthopedic studies, validation of their outcome metrics has been limited to comparison against ex vivo testing. The aim of this study was to validate FE model predictions of the initial cup mechanical environment against patient-matched in vivo measurements of acetabular cup migration using radiostereometric analysis (RSA). Tailored musculoskeletal and FE models were developed using a combination of three-dimensional (3D) motion capture data and clinical computerized tomography (CT) scans for a cohort of eight individuals who underwent primary total hip replacement and were prospectively enrolled in an RSA study. FE models were developed to calculate the mean modulus of cancellous bone, composite peak micromotion (CPM), composite peak strain (CPS) and percentage area of bone ingrowth. The RSA cup migration at 3 months was used to corroborate the FE output metrics. Qualitatively, all FE-predicted metrics followed a similar rank order as the in vivo RSA 3D migration data. The two cases with the lowest predicted CPM (<20 µm), lowest CPS (<0.0041), and high bone modulus (>917 MPa) were confirmed to have the lowest in vivo RSA 3D migration (<0.14 mm). The two cases with the largest predicted CPM (>80 µm), larger CPS (>0.0119) and lowest bone modulus (<472 MPa) were confirmed to have the largest in vivo RSA 3D migration (>0.78 mm). This study enabled the first corroboration between tailored musculoskeletal and FE model predictions with in vivo RSA cup migration. Investigation of additional patient-matched CT, gait, and RSA examinations may allow further development and validation of FE models.

Examining agreement between finite element modelling methodologies in predicting pathological fracture risk in proximal femurs with bone metastases.

BACKGROUND: Finite element modelling methodologies available for assessing femurs with metastases accurately predict strength and pathological fracture risk which has led them to being considered for implementation into the clinic. However, the models available use varying material models, loading conditions, and critical thresholds. The aim of this study was to determine the agreement between finite element modelling methodologies in assessing fracture risk in proximal femurs with metastases. METHODS: CT images of the proximal femur were obtained of 7 patients who presented with a pathologic femoral fracture (fracture group) and the contralateral femur of 11 patients scheduled for prophylactic surgery (non-fracture group). Fracture risk was predicted for each patient following three established finite modelling methodologies which have previously shown to accurately predict strength and determine fracture risk: non-linear isotropic -based model, strain fold ratio -based model, Hoffman failure criteria -based model. FINDINGS: The methodologies demonstrated good diagnostic accuracy in assessing fracture risk (AUC = 0.77, 0.73, and 0.67). There was a stronger monotonic association between the non-linear isotropic and Hoffman -based models (τ = 0.74) than with the strain fold ratio model (τ = -0.24 and - 0.37). There was moderate or low agreement between methodologies in discriminating between individuals at high or low risk of fracture (κ = 0.20, 0.39, and 0.62). INTERPRETATION: The present results suggest there may be a lack of consistency in the management of pathological fractures in the proximal femur based on the finite element modelling methodologies.

Determining the relationship between tibiofemoral geometry and passive motion with partial least squares regression.

Tibiofemoral geometry influences knee passive motion and understanding their relationship can provide insight into knee function and mechanisms of injury. However, the complexity of the geometric constraints has made characterizing the relationship challenging. The aim of this study was to determine the tibiofemoral bone geometries that explain the variation in passive motion using a partial least squares regression (PLSR) model. The PLSR model was developed for 29 healthy cadaver specimens (10 female, 19 male) with femur and tibia geometries retrieved from MRI images and six degree-of-freedom tibiofemoral kinematics determined during a flexion cycle with minimal medial pressure. The first 13 partial least squares (PLS) components explained 90% of the variation in the kinematics and accounted for 89% of the variation in geometry. The first three PLS components which shared geometric changes to particular surface congruencies of the tibial and femoral condyles explained the most amount of variation in the kinematics, primarily in anterior-posterior translation. Meanwhile, variations in femoral condyle width and the intercondylar space, tibia plateau size and conformity, and tibia eminences heights in PLS 2 and 4 explained the greatest amount of variation in internal-external rotation. PLS 4 exhibiting variation in overall size of the knee accounted for greatest amount of variation in geometry (50%) and had the greatest influence on the abduction-adduction motion and some on internal-external rotation but, overall, explained only a small proportion of the kinematics (10%). Elucidating the complex relationship between tibiofemoral bone geometry and passive kinematics may help personalize treatments for improved functional outcomes in patients.

The role of muscle forces and gait cycle discretization when assessing acetabular cup primary stability: A finite element study.

The aim of this study was to investigate the influence of the muscle force contribution and loading cycle discretization on the predicted micromotion and interfacial bone strains in the implanted acetabulum. To this end, a patient specific finite element model of the hemipelvis was developed, based on the CT-scan and gait analysis results, collected as part of the authors' previous work. Outcomes of this study suggests that the acetabular cup micromotion and interfacial bone strains can be predicted just using the joint contact force. This helps to reduce the complexity of the finite element models by ignoring the contribution of muscle forces and the associated challenges of mapping these forces to the pelvis. However, the gait cycle needs to be adequately discretised to capture the micromotion at the bone-implant interface. BACKGROUND AND OBJECTIVE: The Dalstra load case, which includes muscle forces, has been widely adopted in the literature for studying the mechanical environment in the intact and implanted acetabulum. To simplify the modelling approach, some researchers ignore the contribution of muscle forces. The Dalstra load case is also divided into eight separate load steps (five in the stance phase and three in the swing phase), however, it is unclear whether this adequately captures the micromotions, for a cementless acetabular cup, during a simulated activity. The aim of this study was to investigate the influence of the muscle force contribution and loading cycle discretization on the predicted micromotion and interfacial bone strains. METHODS: In this work, a patient specific finite element model of the hemipelvis was developed, based on the CT-scan and gait analysis results, collected as part of the authors' previous work. Finite element simulations were performed using the joint contact and muscle forces derived from two sources. The first approach was used the load case proposed by Dalstra et al. The second approach used joint contact and muscle forces predicted by a musculoskeletal model. Additionally, the musculoskeletal load case was discretised into 50 equal load steps and the results compared with the equivalent Dalstra load steps. RESULTS: The results showed that the contribution of the muscle forces resulted in minor differences in both the magnitude and distribution of the predicted acetabular micromotion (up to 4.01% in the mean acetabular micromotion) and interfacial bone strains (up to 10.34% in the mean interfacial bone strains). The degree of gait cycle discretisation had a significant influence on the acetabular micromotion with a difference of 20.89% in the mean acetabular micromotion. CONCLUSION: Outcomes of this study suggests that the acetabular cup micromotion and interfacial bone strains can be predicted just using the joint contact force. This helps to reduce the complexity of the finite element models by ignoring the contribution of muscle forces and the associated challenges of mapping these forces to the pelvis. However, the gait cycle needs to be adequately discretised to capture the micromotion at the bone-implant interface.

Geometry and bone mineral density determinants of femoral neck strength changes following exercise.

Physical exercise induces spatially heterogeneous adaptation in bone. However, it remains unclear where the changes in BMD and geometry have the greatest impact on femoral neck strength. The aim of this study was to determine the principal BMD-and-geometry changes induced by exercise that have the greatest effect on femoral neck strength. Pre- and post-exercise 3D-DXA images of the proximal femur were collected of male participants from the LIFTMOR-M exercise intervention trial. Meshes with element-by-element correspondence were generated by morphing a template mesh to each bone to calculate changes in BMD and geometry. Finite element (FE) models predicted femoral neck strength changes under single-leg stance and sideways fall load. Partial least squares regression (PLSR) models were developed with BMD-only, geometry-only, and BMD-and-geometry changes to determine the principal modes that explained the greatest variation in neck strength changes. The PLSR models explained over 90% of the strength variation with 3 PLS components using BMD-only (R2 > 0.92, RMSE < 0.06 N) and 8 PLS components with geometry-only (R2 > 0.93, RMSE < 0.06 N). Changes in the superior neck and distal cortex were most important during single-leg stance while the superior neck, medial head, and lateral trochanter were most important during a sideways fall. Local changes in femoral neck and head geometry could differentiate the exercise groups from the control group. Exercise interventions may target BMD changes in the superior neck, inferior neck, and greater trochanter for improved femoral neck strength in single-leg stance and sideways fall.

Predicting rupture locations of ascending aortic aneurysms using CT-based finite element models.

Accurate rupture risk assessment of ascending aortic aneurysms is important for reducing aneurysm-related mortality. More recently, computational models have been shown to better predict rupture risk than diameter-based measurements. However, it remains unclear whether finite element (FE) models of the ascending aorta can predict rupture location, and over what timeframe those predictions are reliable. The aim of this study was to evaluate FE models of the ascending aorta generated from computed tomography (CT) scans in predicting rupture location. Pre- and post-rupture CT scans were obtained of 12 patients who underwent emergency surgical repair for ascending aorta rupture with varying time intervals between scans (20 days - 6 years). A rigid iterative closest point (ICP) registration was used to overlay post-rupture aortic geometries with pre-rupture FE models and directly compare predicted regions of high equivalent strain with actual rupture. The FE model predicted the rupture location in the 5 patients with the shortest time intervals between the pre- and post-rupture scans (20 days - 2 years, 3 months). However, rupture location was not predicted in the 4/5 patients with greater than 3 years between scans. Achieving a physiological equivalent strain distribution in the FE model was highly dependent on the resolution of the pre-rupture scan and whether contrast agent was present. The results suggest there may be a time interval beyond which FE predictions of rupture location may not be reliable. The results in this study may help clinical validation of FE models of ascending aortic aneurysms predicting rupture risk.

A semiautomated method to quantitatively assess osteolytic lesion volume and bone mineral density within acetabular regions of interest from CT.

The objectives of this study were to (1) develop a semiautomated method to obtain lesion volume and bone mineral density (BMD) in terms of Hounsfield units from pelvic computed tomography (CT) scans in three regions of interest, and (2) assess accuracy and reliability of the method based on cadaveric CT scans. Image artefacts due to metal implants reduce CT clarity and are more severe with more than one implant in situ. Therefore, accuracy and reliability tests were performed with varying numbers of total hip arthroplasties implanted. To test the accuracy of lesion size measurements, microcomputed tomography was used as a reference. Mean absolute error ranged from 36 to 284 mm3 after five measurements. Intra- and inter-operator reliability of the entire method was measured for a selection of parameters. All coefficient of variation values were good to excellent for CT scans of the native pelvic anatomy and a CT scans of the same pelvis with one and two implants in situ. Accuracy of quantifying lesion volume decreased with decreasing CT image clarity by 0.6%-3.6% mean absolute relative error. Reliability of lesion volume measurement decreased with decreasing CT clarity. This was also the case for reliability of BMD measurements in the region most disrupted by metal artefact. The presented method proposes an approach for quantifying bone loss which has been proven to be accurate, reliable, and clinically applicable.

Histological regional analysis of the aortic root and thoracic ascending aorta: a complete analysis of aneurysms from root to arch.

BACKGROUND: Although aortic root and ascending aortic aneurysms are treated the same, they differ in embryological development and pathological processes. This study examines the microscopic structural differences between aortic root and ascending aortic aneurysms, correlating these features to the macroscopic pathophysiological processes. METHODS: We obtained surgical samples from ascending aortic aneurysms (n = 11), aortic root aneurysms (n = 3), and non-aneurysmal patients (n = 7), Aortic collagen and elastin content were examined via histological analysis, and immunohistochemistry techniques used to determine collagen I, III, and IV subtypes. Analysis was via observational features, and colour deconvolution quantification techniques. RESULTS: Elastin fiber disruption and fragmentation was the most extensive in the proximal aneurysmal regions. Medial fibrosis and collagen density increased in proximal aneurysmal regions and aortic root aneurysms (p < 0.005). Collagen I was seen in highest quantity in aortic root aneurysms. Collagen I content was greatest in the sinus tissue regions compared to the valvular and ostial regions (p < 0.005) Collagen III and IV quantification did not vary greatly. The most susceptible regions to ultrastructural changes in disease are the proximal ascending aorta and aortic root. CONCLUSIONS: The aortic root differs histologically from the ascending aorta confirming its unique composition in aneurysm pathology. These findings should prompt further evaluation on the influence of this altered structure on function which could potentially guide clinical management.

Assessment of femoral neck strength and bone mineral density changes following exercise using 3D-DXA images.

Physical exercise induces spatially heterogeneous bone changes in the proximal femur. Recent advances have enabled 3D dual-energy X-ray Absorptiometry (DXA)-based finite element (FE) models to estimate bone strength. However, its ability to detect exercise-induced BMD and strength changes is unclear. The aim of this study was to quantify the repeatability of vBMD and femoral neck strength obtained from 3D-DXA images and determine the changes due an exercise intervention. The DXA scans included pairs of same-day repeated scans from ten healthy females and pre- and post-exercise intervention scans of 26 males. FE models with element-by-element correspondence were generated by morphing a template mesh to each bone. BMD and femoral strength under single-leg-stance and sideways fall loading configurations were obtained for both groups and compared. In the repeated images, the total hip vBMD difference was 0.5 ± 2.5%. Element-by-element BMD differences reached 30 ± 50%. The strength difference in single-leg stance was 2.8 ± 13% and in sideways fall was 4.5% ± 19%. In the exercise group, strength changes were 6 ± 19% under single-leg stance and 1 ± 18% under sideways fall. vBMD parameters were weakly correlated to strength (R2 < 0.31). The exercise group had a mean bone accrual exceeding repeatability values in the femoral head and cortical regions. The case with the highest vBMD change (6.4%) caused 18% and -7% strength changes under single-leg stance and sideways fall. 3D-DXA technology can assess the effect of exercise interventions in large cohorts but its validity in individual cases should be interpreted with caution.

The functional limits of the aneurysmal aortic root. A unique pressure testing apparatus.

BACKGROUND: The aortic root has unique embryological development and is a highly sophisticated and complex structure. In studies that report on the biomechanical characteristics of the thoracic aorta, distinction between the aortic root and ascending aorta regions is nonexistent. Our objective is to determine the maximal pressures at which dissection occurs or tissue failure occurs in the aortic root compared to that of the ascending aorta in the presence of aortic aneurysms. This may help guide preoperative monitoring, diagnosis and the decision for operative intervention for aortic root aneurysms in the normal and susceptible populations. METHODS: We developed a simple aortic root and ascending aorta pressure testing unit in series. Ten fresh porcine hearts were obtained from the local abattoir (n = 5 aortic root and n = 5 ascending aorta for comparison). Using a saline filled needle and syringe, artificial fluid-filled aneurysms were created between the intima and medial layers of the aortic root. The aorta lumen was then progressively filled with saline solution. Pressure measurement was taken at time of loss of tissue integrity, obvious tissue dissection or aneurysm rupture, and the tissue structure was then visually examined. RESULTS: In the aortic root, mean maximal pressure (mmHg) at tissue failure was 208 mmHg. Macroscopic examination revealed luminal tears around the coronary ostia in 2/5 specimens, and in all specimens, there was propagation of the dissection in the aortic root in a circumferential direction. In all ascending aorta specimens, the maximal aortic pressures exceeded 300 mmHg without tissue failure or dissection, and eventual apparatus failure. CONCLUSION: Our results indicate that the aneurysmal aortic root tissues are at greater risk of rupture and dissection propagation at lower aortic pressure. With further analysis, this could guide clinical and surgical management.

Patient and surgical variability in the primary stability of cementless acetabular cups: A finite element study.

Aseptic loosening is the most common indication for revision of cementless acetabular cups and often depends on the primary stability achieved following surgery. Cup designs must be capable of achieving primary stability for a wide variety of individuals and surgical conditions to be successful. Typically, preclinical finite element (FE) testing of cups involves assessing the performance in a single patient and under a limited set of idealized conditions. The aim of this study was to assess the effect of patient and surgical parameters on the primary stability of an acetabular cup design in a set of subject-specific FE models. Interference fit was varied in a representative set of 12 patient-specific models of the implanted hemipelvis. Linear mixed models showed a significant association with micromotion for interference fit (P < .0001), acetabular bone elastic modulus (P < .001), native acetabular diameter (P = .03), and the interference fit-elastic modulus interaction (P = .01). There were no significant associations between the polar gap and any of the parameters considered. The significant interference fit-elastic modulus interaction suggests that increasing the interference fit in patients with low bone quality leads to a greater reduction in micromotion than in patients with higher bone quality. However, the significant association between percentage bone yielding and interference fit (P < .0001) suggests a higher periacetabular fracture risk at higher interference fits. This work supports the development of preclinical testing of cup designs for the broad range patients and surgical conditions a cup may face following surgery.

Influence of stems and metaphyseal sleeve on primary stability of cementless revision tibial trays used to reconstruct AORI IIB defects.

Metaphyseal augments, such as sleeves, have been introduced to augment the fixation of revision total knee replacement (rTKR) components, and can be used with or without a stem. The effect of sleeve size in combination with stems on the primary stability and load transfer of a rTKR implant in AORI type IIB defects where the defect involves both condyles are poorly understood. The aim of this study was to examine the primary stability of revision tibial tray augmented with a sleeve in an AORI type IIB defect which involves both condyles with loss of cortical and cancellous bone. Finite element models were generated from computed tomography (CT) scans of nine individuals. All the bones used in the study had an AORI type IIB defect. The cohort included eight females (mean weight: 64 kg, height: 1.6 m). Material properties were sampled from CT data and assigned to the FE model. Joint contact forces for level gait, stair descent, and squat were applied. Stemless sleeved implants under various loading conditions were shown to have adequate primary stability in all AORI type IIB defects investigated. Adding a stem only marginally improved the primary stability of the implant but reduced the strain in the metaphysis compared to stemless implants. Once good initial mechanical stability was established with a sleeve, there was no benefit, in terms of primary stability or bone strains, from increasing sleeve size. This study suggests that metaphyseal sleeves, without a stem, can provide the required primary stability required by a rTKR tibial implant, to reconstruct an AORI type IIB defect. © 2019 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res.

Virtual trial to evaluate the robustness of cementless femoral stems to patient and surgical variation.

Primary stability is essential for the success of cementless femoral stems. In this study, patient specific finite element (FE) models were used to assess changes in primary stability due to variability in patient anatomy, bone properties and stem alignment for two commonly used cementless femoral stems, Corail® and Summit® (DePuy Synthes, Warsaw, USA). Computed-tomography images of the femur were obtained for 8 males and 8 females. An automated algorithm was used to determine the stem position and size which minimized the endo-cortical space, and then span the plausible surgical envelope of implant positions constrained by the endo-cortical boundary. A total of 1952 models were generated and ran, each with a unique alignment scenario. Peak hip contact and muscle forces for stair climbing were scaled to the donor's body weight and applied to the model. The primary stability was assessed by comparing the implant micromotion and peri-prosthetic strains to thresholds (150 µm and 7000 µÎµ, respectively) above which fibrous tissue differentiation and bone damage are expected to prevail. Despite the wide range of implant positions included, FE prediction were mostly below the thresholds (medians: Corail®: 20-74 µm and 1150-2884 µÎµ, Summit®: 25-111 µm and 860-3010 µÎµ), but sensitivity of micromotion and interfacial strains varied across femora, with the majority being sensitive (p < 0.0029) to average bone mineral density, cranio-caudal angle, post-implantation anteversion angle and lateral offset of the femur. The results confirm the relationship between implant position and primary stability was highly dependent on the patient and the stem design used.

Sampling strategies for approximating patient variability in population-based finite element studies of total hip replacement.

Total hip replacements must be robust to patient variability for long-term success in the population. The challenge during the design process is evaluating an implant in a diverse population but the computational cost of simulating a population of subject-specific finite element (FE) models is not practical. We examined five strategies to generate representative subsets of subjects from a cohort of 103 implanted hip joint FE models to approximate the variability in output metrics. Comparing with the median and distribution of the 95th percentile composite peak micromotion (CPM) and polar gap in the full cohort (CPM median: 136 µm, interquartile range [IQR]: 74-230 µm) (Polar Gap median: 467 µm, IQR: 434-548 µm), the Anatomic Sampling strategy (12 subjects) achieved the best balance of computational cost and approximation of the output metrics (CPM median: 169 µm, IQR: 78-236 µm) (Polar Gap median: 469 µm, IQR: 448-537 µm). Convex hull sampling (41 subjects) more closely captured the output metrics (CPM median: 99 µm, IQR: 70-191 µm) (Polar Gap median: 456 µm, IQR: 418-533 µm) but required over three times the number of subjects. Volume reduction of the convex hull captured the extremes of variability with subsets of 5 to 20 subjects, while the largest minimum-distance strategy captured the variability toward the middle of the cohort. These strategies can estimate the level of variability in FE model output metrics with a low computational cost when evaluating implants during the design process.

Biomechanical Robustness of a Contemporary Cementless Stem to Surgical Variation in Stem Size and Position.

Successful designs of total hip replacement (THR) need to be robust to surgical variation in sizing and positioning of the femoral stem. This study presents an automated method for comprehensive evaluation of the potential impact of surgical variability in sizing and positioning on the primary stability of a contemporary cementless femoral stem (Corail®, DePuy Synthes). A patient-specific finite element (FE) model of a femur was generated from computed tomography (CT) images from a female donor. An automated algorithm was developed to span the plausible surgical envelope of implant positions constrained by the inner cortical boundary. The analysis was performed on four stem sizes: oversized, ideal (nominal) sized, and undersized by up to two stem sizes. For each size, Latin hypercube sampling was used to generate models for 100 unique alignment scenarios. For each scenario, peak hip contact and muscle forces published for stair climbing were scaled to the donor's body weight and applied to the model. The risk of implant loosening was assessed by comparing the bone-implant micromotion/strains to thresholds (150 µm and 7000 µÎµ) above which fibrous tissue is expected to prevail and the periprosthetic bone to yield, respectively. The risk of long-term loosening due to adverse bone resorption was assessed using bone adaptation theory. The range of implant positions generated effectively spanned the available intracortical space. The Corail stem was found stable and robust to changes in size and position, with the majority of the bone-implant interface undergoing micromotion and interfacial strains that are well below 150 µm and 7000 µÎµ, respectively. Nevertheless, the range of implant positions generated caused an increase of up to 50% in peak micromotion and up to 25% in interfacial strains, particularly for retroverted stems placed in a medial position.

Influence of varying stem and metaphyseal sleeve size on the primary stability of cementless revision tibial trays used to reconstruct AORI IIA defects. A simulation study.

Traditionally, diaphyseal stems have been utilized to augment the stability of revision total knee replacement (rTKR) implants. More recently metaphyseal augments, such as sleeves, have been introduced to further augment component fixation. The effect of augments such as stems and sleeves have on the primary stability of a rTKR implant is poorly understood, however it has important implications on the complexity, costs and survivorship of the procedure. Finite element analysis was used to investigate the primary stability and strain distribution of various size stems and sleeves used in conjunction with a cementless revision tibial tray. The model was built from computer tomography images of a single healthy tibia obtained from an 81-year-old patient to which an Anderson Orthopaedic Research Institute (AORI) IIA defect was virtually added. The influences of varying body mass index (BMI) and bone modulus were also investigated. Stemless sleeves were found to provided adequate primary implant stability (average implant micro-motion <50 µm) for the studied defect. Addition of a stem did not enhance the primary stability. Furthermore, this study found that varying BMI and bone modulus had a considerable effect on strain distribution but negligible effect on micro-motion in the sleeve area. In conclusion, the addition of diaphyseal stem to a metaphyseal sleeve had little benefit in enhancing the primary stability of tibial trays augmented when simulating reconstructions of AORI IIA tibial defects. Additional studies are required to determine the relative benefit of the diaphyseal stem when using metaphyseal sleeves defects with more extensive bone loss. © 2018 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:1876-1886, 2018.

Primary stability of a cementless acetabular cup in a cohort of patient-specific finite element models.

The primary stability achieved during total hip arthroplasty determines the long-term success of cementless acetabular cups. Pre-clinical finite element testing of cups typically use a model of a single patient and assume the results can be extrapolated to the general population. This study explored the variability in predicted primary stability of a Pinnacle® cementless acetabular cup in 103 patient-specific finite element models of the hemipelvis and examined the association between patient-related factors and the observed variability. Cups were inserted by displacement-control into the FE models and then a loading configuration simulating a complete level gait cycle was applied. The cohort showed a range of polar gap of 284-1112 µm and 95th percentile composite peak micromotion (CPM) of 18-624 µm. Regression analysis was not conclusive on the relationship between patient-related factors and primary stability. No relationship was found between polar gap and micromotion. However, when the patient-related factors were categorised into quartile groups, trends suggested higher polar gaps occurred in subjects with small and shallow acetabular geometries and cup motion during gait was affected most by low elastic modulus and high bodyweight. The variation in primary stability in the cohort for an acetabular cup with a proven clinical track record may provide benchmark data when evaluating new cup designs. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:1012-1023, 2018.

A Computational Efficient Method to Assess the Sensitivity of Finite-Element Models: An Illustration With the Hemipelvis.

Assessing the sensitivity of a finite-element (FE) model to uncertainties in geometric parameters and material properties is a fundamental step in understanding the reliability of model predictions. However, the computational cost of individual simulations and the large number of required models limits comprehensive quantification of model sensitivity. To quickly assess the sensitivity of an FE model, we built linear and Kriging surrogate models of an FE model of the intact hemipelvis. The percentage of the total sum of squares (%TSS) was used to determine the most influential input parameters and their possible interactions on the median, 95th percentile and maximum equivalent strains. We assessed the surrogate models by comparing their predictions to those of a full factorial design of FE simulations. The Kriging surrogate model accurately predicted all output metrics based on a training set of 30 analyses (R2 = 0.99). There was good agreement between the Kriging surrogate model and the full factorial design in determining the most influential input parameters and interactions. For the median, 95th percentile and maximum equivalent strain, the bone geometry (60%, 52%, and 76%, respectively) was the most influential input parameter. The interactions between bone geometry and cancellous bone modulus (13%) and bone geometry and cortical bone thickness (7%) were also influential terms on the output metrics. This study demonstrates a method with a low time and computational cost to quantify the sensitivity of an FE model. It can be applied to FE models in computational orthopaedic biomechanics in order to understand the reliability of predictions.