Accounting for patient variability in finite element analysis of the intact and implanted hip and knee: a review

Stem to prevent periprosthetic fracture after notching in total knee arthroplasty

Improper osteotomy during total knee arthroplasty (TKA) can lead to anterior femoral notching, which increases the risk of periprosthetic fractures due to stress concentration. One potential solution is the addition of an intramedullary stem to the femoral component. However, the optimal stem length remains unclear. In this study, we aimed to determine the optimal stem length using finite element models. Finite element models of femurs were developed with unstemmed prostheses and prostheses with stem lengths of 50, 75, and 100 mm. Under squat loading conditions, the von Mises stress at the notch and stress distribution on four transversal sections of the femur were analyzed. Additionally, micromotion of the prosthesis-bone interface was evaluated to assess initial stability. The unstemmed prosthesis exhibited a von Mises stress of 191.8 MPa at the notch, which decreased to 43.1, 8.8, and 23.5 MPa for stem lengths of 50, 75, and 100 mm, respectively. The stress reduction on four selected femoral transversal sections compared with the unstemmed prosthesis was 40.0%, 84.4%, and 67.1% for stem lengths of 50, 75, and 100 mm, respectively. Micromotion analysis showed a maximum of 118.8 μm for the unstemmed prosthesis, which decreased significantly with the application of stems, particularly at the anterior flange. Intramedullary stems effectively reduced stress concentration at the femoral notch. The 50-mm stem length provided the optimal combination of reduced notch stress, minimized stress-shielding effect, and decreased micromotion at the anterior flange.

Accuracy tradeoffs between individual bone and joint-level statistical shape models of knee morphology

Statistical shape models (SSMs) are useful tools in evaluating variation in bony anatomy to assess pathology, plan surgical interventions, and inform the design of orthopaedic implants and instrumentation. Recently, by considering multiple bones spanning a joint or the whole lower extremity, SSMs can support studies investigating articular conformity and joint mechanics. The objective of this study was to assess tradeoffs in accuracy between SSMs of the femur or tibia individually versus a combined joint-level model. Three statistical shape models were developed (femur-only, tibia-only, and joint-level) for a training set of 179 total knee arthroplasty (TKA) patients with osteoarthritis representing both genders and several ethnicities. Bone geometries were segmented from preoperative CT scans, meshed with triangular elements, and registered to a template for each SSM. Principal component analysis was performed to determine modes of variation. The statistical shape models were compared using measures of compactness, accuracy, generalization, and specificity. The generalization evaluation, assessing the ability to describe an unseen instance in a leave-one-out analysis, showed that errors were consistently smaller for the individual femur and tibia SSMs than for the joint-level model. However, when additional modes were included in the joint-level model, the errors were comparable to the individual bone results, with minimal additional computational expense. When developing more complex SSMs at the joint, lower limb, or whole-body level, the use of an error threshold to inform the number of included modes, instead of 95 % of the variation explained, can help to ensure accurate representations of anatomy.

Using a finite element model of the thumb to study Trapeziometacarpal joint contact during lateral pinch

BACKGROUND

Finite element (FE) analysis is widely used in different fields of orthopaedic surgery, however, its application to the trapeziometacarpal joint has been limited due to the small size, complex biconcave-convex joint geometry, and complex musculature. The goal of this study was to improve upon existing models by creating a muscle-driven FE thumb model and use the model to simulate the biomechanical effect of hand therapy exercises and ligament reconstructive surgeries.

METHODS

Bone and cartilage geometry were based on a CT dataset of a subject performing a static lateral pinch task. A previously validated musculoskeletal model was utilized to extract electromyography (EMG)-driven muscle forces. Five ligaments with biomechanical significance were modeled as springs using literature values and attached according to their anatomical landmarks.

FINDINGS

The biomechanical consequence of various interventions was proxied as a change in the maximum cartilage stress. The result shows tightening the dorsal ligament complex (dorsal radial ligament, dorsal central ligament, posterior oblique ligament) is the most effective, achieving a stress reduction of 4.8%. Five exercises used in hand therapies were modeled, among which thenar eminence strengthening showed the most prominent stress reduction of 4.0%. Four ligament reconstructive surgeries were modeled, with Eaton-Littler reconstruction showed the most significant stress reduction of 25.0%.

INTERPRETATION

Among the routinely utilized treatment options for early thumb osteoarthritis, we found that three methods: dorsal ligament imbrication, thenar eminence exercise, and the Eaton-Littler method may confer biomechanical advantages cartilage loading. These advantages align with the clinically observed favorable outcomes.

LibHip: An open-access hip joint model repository suitable for finite element method simulation

BACKGROUND AND OBJECTIVE

population-based finite element analysis of hip joints allows us to understand the effect of inter-subject variability on simulation results. Developing large subject-specific population models is challenging and requires extensive manual effort. Thus, the anatomical representations are often subjected to simplification. The discretized geometries do not guarantee conformity in shared interfaces, leading to complications in setting up simulations. Additionally, these models are not openly accessible, challenging reproducibility. Our work provides multiple subject-specific hip joint finite element models and a novel semi-automated modeling workflow.

METHODS

we reconstruct 11 healthy subject-specific models, including the sacrum, the paired pelvic bones, the paired proximal femurs, the paired hip joints, the paired sacroiliac joints, and the pubic symphysis. The bones are derived from CT scans, and the cartilages are generated from the bone geometries. We generate the whole complex's volume mesh with conforming interfaces. Our models are evaluated using both mesh quality metrics and simulation experiments.

RESULTS

the geometry of all the models are inspected by our clinical expert and show high-quality discretization with accurate geometries. The simulations produce smooth stress patterns, and the variance among the subjects highlights the effect of inter-subject variability and asymmetry in the predicted results.

CONCLUSIONS

our work is one of the largest model repositories with respect to the number of subjects and regions of interest in the hip joint area. Our detailed research data, including the clinical images, the segmentation label maps, the finite element models, and software tools, are openly accessible on GitHub and the link is provided in Moshfeghifar et al.(2022)[1]. Our aim is to empower clinical researchers to have free access to verified and reproducible models. In future work, we aim to add additional structures to our models.

A Method to Compare Heterogeneous Types of Bone and Cartilage Meshes

Accurately capturing the bone and cartilage morphology and generating a mesh remains a critical step in the workflow of computational knee joint modeling. Currently, there is no standardized method to compare meshes of different element types and nodal densities, making comparisons across research teams a significant challenge. The aim of this paper is to describe a method to quantify differences in knee joint bone and cartilages meshes, independent of bone and cartilage mesh topology. Bone mesh-to-mesh distances, subchondral bone boundaries, and cartilage thicknesses from meshes of any type of mesh are obtained using a series of steps involving registration, resampling, and radial basis function fitting after which the comparisons are performed. Subchondral bone boundaries and cartilage thicknesses are calculated and visualized in a common frame of reference for comparison. The established method is applied to models developed by five modeling teams. Our approach to obtain bone mesh-to-mesh distances decreased the divergence seen in selecting a reference mesh (i.e., comparing mesh A-to-B versus mesh B-to-A). In general, the bone morphology was similar across teams. The cartilage thicknesses for all models were calculated and the mean absolute cartilage thickness difference was presented, the articulating areas had the best agreement across teams. The teams showed disagreement on the subchondral bone boundaries. The method presented in this paper allows for objective comparisons of bone and cartilage geometry that is agnostic to mesh type and nodal density.

Application of statistical shape modeling to the human hip joint: a scoping review protocol

OBJECTIVE

This review aims to identify all examples of the application of statistical shape models to the human hip joint, with a focus on methodology, validation, and applications.

INTRODUCTION

Abnormal hip joint morphology (eg, deformity secondary to Legg-Calvé-Perthes disease) is an important precursor to osteoarthritis. Clinical radiographs are often used to characterize deformity and provide indication for surgical correction, but it is unclear whether radiographs can adequately describe three-dimensional deformity. Statistical shape modeling, a method of describing a population of shapes using a small number of variables, has been identified as a potential tool that will allow clinicians and researchers to validate current and novel radiographic measurements of hip deformity. In identifying all previous examples of statistical shape modeling applied to the hip joint, this review will determine its prevalence, strengths, and weaknesses, and identify gaps in the literature.

INCLUSION CRITERIA

Peer-reviewed and gray literature focusing on the development and/or application of statistical shape models to the human hip joint will be included.

METHODS

Several relevant databases, including Ovid MEDLINE, Embase, and IEEE, will be searched for literature published from 1992, and for a title and abstract that can be searched in English. After removal of duplicates, two reviewers will independently screen papers by title and abstract, then screen the full text of selected or uncertain papers. The same reviewers will then independently chart data from the final selection. At each stage, disagreements will be resolved through discussion or third-party arbitration.

Computational framework for population-based evaluation of TKR-implanted patellofemoral joint mechanics

Differences in patient anatomy are known to influence joint mechanics. Accordingly, intersubject anatomical variation is an important consideration when assessing the design of joint replacement implants. The objective of this study was to develop a computational workflow to perform population-based evaluations of total knee replacement implant mechanics considering variation in patient anatomy and to assess the potential for an efficient sampling strategy to support design phase screening analyses. The approach generated virtual subject anatomies using a statistical shape model of the knee and performed virtual implantation to size and align the implants. A finite-element analysis simulated a deep knee bend activity and predicted patellofemoral (PF) mechanics. The study predicted bounds of performance for kinematics and contact mechanics and investigated relationships between patient factors and outputs. For example, the patella was less flexed throughout the deep knee bend activity for patients with an alta patellar alignment. The results also showed the PF range of motions in AP and ML were generally larger with increasing femoral component size. Comparison of the 10-90% bounds between sampling strategies agreed reasonably, suggesting that Latin Hypercube sampling can be used for initial screening evaluations and followed up by more intensive Monte Carlo simulation for refined designs. The platform demonstrated a functional workflow to consider variation in joint anatomy to support robust implant design.

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.

Computational efficient method for assessing the influence of surgical variability on primary stability of a contemporary femoral stem in a cohort of subjects

Finite element (FE) modelling can provide detailed information on implant stability; however, its computational cost limits the possibility of completing large numerical analyses into the effect of surgical variability in a cohort of patients. The aim of this study was to develop an efficient surrogate model for a cohort of patients implanted using a common cementless hip stem. FE models of implanted femora were generated from computed tomography images for 20 femora (11 males, 9 females; 50-80 years; 52-116 kg). An automated pipeline generated FE models for 61 different unique scenarios that span the femur-specific range of implant positions. Peak hip contact and muscle forces for stair climbing were scaled to the donors' body weight and applied to the models. A cohort-specific surrogate for implant micromotion was constructed from Gaussian process models trained using data from FE simulations representing the median and extreme implant positions for each femur. A convergence study confirmed suitability of the sampling method for cohorts with 10+ femora. The final model was trained using data from the 20 femora. Results showed very good agreement between the FE and the surrogate predictions for a total of 1036 alignment scenarios [root mean squared error (RMSE) < 20 µm; [Formula: see text] = 0.81]. The total time required for the surrogate model to predict the micromotion range associated with surgical variability was approximately one-eighth of the corresponding full FE analysis. This confirms that the developed model is an accurate yet computationally cheaper alternative to full FE analysis when studying the implant robustness in a cohort of 10+ femora.

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 95^th 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.

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.

Investigating gender and ethnicity differences in proximal humeral morphology using a statistical shape model

Morphological variability in the shoulder influences the joint biomechanics and is an important consideration in arthroplasty and implant design. The objectives of this study were to quantify cortical and cancellous proximal humeral morphology and to assess whether shape variation was influenced by gender and ethnicity, with the overarching goal of informing implant design and treatment. A statistical shape model of the proximal humeral cortical and cancellous regions was developed for a training set of 84 subjects of both genders and different ethnicities. Cortical and cancellous bone geometries were reconstructed from CT scans, meshed with triangular elements, and registered to a template. Principal component analysis was applied to quantify modes of variation. Anatomical measurements were computed on the registered geometries to assess correlation with modes of variation. Parallel analysis identified six significant modes of variation, which accounted for 93% of variation in the training set and described scaling (Mode 1), inclination of the head (Modes 2 and 5), and shape of the greater tuberosity and neck region (Modes 3, 4, and 6). Size differences as described by Mode 1 were statistically significant for gender and ethnicity, where female and Asian subjects were smaller than male and Caucasian subjects, respectively; however, differences in other modes were not significant. Cortical thickness of the shaft after normalization by outer diameter was significantly larger for Asian subjects compared to Caucasian subjects. The statistical shape model quantified cortical and cancellous humeral morphology considering gender and ethnicity, providing descriptive data to support surgical planning, and implant design. © 2018 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:3043-3052, 2018.

Evaluating the primary stability of standard vs lateralised cementless femoral stems - A finite element study using a diverse patient cohort

BACKGROUND

Restoring the original femoral offset is desirable for total hip replacements as it preserves the original muscle lever arm and soft tissue tensions. This can be achieved through lateralised stems, however, the effect of variation in the hip centre offset on the primary stability remains unclear.

METHODS

Finite element analysis was used to compare the primary stability of lateralised and standard designs for a cementless femoral stem (Corail®) across a representative cohort of male and female femora (N = 31 femora; age from 50 to 80 years old). Each femur model was implanted with three designs of the Corail® stem, each designed to achieve a different degree of lateralisation. An automated algorithm was used to select the size and position that achieve maximum metaphyseal fit for each of the designs. Joint contact and muscle forces simulating the peak forces during level gait and stair climbing were scaled to the body mass of each subject.

FINDINGS

The study found that differences in restoring the native femoral offset introduce marginal differences in micromotion (differences in peak micromotion <21 μm), for most cases. Nonetheless, significant reduction in the interfacial strains (>3000 με) was achieved for some subjects when lateralized stems were used.

INTERPRETATION

Findings of this study suggest that, with the appropriate size and alignment, the standard offset design is likely to be sufficient for primary stability, in most cases. Nonetheless, appropriate use of lateralised stems has the potential reduce the risk of peri-prosthetic bone damage. This highlights the importance of appropriate implant selection during the surgical planning stage.

How can a short stem hip implant preserve the natural, pre-surgery force flow? A finite element analysis on a collar cortex compression concept (CO4)

The present work proposes a simple, novel fixation concept for short stem hip endoprostheses, which preserves the pre-surgery force flow through femoral bone to an unprecedented extent. It is demonstrated by finite element analyses that a standard implant model endowed with minor geometrical changes can overcome bone loading reduction and can achieve almost physiological conditions. The numerical results underpin that the key aspect of the novel, so-called "collar cortex compression concept CO4" is the direct, almost full load transmission from the implant collar to the resected femur cortex, which implies that the implant stem must be smooth and therefore interacts mainly by normal contact with the surrounding bone. For a stem endowed with surface porosity at already small areas, it is mainly the stem which transmits axial forces by shear, whereas the collar shows considerable unloading, which is the standard metaphyseal fixation. Only in the latter case the implant-bone stiffness contrast induces stress shielding, whereas for CO4 stress shielding is avoided almost completely, although the implant is made of a stiff Ti-alloy. CO4 is bionics-inspired in that it mimics force transmission at implant-bone interfaces following the natural conditions and it thereby preserves pre-surgery bone architecture as an optimized solution of nature.

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 collars on the primary stability of cementless femoral stems: A finite element study using a diverse patient cohort

For cementless femoral stems, there is debate as to whether a collar enhances primary stability and load transfer compared to collarless designs. Finite Element (FE) analysis has the potential to compare stem designs within the same cohort, allowing for subtle performance differences to be identified, if present. Subject-specific FE models of intact and implanted femora were run for a diverse cohort (21 males, 20 females; BMI 16.4-41.2 kg/m² , age 50-80 yrs). Collared and collarless versions of Corail^® (DePuy Synthes, Warsaw, IN) were sized and positioned using an automated algorithm that aligns the femoral/stem axes, preserves the head-center location, and maximizes metaphyseal fit. Joint contact and muscle forces simulating peak forces in level gait and stair climbing and were scaled to the body mass and applied to each subject. Three failure scenarios were assessed: Potential for peri-prosthetic fibrous tissue formation (stem micromotion), potential for peri-prosthetic bone damage (equivalent strains), and calcar bone remodeling (changes in strain-energy density). Comparisons were performed using paired t-tests. Only subtle differences were found (mean 90th percentile micromotion: Collared = 86 µm, collarless = 92.5 µm, mean 90th percentile interface strains: Collared = 733 µϵ, collarless = 767 µϵ, and similar remodeling stimuli were predicted). The slight differences observed were small in comparison with the inter-patient variability. Statement of clinical significance: Our results suggest that the presence/absence of a collar is unlikely to substantially alter the bone-implant biomechanics nor the initial mechanical environment. Hence, a collar is likely to have minimal clinical impact. Analysis using different femoral stem designs is recommended before generalising these findings. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:1185-1195, 2018.

Assessment of fracture risk in proximal tibia with tumorous bone defects by a finite element method

There has not been a satisfying method to predict the fracture risk in tumorous bone lesions. To tackle this challenge, we used a finite element method to assess the fracture risk in the proximal tibia (pT) when the size and location of the tumorous defects are varied in bone. Towards this end, the circular cortical defects, mimicking the tumorous lesions by forming cortical window defects, with a diameter (Ф) of 20, 30, 40, or 50 mm, are structured on the anteromedial, lateral, posterior wall of pT, which is located 5, 15, and 25 mm below articular margin, respectively. We found that under walking conditions, the Von Mises Stress of each defective tibia model was larger than that of the intact tibia model and also showed a positive linear correlation with the sizes of the defects. A notable fracture risk was not observed until the defect was Ф30 mm or larger. When the defect emerged, the anteromedial wall resisted fracture risk more than the rest of wall. Our results show that the size and location of the bone tumors are important factors affecting the fracture risk of pT. Our findings will be beneficial to clinicians when deciding what treatment to use for pT lesions.

Multiobjective optimization of cartilage stress for non-invasive, patient-specific recommendations of high tibial osteotomy correction angle - a novel method to investigate alignment correction

Medial opening wedge high tibial osteotomy (MOWHTO) is a surgical procedure to treat knee osteoarthritis associated with varus deformity. However, the ideal final alignment of the Hip-Knee-Ankle (HKA) angle in the frontal plane, that maximizes procedural success and post-operative knee function, remains controversial. Therefore, the purpose of this study was to introduce a subject-specific modeling procedure in determining the biomechanical effects of MOWHTO alignment on tibiofemoral cartilage stress distribution. A 3D finite element knee model derived from magnetic resonance imaging of a healthy participant was manipulated in-silico to simulate a range of final HKA angles (i.e. 0.2°, 2.7°, 3.9° and 6.6° valgus). Loading and boundary conditions were assigned based on subject-specific kinematic and kinetic data from gait analysis. Multiobjective optimization was used to identify the final alignment that balanced compressive and shear forces between medial and lateral knee compartments. Peak stresses decreased in the medial and increased in the lateral compartment as the HKA was shifted into valgus, with balanced loading occurring at angles of 4.3° and 2.9° valgus for the femoral and tibial cartilage respectively. The concept introduced here provides a platform for non-invasive, patient-specific preoperative planning of the osteotomy for medial compartment knee osteoarthritis.

The effects of musculoskeletal loading regimes on numerical evaluations of acetabular component

The importance of clinical studies notwithstanding, the failure assessment of implant–bone structure has alternatively been carried out using finite element analysis. However, the accuracy of the finite element predicted results is dependent on the applied loading and boundary conditions. Nevertheless, most finite element–based evaluations on acetabular component used a few selective load cases instead of the eight load cases representing the entire gait cycle. These in silico evaluations often suffer from limitations regarding the use of simplified musculoskeletal loading regimes. This study attempts to analyse the influence of three different loading regimes representing a gait cycle, on numerical evaluations of acetabular component. Patient-specific computer tomography scan-based models of intact and resurfaced pelvises were used. One such loading regime consisted of the second load case that corresponded to peak hip joint reaction force. Whereas the other loading regime consisted of the second and fifth load cases, which corresponded to peak hip joint reaction force and peak muscle forces, respectively. The third loading regime included all the eight load cases. Considerable deviations in peri-acetabular strains, standard error ranging between 115 and 400 µε, were observed for different loading regimes. The predicted bone strains were lower when selective loading regimes were used. Despite minor quantitative variations in bone density changes (less than 0.15 g cm−3), the final bone density pattern after bone remodelling was found to be similar for all the loading regimes. Underestimations in implant–bone micromotions (40–50 µm) were observed for selective loading regimes after bone remodelling. However, at immediate post-operative condition, such underestimations were found to be less (less than 5 µm). The predicted results highlight the importance of inclusion of eight load cases representing the gait cycle for in silico evaluations of resurfaced pelvis.

Comparison of proximal humeral bone stresses between stemless, short stem, and standard stem length: a finite element analysis

BACKGROUND

The stem lengths of humeral components used in shoulder arthroplasty vary; however, the literature on these devices is limited. This finite element study investigates the effect of humeral component length on stresses in the proximal humerus.

METHODS

Intact and 3 reconstructed (standard length, short, and stemless implants) finite element models were created from shoulder computed tomography scan data (N = 5). Loading was simulated at varying abduction angles (15°, 45°, and 75°). The average bone stress (represented as a percentage of intact values) was reported at 8 transverse slices. In addition, the overall average change in cortical and trabecular bone stresses was quantified.

RESULTS

Cortical bone stresses in the most proximal slice for the standard (58% ± 12%) and short (78% ± 10%) stem models were significantly reduced compared with the intact (100%) and stemless (101% ± 6%) models (P = .005). These reductions persisted in the second cortical slice for the standard stem compared with the intact, stemless, and short models (P = .025). Interestingly, stresses in the trabecular bone within these proximal slices were significantly elevated when stemless implants were used compared with all other implants (P < .001), regardless of abduction angle.

CONCLUSION

Reducing stem length produced humeral stresses that more closely matched the intact stress distribution in proximal cortical bone. Opposing trends presented in the proximal trabecular bone, probably because of differences in load transfer when shorter stems are used. Accordingly, the results suggest that implant stem length is 1 variable that can be modified in an attempt to better mimic intact bone stresses during humeral component insertion, provided stem fixation is adequate.

Comparison of proximal femur and vertebral body strength improvements in the FREEDOM trial using an alternative finite element methodology

Denosumab reduced the incidence of new fractures in postmenopausal women with osteoporosis by 68% at the spine and 40% at the hip over 36 months compared with placebo in the FREEDOM study. This efficacy was supported by improvements from baseline in vertebral (18.2%) strength in axial compression and femoral (8.6%) strength in sideways fall configuration at 36 months, estimated in Newtons by an established voxel-based finite element (FE) methodology. Since FE analyses rely on the choice of meshes, material properties, and boundary conditions, the aim of this study was to independently confirm and compare the effects of denosumab on vertebral and femoral strength during the FREEDOM trial using an alternative smooth FE methodology. Unlike the previous FE study, effects on femoral strength in physiological stance configuration were also examined. QCT data for the proximal femur and two lumbar vertebrae were analyzed by smooth FE methodology at baseline, 12, 24, and 36 months for 51 treated (denosumab) and 47 control (placebo) subjects. QCT images were segmented and converted into smooth FE models to compute bone strength. L1 and L2 vertebral bodies were virtually loaded in axial compression and the proximal femora in both fall and stance configurations. Denosumab increased vertebral body strength by 10.8%, 14.0%, and 17.4% from baseline at 12, 24, and 36 months, respectively (p<0.0001). Denosumab also increased femoral strength in the fall configuration by 4.3%, 5.1%, and 7.2% from baseline at 12, 24, and 36 months, respectively (p<0.0001). Similar improvements were observed in the stance configuration with increases of 4.2%, 5.2%, and 5.2% from baseline (p≤0.0007). Differences between the increasing strengths with denosumab and the decreasing strengths with placebo were significant starting at 12 months (vertebral and femoral fall) or 24 months (femoral stance). Using an alternative smooth FE methodology, we confirmed the significant improvements in vertebral body and proximal femur strength previously observed with denosumab. Estimated increases in strength with denosumab and decreases with placebo were highly consistent between both FE techniques.

A Hertzian Integrated Contact Model of the Total Knee Replacement Implant for the Estimation of Joint Contact Forces

The prediction of lower limb muscle and contact forces may provide useful knowledge to assist the clinicians in the diagnosis as well as in the development of appropriate treatment for musculoskeletal disorders. Research studies have commonly estimated joint contact forces using model-based muscle force estimation due to the lack of a reliable contact model and material properties. The objective of this present study was to develop a Hertzian integrated contact model. Then, in vivo elastic properties of the Total Knee Replacement (TKR) implant were identified using in vivo contact forces leading to providing reliable material properties for modeling purposes. First, a patient specific rigid musculoskeletal model was built. Second, a STL-based implant model was designed to compute the contact area evolutions during gait motions. Finally, a Hertzian integrated contact model was defined for the in vivo identification of elastic properties (Young’s modulus and Poisson coefficient) of the instrumented TKR implant. Our study showed a potential use of a new approach to predict the contact forces without knowledge of muscle forces. Thus, the outcomes may lead to accurate and reliable prediction of human joint contact forces for new case study.

Four decades of finite element analysis of orthopaedic devices: where are we now and what are the opportunities?

Finite element has been used for more than four decades to study and evaluate the mechanical behaviour total joint replacements. In Huiskes seminal paper "Failed innovation in total hip replacement: diagnosis and proposals for a cure", finite element modelling was one of the potential cures to avoid poorly performing designs reaching the market place. The size and sophistication of models has increased significantly since that paper and a range of techniques are available from predicting the initial mechanical environment through to advanced adaptive simulations including bone adaptation, tissue differentiation, damage accumulation and wear. However, are we any closer to FE becoming an effective screening tool for new devices? This review contains a critical analysis of currently available finite element modelling techniques including (i) development of the basic model, the application of appropriate material properties, loading and boundary conditions, (ii) describing the initial mechanical environment of the bone-implant system, (iii) capturing the time dependent behaviour in adaptive simulations, (iv) the design and implementation of computer based experiments and (v) determining suitable performance metrics. The development of the underlying tools and techniques appears to have plateaued and further advances appear to be limited either by a lack of data to populate the models or the need to better understand the fundamentals of the mechanical and biological processes. There has been progress in the design of computer based experiments. Historically, FE has been used in a similar way to in vitro tests, by running only a limited set of analyses, typically of a single bone segment or joint under idealised conditions. The power of finite element is the ability to run multiple simulations and explore the performance of a device under a variety of conditions. There has been increasing usage of design of experiments, probabilistic techniques and more recently population based modelling to account for patient and surgical variability. In order to have effective screening methods, we need to continue to develop these approaches to examine the behaviour and performance of total joint replacements and benchmark them for devices with known clinical performance. Finite element will increasingly be used in the design, development and pre-clinical testing of total joint replacements. However, simulations must include holistic, closely corroborated, multi-domain analyses which account for real world variability.

Development of mapped stress-field boundary conditions based on a Hill-type muscle model

Forces generated in the muscles and tendons actuate the movement of the skeleton. Accurate estimation and application of these musculotendon forces in a continuum model is not a trivial matter. Frequently, musculotendon attachments are approximated as point forces; however, accurate estimation of local mechanics requires a more realistic application of musculotendon forces. This paper describes the development of mapped Hill-type muscle models as boundary conditions for a finite volume model of the hip joint, where the calculated muscle fibres map continuously between attachment sites. The applied muscle forces are calculated using active Hill-type models, where input electromyography signals are determined from gait analysis. Realistic muscle attachment sites are determined directly from tomography images. The mapped muscle boundary conditions, implemented in a finite volume structural OpenFOAM (ESI-OpenCFD, Bracknell, UK) solver, are employed to simulate the mid-stance phase of gait using a patient-specific natural hip joint, and a comparison is performed with the standard point load muscle approach. It is concluded that physiological joint loading is not accurately represented by simplistic muscle point loading conditions; however, when contact pressures are of sole interest, simplifying assumptions with regard to muscular forces may be valid.

A physically motivated constitutive model for 3D numerical simulation of skeletal muscles

A detailed numerical implementation within the FEM is presented for a physically motivated three-dimensional constitutive model describing the passive and active mechanical behaviors of the skeletal muscle. The derivations for the Cauchy stress tensor and the consistent material tangent are provided. For nearly incompressible skeletal muscle tissue, the strain energy function may be represented either by a coupling or a decoupling of the distortional and volumetric material response. In the present paper, both functionally different formulations are introduced allowing for a direct comparison between the coupled and decoupled isochoric-volumetric approach. The numerical validation of both implementations revealed significant limitations for the decoupled approach. For an extensive characterization of the model response to different muscle contraction modes, a benchmark model is introduced. Finally, the proposed implementation is shown to provide a reliable tool for the analysis of complex and highly nonlinear problems through the example of the human mastication system by studying bite force and three-dimensional muscle shape changes during mastication.

Finite element analysis on longitudinal and radial functionally graded femoral prosthesis

This study focused on developing a 3D finite element model of functionally graded femoral prostheses to decrease stress shielding and to improve total hip replacement performance. The mechanical properties of the modeled functionally graded femoral prostheses were adjusted in the sagittal and transverse planes by changing the volume fraction gradient exponent. Prostheses with material changes in the sagittal and transverse planes were considered longitudinal and radial prostheses, respectively. The effects of cemented and noncemented implantation methods were also considered in this study. Strain energy and von Mises stresses were determined at the femoral proximal metaphysis and interfaces of the implanted femur components, respectively. Results demonstrated that the strain energy increased proportionally with increasing volume fraction gradient exponent, whereas the interface stresses decreased on the prostheses surfaces. A limited increase was also observed at the surfaces of the bone and cement. The periprosthetic femur with a noncemented prosthesis exhibited higher strain energy than with a cemented prosthesis. Radial prostheses implantation displayed more strain energy than longitudinal prostheses implantation in the femoral proximal part. Functionally graded materials also increased strain energy and exhibited promising potentials as substitutes of conventional materials to decrease stress shielding and to enhance total hip replacement lifespan.