Tissue differentiation around a short stemmed metaphyseal loading implant employing a modified mechanoregulatory algorithm: a finite element study

Bone ingrowth in randomly distributed porous interbody cage during lumbar spinal fusion

Porous interbody cages are often used in spinal fusion surgery since they allow bone ingrowth which facilitates long-term stability. However, the extent of bone ingrowth in and around porous interbody cages has scarcely been investigated. Moreover, tissue differentiation might not be similar around the superior and inferior cage-bone interfaces. Using mechanobiology-based numerical framework and physiologic loading conditions, the study investigates the spatial distribution of evolutionary bone ingrowth within randomly distributed porous interbody cages, having varied porosities. Finite Element (FE) microscale models, corresponding to cage porosities of 60 %, 72 %, and 83 %, were developed for the superior and inferior interfacial regions of the cage, along with the macroscale model of the implanted lumbar spine. The implant-bone relative displacements of different porosity models were mapped from macroscale to microscale model. Bone formation of 10-40 % was predicted across the porous cage models, resulting in an average Young's modulus ranging between 765 MPa and 915 MPa. Maximum bone ingrowth of ∼34 % was observed for the 83 % porous cage, which was subject to low implant-bone relative displacements (maximum 50μm). New bone formation was found to be greater at the superior interface (∼34 %) as compared to the inferior interface (∼30 %) for P83 model. Relatively greater volume of fibrous tissue was formed at the implant-bone interface for the cage with 60 % and 72 % porosities, which might lead to cage migration and eventual failure of the implant. Hence, the interbody cage with 83 % porosity appears to be most favorable for bone ingrowth, provided sufficient mechanical strength is offered.

Prediction of bone ingrowth into a porous novel hip-stem: A finite element analysis integrated with mechanoregulatory algorithm

Bone ingrowth into a porous implant is necessary for its long-term fixation. Although attempts have been made to quantify the peri-implant bone growth using finite element (FE) analysis integrated with mechanoregulatory algorithms, bone ingrowth into a porous cellular hip stem has scarcely been investigated. Using a three-dimensional (3D) FE model and mechanobiology-based numerical framework, the objective of this study was to predict the spatial distribution of evolutionary bone ingrowth into an uncemented novel porous hip stem proposed earlier by the authors. A CT-based FE macromodel of the implant-bone structure was developed. The bone material properties were assigned based on CT grey value. Peak musculoskeletal loading conditions, corresponding to level walking and stair climbing, were applied. The geometry of the implant-bone macromodel was divided into multiple submodels. A suitable mapping framework was used to transfer maximum nodal displacements from the FE macromodel to the cut boundaries of the FE submodels. CT grey value-based bone materials properties were assigned to the submodels. Thereafter, the submodels were solved and simulations of bone ingrowth were carried out using mechanoregulatory principle. A gradual increase in the average Young's modulus, from 1200 to 1500 MPa, of the bone tissue layer was observed considering all the submodels. The distal submodel exhibited 82% of bone ingrowth, whereas the proximal submodel experienced 65% bone ingrowth. Equilibrium in the bone ingrowth process was achieved in 7 weeks postoperatively, with a notable amount of bone ingrowth that should lead to biological fixation of the novel hip stem.

Assessment of bone ingrowth around beaded coated tibial implant for total ankle replacement using mechanoregulatory algorithm

The long-term performance of porous coated tibial implants for total ankle replacement (TAR) primarily depends on the extent of bone ingrowth at the bone-implant interface. Although attempts were made for primary fixation for immediate post-operative stability, no investigation was conducted on secondary fixation. The aim of this study is to assess bone ingrowth around the porous beaded coated tibial implant for TAR using a mechanoregulatory algorithm. A realistic macroscale finite element (FE) model of the implanted tibia was developed based on computer tomography (CT) data to assess implant-bone micromotions and coupled with microscale FE models of the implant-bone interface to predict bone ingrowth around tibial implant for TAR. The macroscale FE model was subjected to three near physiological loading conditions to evaluate the site-specific implant-bone micromotion, which were then incorporated into the corresponding microscale model to mimic the near physiological loading conditions. Results of the study demonstrated that the implant experienced tangential micromotion ranged from 0 to 71 μm with a mean of 3.871 μm. Tissue differentiation results revealed that bone ingrowth across the implant ranged from 44 to 96 %, with a mean of around 70 %. The average Young's modulus of the inter-bead tissue layer varied from 1444 to 4180 MPa around the different regions of the implant. The analysis postulates that when peak micromotion touches 30 μm around different regions of the implant, it leads to pronounced fibrous tissues on the implant surface. The highest amount of bone ingrowth was observed in the central regions, and poor bone ingrowth was seen in the anterior parts of the implant, which indicate improper osseointegration around this region. This macro-micro mechanical FE framework can be extended to improve the implant design to enhance the bone ingrowth and in future to develop porous lattice-structured implants to predict and enhance osseointegration around the implant.

Computational assessment of growth of connective tissues around textured hip stem subjected to daily activities after THA

Longer-term stability of uncemented femoral stem depends on ossification at bone-implant interface. Although attempts have been made to assess the amount of bone growth using finite element (FE) analysis in combination with a mechanoregulatory algorithm, there has been little research on tissue differentiation patterns on hip stems with proximal macro-textures. The primary goal of this investigation is to qualitatively compare the formation of connective tissues around a femoral implant with/without macro-textures on its proximal surfaces. This study also predicts formation of different tissue phenotypes and their spatio-temporal distribution around a macro-textured femoral stem under routine activities. Results from the study show that non-textured implants (80 to 94%) encourage fibroplasia compared to that in textured implants (71 to 85.38%) under similar routine activity, which might trigger aseptic loosening of implant. Formation of bone was more on medio-lateral sides and towards proximal regions of Gruen zones 2 and 6, which was found to be in line with clinical observations. Fibroplasia was higher under stair climbing (85 to 91%) compared to that under normal walking (71 to 85.38%). This study suggests that stair climbing, although falls under recommended activity, might be detrimental to patient compared to normal walking in the initial rehabilitation period.

Application of finite element analysis to tissue differentiation and bone remodelling approaches and their use in design optimization of orthopaedic implants: A review

Post-operative bone growth and long-term bone adaptation around the orthopaedic implants are simulated using the mechanoregulation based tissue-differentiation and adaptive bone remodelling algorithms, respectively. The primary objective of these algorithms was to assess biomechanical feasibility and reliability of orthopaedic implants. This article aims to offer a comprehensive review of the developments in mathematical models of tissue-differentiation and bone adaptation and their applications in studies involving design optimization of orthopaedic implants over three decades. Despite the different mechanoregulatory models developed, existing literature confirm that none of the models can be highly regarded or completely disregarded over each other. Not much development in mathematical formulations has been observed from the current state of knowledge due to the lack of in vivo studies involving clinically relevant animal models, which further retarded the development of such models to use in translational research at a fast pace. Future investigations involving artificial intelligence (AI), soft-computing techniques and combined tissue-differentiation and bone-adaptation studies involving animal subjects for model verification are needed to formulate more sophisticated mathematical models to enhance the accuracy of pre-clinical testing of orthopaedic implants.

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.

The influence of macro-textural designs over implant surface on bone on-growth: A computational mechanobiology based study

The longerterm secondary stability of an uncemented implant depends primarily on the quality and extent of bone in-growth or on-growth at the bone-implant interface. Investigations are warranted to predict the influences of implant macro-textures on bone on-growth pattern. Mechanoregulatory tissue differentiation algorithms can predict such patterns effectively. There is, however, a dearth of volumetric in silico study to assess the influence of macro-textures on bone growth. The present study investigated the influence of macro-textural grooves/ribs on changes in tissue formation at the bone-implant interface by carrying out a 3D finite element (FE) analysis. Three distinct macro-textures, loosely based on commercially viable hip stem models, were comparatively assessed for varying levels of interfacial micromotion. The study predicted elevated fibrogenesis and chondrogenesis, followed by a suppressed osteogenesis for higher levels of micromotion (60 μm and 100 μm), resulting in weak bone-implant interface strength. However, small judicious modifications in implant surface texture may enhance bone growth to a considerable extent. The numerical scheme can further be used as a template for more rigorous parametric and multi-scale studies.

A comprehensive analysis of bio-inspired design of femoral stem on primary and secondary stabilities using mechanoregulatory algorithm

The coated porous section of stem surface is initially filled with callus that undergoes osseointegration process, which develops a bond between stem and bone, lessens the micromotions and transfers stresses to the bone, proximally. This phenomenon attributes to primary and secondary stabilities of the stems that exhibit trade-off the stem stiffness. This study attempts to ascertain the influence of stem stiffness on peri-prosthetic bone formation and stress shielding when in silico models of solid CoCr alloy and Ti alloy stems, and porous Ti stems (53.8 GPa and 31.5 GPa Young's moduli) were implanted. A tissue differentiation predictive mechanoregulation algorithm was employed to estimate the evolutionary bond between bone and stem interfaces with 0.5-mm- and 1-mm-thick calluses. The results revealed that the high stiffness stems yielded higher stress shielding and lower micromotions than that of low stiffness stems. Contrarily, bone formation around solid Ti alloy stem and porous Ti 53.8 GPa stem was augmented in 0.5-mm- and 1-mm-thick calluses, respectively. All designs of stems exhibited different rates of bone formation, diverse initial micromotions and stress shielding; however, long-term bone formation was coherent with different stress shielding. Therefore, contemplating the secondary stability of the stems, low stiffness stem (Ti 53.8 GPa) gave superior biomechanical performance than that of high stiffness stems.

Up to 10-year follow-up of the Symax stem in THA: a Danish single-centre study

INTRODUCTION

The design of the cementless Symax-HA femoral stem is based on geometrical analysis of human femoral anatomy to optimise the fit within the femur. The stem combines an anatomical proximal section enabling a metaphyseal anchorage with a straight distal section. This results in an improved proximal load transfer. We aimed to determine the overall survival of total hip arthroplasty (THA) in patients operated at our institution with the Symax stem. Secondly, we investigated the cause of revision for patients having THA including the Symax stem.

METHODS

The study is a registry study with results from a single centre, Vejle Hospital. From the Danish Hip -Arthroplasty Registry, we identified all THAs operated with the Symax stem. The primary outcome was revision. The secondary outcomes were aseptic loosening, periprosthetic fracture, and all other causes for revision.

RESULTS

In total, 1,055 THAs were performed in 875 patients. The median 6.5-year survival rate was 97.5% (95% confidence interval [CI]: 96.6%-98.3%). 29 THAs had undergone revision due to the following causes: periprosthetic fracture (n = 11); aseptic loosening of the cup (n = 1); infection (n = 5); dislocation (n = 10); and component failure (n = 2). There was no revision due to aseptic loosening of the stem.

CONCLUSIONS

From the Danish Hip Arthroplasty Registry, we identified the overall 6.5-year survival rate to be 95% for cementless THAs. In the present study, the Symax stem has a median 6.5-year survival rate at 97.5% (96.6%-98.3%). No revisions were due to aseptic loosening but the stem had a relative high prevalence of periprosthetic fractures.

Influence of functionally graded pores on bone ingrowth in cementless hip prosthesis: a finite element study using mechano-regulatory algorithm

Cementless hip prostheses with porous outer coating are commonly used to repair the proximally damaged femurs. It has been demonstrated that stability of prosthesis is also highly dependent on the bone ingrowth into the porous texture. Bone ingrowth is influenced by the mechanical environment produced in the callus. In this study, bone ingrowth into the porous structure was predicted by using a mechano-regulatory model. Homogenously distributed pores (200 and 800 [Formula: see text]m in diameter) and functionally graded pores along the length of the prosthesis were introduced as a porous coating. Bone ingrowth was simulated using 25 and 12 [Formula: see text]m micromovements. Load control simulations were carried out instead of traditionally used displacement control. Spatial and temporal distributions of tissues were predicted in all cases. Functionally graded pore decreasing models gave the most homogenous bone distribution, the highest bone ingrowth (98%) with highest average Young's modulus of all tissue phenotypes approximately 4.1 GPa. Besides this, the volume of the initial callus increased to 8.33% in functionally graded pores as compared to the 200 [Formula: see text]m pore size models which increased the bone volume. These findings indicate that functionally graded porous surface promote bone ingrowth efficiently which can be considered to design of surface texture of hip prosthesis.

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.

Influence of Implant Surface Texture Design on Peri-Acetabular Bone Ingrowth: A Mechanobiology Based Finite Element Analysis

The fixation of uncemented acetabular components largely depends on the amount of bone ingrowth, which is influenced by the design of the implant surface texture. The objective of this numerical study is to evaluate the effect of these implant texture design factors on bone ingrowth around an acetabular component. The novelty of this study lies in comparative finite element (FE) analysis of 3D microscale models of the implant-bone interface, considering patient-specific mechanical environment, host bone material property and implant-bone relative displacement, in combination with sequential mechanoregulatory algorithm and design of experiment (DOE) based statistical framework. Results indicated that the bone ingrowth process was inhibited due to an increase in interbead spacing from 200 μm to 600 μm and bead diameter from 1000 μm to 1500 μm and a reduction in bead height from 900 μm to 600 μm. Bead height, a main effect, was found to have a predominant influence on bone ingrowth. Among the interaction effects, the combination of bead height and bead diameter was found to have a pronounced influence on bone ingrowth process. A combination of low interbead spacing (P = 200 μm), low bead diameter (D = 1000 μm), and high bead height (H = 900 μm) facilitated peri-acetabular bone ingrowth and an increase in average Young's modulus of newly formed tissue layer. Hence, such a surface texture design seemed to provide improved fixation of the acetabular component.

Bone ingrowth around porous-coated acetabular implant: a three-dimensional finite element study using mechanoregulatory algorithm

Fixation of uncemented implant is influenced by peri-prosthetic bone ingrowth, which is dependent on the mechanical environment of the implant-bone structure. The objective of the study is to gain an insight into the tissue differentiation around an acetabular component. A mapping framework has been developed to simulate appropriate mechanical environment in the three-dimensional microscale model, implement the mechanoregulatory tissue differentiation algorithm and subsequently assess spatial distribution of bone ingrowth around an acetabular component, quantitatively. The FE model of implanted pelvis subjected to eight static load cases during a normal walking cycle was first solved. Thereafter, a mapping algorithm has been employed to include the variations in implant-bone relative displacement and host bone material properties from the macroscale FE model of implanted pelvis to the microscale FE model of the beaded implant-bone interface. The evolutionary tissue differentiation was observed in each of the 13 microscale models corresponding to 13 acetabular regions. The total implant-bone relative displacements, averaged over each region of the acetabulum, were found to vary between 10 and 60 μm. Both the linear elastic and biphasic poroelastic models predicted similar mechanoregulatory peri-prosthetic tissue differentiation. Considerable variations in bone ingrowth (13-88%), interdigitation depth (0.2-0.82 mm) and average tissue Young's modulus (970-3430 MPa) were predicted around the acetabular cup. A progressive increase in the average Young's modulus, interdigitation depth and decrease in average radial strains of newly formed tissue layer were also observed. This scheme can be extended to investigate tissue differentiation for different surface texture designs on the implants.

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.

Bone remodeling after total hip arthroplasty with a short stemmed metaphyseal loading implant: finite element analysis validated by a prospective DEXA investigation

In total hip arthroplasty (THA), short stemmed cementless implants are used because they are thought to stimulate physiological bone remodeling and reduce stress shielding. We performed a numerical investigation on bone remodeling after implantation of a specific short stemmed implant using finite element analysis (FEA). Overall bone mass loss was 2.8% in the entire femur. Bone mass decrease was mostly found in the proximal part of the calcar and in the greater trochanter due to the vast cross section of the implant, probably leading to stress shielding. In the diaphysis, no change in the apparent bone density was proven. The assumptions made agreed well with bone remodeling data from THA recipients who underwent dual-energy X-ray absorptiometry. However, the clinical investigation revealed a bone mass increase in the minor trochanter region that was less pronounced in the FEA. Further comparisons to other stem designs must be done to verify if the relative advantages of the investigated implant can be accepted.