The effect of primary stability on load transfer and bone remodelling within the uncemented resurfaced femur

Custom design and biomechanical clinical trials of 3D-printed polyether ether ketone femoral shaft prosthesis

During the surgical resection and reconstruction of a pathological femoral fracture, the removal of the femoral tumor leaves a large bone defect. Thus, it is necessary to reconstruct the defect and perform internal fixation. Polyether ether ketone (PEEK) has been widely used in spinal fusion and cranioplasty given its excellent biomechanical properties, biocompatibility, and stability. The typical design method of femoral prosthesis is based on the contralateral mirror image model (M-model), and we propose a novel method for designing femoral prosthesis, which is based on the cross section and centerline of the mirrored femur (C-model). In this study, the femoral shaft prostheses based on two models were manufactured using fused deposition modeling technology, and we use mechanical test and finite element analysis (FEA) to reveal the differences in mechanical properties of the two models. The mechanical results showed that the maximum loading force and yield strength were increased by 3% and 6% in the C-model prosthesis compared with the M-model prosthesis, respectively. In FEA, the results indicate that the C-model prosthesis could reduce the stress concentration by 5.4%-10.9% compared to the M-model prosthesis. Finally, the 3D-printed PEEK femoral shaft prosthesis based on C-model was implanted, no early complications occurred. Postoperative radiological examination indicated that the prosthesis and the femoral osteotomy end were closely matched and fixed well.

Design optimization of functionally graded lattice infill total hip arthroplasty stem for stress shielding reduction

Reducing stress shielding of stem-inserted femurs in total hip arthroplasty caused by the high stiffness of the stem is an emerging medical engineering issue. In this study, a numerical design optimization methodology lattice infill stem was developed to realize a stem, balancing the low stiffness and strength requirements. Two pairs of models and loading conditions were introduced for the stress shielding and strength criteria. The objective function was set as the weighted sum of the criteria. Its effective density distribution was optimized by handling the representative size of the lattice as a design variable, assuming that the so-called body-centered cubic lattice was the base shape of the lattice. In the optimization, the approximated model of the lattice was handled as a solid material with the effective physical properties of the lattice derived by the homogenization method. After optimization, the detailed lattice stem geometry was modeled based on the obtained optimal lattice distribution, and the actual performance was numerically evaluated. The developed stem increased the stress applied to the remaining femur by 32.4% compared with the conventional stem.

Bone Ingrowth Around an Uncemented Femoral Implant Using Mechanoregulatory Algorithm: A Multiscale Finite Element Analysis

The primary fixation and long-term stability of a cementless femoral implant depend on bone ingrowth within the porous coating. Although attempts were made to quantify the peri-implant bone ingrowth using the finite element (FE) analysis and mechanoregulatory principles, the tissue differentiation patterns on a porous-coated hip stem have scarcely been investigated. The objective of this study is to predict the spatial distribution of evolutionary bone ingrowth around an uncemented hip stem, using a three-dimensional (3D) multiscale mechanobiology-based numerical framework. Multiple load cases representing a variety of daily living activities, including walking, stair climbing, sitting down, and standing up from a chair, were used as applied loading conditions. The study accounted for the local variations in host bone material properties and implant-bone relative displacements of the macroscale implanted FE model, in order to predict bone ingrowth in microscale representative volume elements (RVEs) of 12 interfacial regions. In majority RVEs, 20-70% bone tissue (immature and mature) was predicted after 2 months, contributing toward a progressive increase in average Young's modulus (1200-3000 MPa) of the interbead tissue layer. Higher bone ingrowth (mostly greater than 60%) was predicted in the anterolateral regions of the implant, as compared to the posteromedial side (20-50%). New bone tissue was formed deeper inside the interbead spacing, adhering to the implant surface. The study helps to gain an insight into the degree of osseointegration of a porous-coated femoral implant.

Qualitative predictions of bone growth over optimally designed macro-textured implant surfaces obtained using NN-GA based machine learning framework

The surface features on implant surface can improve biologic fixation of the implant with the host bone leading to improved secondary (biological) implant stability. Application of finite element (FE) based mechanoregulatory schemes to estimate the amount of bone growth for a wide range of implant surface features is either manually intensive or computationally expensive. This study adopts an integrated approach combining FE, back-propagation neural network (BPNN) and genetic algorithm (GA) based search to evaluate optimum surface macro-textures from three representative implant models so as to enhance bone growth. Initial surface textures chosen for the implant models were based on an earlier investigation. Based on FE predicted dataset, a BPNN was formulated for faster prediction of bone growth. Using the BPNN predicted output, a GA-based search was carried out to maximize bone growth subject to clinically admissible micromotion at the bone-implant interface. The results from FE analysis and bone growth predictions from the BPNN were found to have strong correlation. The optimal osseointegration-maximized-textures (OMTs) obtained were found to offer enhanced biological fixation, as compared to that offered by the textures in the initial models. Results from the present study reveal that certain reduction in the dimension of ribs/grooves promotes bone growth. However, periodic patterns of ribs with higher and lower rib dimensions provide uniform stress environment at the interface thus promoting osseointegration.

Numerical simulations on periprosthetic bone remodeling: a systematic review

BACKGROUND AND OBJECTIVE

The aim of the present study was to review the literature concerning the analysis of periprosthetic bone remodeling through finite element (FE) simulation.

METHODS

A systematic review was conducted on 9 databases, taking into account a ten-year time period (from 2009 until 2020). The inclusion criteria were: articles published in English, publication date after 2009, full text articles, articles containing the keywords both in the abstract and in the title. The articles were classified through the following parameters: dimensionality of the simulation, modelling of the bone-prosthesis interface, output parameters, type of simulated prosthesis, bone remodeling algorithm.

RESULTS

Sixty-seven articles were included in the study. Femur and tooth were the most evaluated bone segment (respectively 41.8% and 29.9%). The 55.2% of the evaluated articles used a bonded bone-prosthesis interface, 73% used 3D simulations, 67.2% of the articles (45 articles) evaluate the bone remodeling by the bone density variation. At last, 59.7% of the articles employed algorithms based on a specific remodeling function.

CONCLUSIONS

Increasing interest in the bone remodeling FE analysis in different bone segments emerged from the review, and heterogeneous solutions were adopted. An optimal balance between computational cost and accuracy is needed to accurately simulate the bone remodeling phenomenon in the post-operative period.

A comparative assessment of two designs of hip stem using rule-based simulation of combined osseointegration and remodelling

The stem–bone interface of cementless total hip arthroplasty undergoes an adaptive process of bone ingrowth until the two parts become osseointegrated. Another important phenomenon associated with aseptic loosening of hip stem is stress-shielding induced adverse bone remodelling. The objective of this study was to preclinically assess the relative performances of two distinct designs of hip stems by addressing the combined effect of bone remodelling and osseointegration, based on certain rule-based criteria obtained from the literature. Premised upon non-linear finite element analyses of patient-specific implanted femur models, the study attempts to ascertain in silico outcome of the hip stem designs based on an evolutionary interfacial condition, and to further comment on the efficacy of the rule-based technique on the prediction of peri-prosthetic osseointegration. One of the two hip stem models was a trade-off design obtained from an earlier shape optimization study, and the other was based on TriLock stem (DePuy). Both designs predicted similar long-term osseointegration (∼89% surface), although trade-off stem predicted higher post-operative osseointegration. Proximal bone resorption was found higher for TriLock (by ∼110%) as compared to trade-off model. The rule-based technique predicted clinically coherent osseointegration around both stems and appears to be an alternative to expensive mechanobiology-based schemes.

Bone remodelling of the humerus after a resurfacing and a stemless shoulder arthroplasty

BACKGROUND

New implant designs, such as resurfacing and stemless implants, have been developed to improve the long-term outcomes of the shoulder arthroplasty. However, it is not yet fully understood if their influence on the bone load distribution can compromise the long-term stability of the implant due to bone mass changes. Using three-dimensional finite element models, the aim of the present study was to analyse the bone remodelling process of the humerus after the introduction of resurfacing and stemless implants based on the Global C.A.P. and Sidus Stem-Free designs, respectively.

METHODS

The 3D geometric model of the humerus was generated from the CT data of the Visible Human Project and the resurfacing and stemless implants were modelled in Solidworks. Considering a native humerus model, a humerus model with the resurfacing implant, and a humerus model with the stemless implant, three finite element models were developed in Abaqus. Bone remodelling simulations were performed considering healthy and poor bone quality conditions. The loading condition considered comprised 6 load cases of standard shoulder movements, including muscle and joint reaction forces estimated by a multibody model of the upper limb.

FINDINGS

The results showed similar levels of bone resorption for the resurfacing and stemless implants for common humeral regions. The regions underneath the head of the resurfacing implant, unique to this design, showed the largest bone loss. For both implants, bone resorption was more pronounced for the poor bone quality condition than for the healthy bone quality condition.

INTERPRETATION

The stemless implant lost less density at the fixation site, which might suggest that these implants may be better supported in the long-term than the resurfacing implants. However, further investigation is necessary to allow definite recommendations.

THE EFFECT OF TRANSIENT LOADING ON A FOOT-ORTHOTIC USING TEMPORAL PARAMETERS OF GAIT

This study describes a method for performing transient finite element analysis (FEA) of an assistive device using experimental parameters obtained from gait analysis. A subject displaying pathologic gait, owing to lower limb deformity, was chosen for gait study. Using CAD tools, a remedial orthotic device was designed, which is expected to improve the gait of the subject. The orthotic model was subjected to static and transient loading conditions obtained from gait study, using an FEA tool. The stress ‘hot’ zones between the two modes of analysis are studied. In addition, the experimental gait data of a healthy control group were recorded to perform univariate regression studies for predicting the peak values of the normal forces, and validated by comparing with those available in the literature. The values thus obtained may be used for static behavioral analysis of assistive devices. From the FEA results, it can be conclusively said that the orthotic model is capable of sustaining gait cycle loading. The regression studies suggest the possibility of using anthropometric data to predict gait forces and subsequently perform static and transient loading analysis of assistive devices.

Combined Bone Ingrowth and Remodeling Around Uncemented Acetabular Component: A Multiscale Mechanobiology-Based Finite Element Analysis

Bone ingrowth and remodeling are two different evolutionary processes which might occur simultaneously. Both these processes are influenced by local mechanical stimulus. However, a combined study on bone ingrowth and remodeling has rarely been performed. This study is aimed at understanding the relationship between bone ingrowth and adaptation and their combined influence on fixation of the acetabular component. Based on three-dimensional (3D) macroscale finite element (FE) model of implanted pelvis and microscale FE model of implant–bone interface, a multiscale framework has been developed. The numerical prediction of peri-acetabular bone adaptation was based on a strain-energy density-based formulation. Bone ingrowth in the microscale models was simulated using the mechanoregulatory algorithm. An increase in bone strains near the acetabular rim was observed in the implanted pelvis model, whereas the central part of the acetabulum was observed to be stress shielded. Consequently, progressive bone apposition near the acetabular rim and resorption near the central region were observed. Bone remodeling caused a gradual increase in the implant–bone relative displacements. Evolutionary bone ingrowth was observed around the entire acetabular component. Poor bone ingrowth of 3–5% was predicted around the centro-inferio and inferio-posterio-superio-peripheral regions owing to higher implant–bone relative displacements, whereas the anterio-inferior and centro-superior regions exhibited improved bone ingrowth of 35–55% due to moderate implant–bone relative displacement. For an uncemented acetabular CoCrMo component, bone ingrowth had hardly any effect on bone remodeling; however, bone remodeling had considerable influence on bone ingrowth.

EVALUATION AND PREDICTION OF HUMAN GAIT PARAMETERS USING UNIVARIATE, MULTIVARIATE AND STEPWISE STATISTICAL METHODS

This research was carried out to establish the relationship between human anthropometric data and corresponding gait variables. A group comprising 35 participants (18 male and 17 female) was selected for the current study. The study consisted of trials in which each participant was asked to walk the length of the instrumented walkway (Kistler’s force platform inset) at a self-selected speed. Using a four-camera motion analysis system, the kinematic and kinetic parameters of each trial were calculated. The peak values obtained from the data curves were used to generate the necessary regression fits. In order to establish the correlation between the anthropometric data of human and the gait parameters, the univariate, multivariate and stepwise fits were generated. Further, the statistical methods were employed to evaluate the [Formula: see text], [Formula: see text] and [Formula: see text]-values for each fit. The current multivariate study indicates an increasing trend in [Formula: see text] values and decreasing trend for [Formula: see text]-values when compared with the univariate fits and the results follow the expected line.

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.

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.

Stress-shielding induced bone remodeling in cementless shoulder resurfacing arthroplasty: a finite element analysis and in vivo results

Cementless surface replacement arthroplasty (CSRA) of the shoulder was designed to preserve the individual anatomy and humeral bone stock. A matter of concern in resurfacing implants remains the stress shielding and bone remodeling processes. The bone remodeling processes of two different CSRA fixation designs, conical-crown (Epoca RH) and central-stem (Copeland), were studied by three-dimensional (3-D) finite element analysis (FEA) as well as evaluation of contact radiographs from human CSRA retrievals. FEA included one native humerus model with a normal and one with a reduced bone stock quality. Compressive strains were evaluated before and after virtual CSRA implantation and the results were then compared to the bone remodeling and stress-shielding pattern of eight human CSRA retrievals (Epoca RH n=4 and Copeland n=4). FEA revealed for both bone stock models increased compressive strains at the stem and outer implant rim for both CSRA designs indicating an increased bone formation at those locations. Unloading of the bone was seen for both designs under the central implant shell (conical-crown 50-85%, central-stem 31-93%) indicating high bone resorption. Those effects appeared more pronounced for the reduced than for the normal bone stock model. The assumptions of the FEA were confirmed in the CSRA retrieval analysis which showed bone apposition at the outer implant rim and stems with highly reduced bone stock below the central implant shell. Overall, clear signs of stress shielding were observed for both CSRAs designs in the in vitro FEA and human retrieval analysis. Especially in the central part of both implant designs the bone stock was highly resorbed. The impact of these bone remodeling processes on the clinical outcome as well as long-term stability requires further evaluation.