Multi-factorial analysis of variables influencing the bone loss of an implant placed in the maxilla: Prediction using FEA and SED bone remodeling algorithm

The aim of this study was to investigate the interactions of implant position, implant–abutment connection and loading condition influencing bone loss of an implant placed in the maxilla using finite element (FE) analysis and mathematical bone remodeling theory. The maxilla section contours were acquired using CT images to construct FE models containing RS (internal retaining-screw) and the TIS (taper integrated screwed-in) implants placed in SC (along the axis of occlusal force) and RA (along the axis of residual ridge) positions. The adaptive strain energy density (SED) algorithm was combined with FE approach to study the preliminary bone remodeling around implant systems under different load conditions. The simulated results showed that the implant position obviously influenced the bone loss. An implant placed in the RA position resulted in substantially increased bone loss. Implant receiving a lateral load slightly increased bone loss compared with an axial load. The implant type did not significantly influence bone loss. It was found that buccal site suffered the most bone loss around the implant, followed by distal, lingual and mesial sites. The implant position primarily influenced bone loss and it was found most obviously at the buccal site. Implant placed along the axial load direction of a proposed prosthesis could obtain less bone loss around the implant. Attaining proper occlusal adjustments to reduce the lateral occlusal force is recommended in implant–bone–prosthesis system. Abutments of internal engagement with or without taper-fit did not affect the bone loss in the surrounding bone.     

Comparison of different hip prosthesis shapes considering micro-level bone remodeling and stress-shielding criteria using three-dimensional design space topology optimization

Since the late 1980s, computational analysis of total hip arthroplasty (THA) prosthesis components has been completed using macro-level bone remodeling algorithms. The utilization of macro-sized elements requires apparent bone densities to predict cancellous bone strength, thereby, preventing visualization and analysis of realistic trabecular architecture. In this study, we utilized a recently developed structural optimization algorithm, design space optimization (DSO), to perform a micro-level three-dimensional finite element bone remodeling simulation on the human proximal femur pre- and post-THA. The computational simulation facilitated direct performance comparison between two commercially available prosthetic implant stems from Zimmer Inc.: the Alloclassic and the Mayo conservative. The novel micro-level approach allowed the unique ability to visualize the trabecular bone adaption post-operation and to quantify the changes in bone mineral content by region. Stress-shielding and strain energy distribution were also quantified for the immediate post-operation and the stably fixated, post-remodeling conditions. Stress-shielding was highest in the proximal region and remained unchanged post-remodeling; conversely, the mid and distal portions show large increases in stress, suggesting a distal shift in the loadpath. The Mayo design conserves bone mass, while simultaneously reducing the incidence of stress-shielding compared to the Alloclassic, revealing a key benefit of the distinctive geometry. Several important factors for stable fixation, determined in clinical evaluations from the literature, were evident in both designs: high levels of proximal bone loss and distal bone densification. The results suggest this novel computational framework can be utilized for comparative hip prosthesis shape, uniquely considering the post-operation bone remodeling as a design criterion.     

Analysis of a femoral hip prosthesis designed to reduce stress shielding

The natural stress distribution in the femur is significantly altered after total hip arthroplasty (THA). When an implant is introduced, it will carry a portion of the load, causing a reduction of stress in some regions of the remaining bone. This phenomenon is commonly known as stress shielding. In response to the changed mechanical environment the shielded bone will remodel according to Wolff's law, resulting in a loss of bone mass through the biological process called resorption. Resorption can, in turn, cause or contribute to loosening of the prosthesis. The problem is particularly common among younger THA recipients. This study explores the hypothesis that through redesign, a total hip prosthesis can be developed to substantially reduce stress shielding. First, we describe the development of a new femoral hip prosthesis designed to alleviate this problem through a new geometry and system of proximal fixation. A numerical comparison with a conventional intramedullary prosthesis as well as another proximally fixed prosthesis, recently developed by Munting and Verhelpen (1995. Journal of Biomechanics 28(8), 949–961) is presented. The results show that the new design produces a more physiological stress state in the proximal femur.     

骨科假体表面涂层技术现状

由植入物界面处的相互作用可能引起无菌性松动和假体周围感染。而无菌性松动和假体周围感染仍然是一个难以治疗的问题,并且最终可能导致假体植入失败,引起严重后果。理想的植入物应能促进骨整合,防止细菌粘附,减少细菌感染。骨科植入技术主要基于生物材料的开发和使用,随着材料科学和细胞生物学的发展,已可以用新的植入物表面涂层的进展来解决这些问题。本文回顾总结了时下骨科常见的假体涂层设计和相关问题,以期为进一步研究提供借鉴。     

Adaptive glenoid bone remodeling simulation

Glenoid prosthesis loosening is the most common cause for revision total shoulder arthroplasty. Stress-induced bone remodeling may compromise long-term prosthesis fixation and significantly contribute to loosening. Realistic, robust analysis of bone-prosthesis constructs need to look beyond initial post-implantation mechanics provided by static finite element (FE) simulation. Adaptive bone remodeling simulations based on Wolff's law are needed for evaluating long-term glenoid prostheses fixation. The purpose of this study was to take a first step towards this goal and create and validate two-dimensional FE simulations, using the intact glenoid, for computing subject-specific adaptive glenoid remodeling. Two-dimensional glenoid FE models were created from scapulae computed tomography images. Two distinct processes, “element” and “node” simulations, used the forward-Euler method to compute bone remodeling. Initial bone density was homogeneous. Center and offset load combinations were iteratively applied. To validate the simulations we performed location-specific statistical comparisons between predicted and actual bone density, load combinations, and “element” and “node” processes. Visually and quantitatively “element” simulations produced better results (p>0.22), and correlation coefficients ranged 0.51–0.69 (p<0.001). Having met this initial work's goals, we expect subject-specific FE glenoid bone remodeling simulations together with static FE stress analyses to be effective tools for designing and evaluating glenoid prostheses.     

A model of bone adaptation as an optimization process

The internal bone adaptation of the proximal femur is considered. A three-dimensional finite element model of the proximal femur is used. The bone remodeling in this work is numerically described by an evolutionary remodeling scheme with anisotropic material parameters and time-dependent loading. The memory of past loading is included in the model to account for the delay in the bone response from the load changes. The remodeling rate equation is derived from the structural optimization task of maximizing the stiffness in each time step. Additional information can be extracted from the optimization process; the remodeling equilibrium parameter where no apposition or resorption takes place, is defined as the element optimality conditions and the optimal design is used as an initial design for the onset of the remodeling simulation. Two examples of bone adaptation resulting from load changes are given, and the irreversible nature of the model is illustrated.     

Predicting bone remodeling around tissue- and bone-level dental implants used in reduced bone width

The objective of this study was to predict time-dependent bone remodeling around tissue- and bone-level dental implants used in patients with reduced bone width. The remodeling of bone around titanium tissue-level, and titanium and titanium–zirconium alloy bone-level implants was studied under 100 N oblique load for one month by implementing the Stanford theory into three-dimensional finite element models. Maximum principal stress, minimum principal stress, and strain energy density in peri-implant bone and displacement in x- and y- axes of the implant were evaluated. Maximum and minimum principal stresses around tissue-level implant were higher than bone-level implants and both bone-level implants experienced comparable stresses. Total strain energy density in bone around titanium implants slightly decreased during the first two weeks of loading followed by a recovery, and the titanium–zirconium implant showed minor changes in the axial plane. Total strain energy density changes in the loading and contralateral sides were higher in tissue-level implant than other implants in the cortical bone at the horizontal plane. The displacement values of the implants were almost constant over time. Tissue-level implants were associated with higher stresses than bone-level implants. The time-dependent biomechanical outcome of titanium–zirconium alloy bone-level implant was comparable to the titanium implant.     

Prediction of bone density around orthopedic implants delivering bisphosphonate

The fixation of an orthopedic implant depends strongly upon its initial stability. Peri-implant bone may resorb shortly after the surgery. This resorption is directly followed by new bone formation and implants fixation strengthening, the so-called secondary fixation. If the initial stability is not reached, the resorption continues and the implant fixation weakens, which leads to implant loosening. Studies with rats and dogs have shown that a solution to prevent peri-implant resorption is to deliver bisphosphonate from the implant surface.The aims of the study were, first, to develop a model of bone remodeling around an implant delivering bisphosphonate, second, to predict the bisphosphonate dose that would induce the maximal peri-implant bone density, and third to verify in vivo that peri-implant bone density is maximal with the calculated dose.The model consists of a bone remodeling equation and a drug diffusion equation. The change in bone density is driven by a mechanical stimulus and a drug stimulus. The drug stimulus function and the other numerical parameters were identified from experimental data. The model predicted that a dose of 0.3 μg of zoledronate on the implant would induce a maximal bone density. Implants with 0.3 μg of zoledronate were then implanted in rat femurs for 3, 6 and 9 weeks. We measured that peri-implant bone density was 4% greater with the calculated dose compared to the dose empirically described as best.The approach presented in this paper could be used in the design and analysis processes of experiments in local delivery of drug such as bisphosphonate.     

Mechanical analysis of a rodent segmental bone defect model: The effects of internal fixation and implant stiffness on load transfer

Segmental bone defect animal models are often used for evaluating the bone regeneration performance of bone substituting biomaterials. Since bone regeneration is dependent on mechanical loading, it is important to determine mechanical load transfer after stabilization of the defect and to study the effects of biomaterial stiffness on the transmitted load. In this study, we assess the mechanical load transmitted over a 6 mm femur defect that is stabilized with an internal PEEK fixation plate. Subsequently, three types of selective laser melted porous titanium implants with different stiffness values were used to graft the defect (five specimens per group). In one additional group, the defect was left empty. Micro strain gauges were used to measure strain values at four different locations of the fixation plate during external loading on the femoral head. The load sharing between the fixation plate and titanium implant was highly variable with standard deviations of measured strain values between 31 and 93% of the mean values. As a consequence, no significant differences were measured between the forces transmitted through the titanium implants with different elastic moduli. Only some non-significant trends were observed in the mean strain values that, consistent with the results of a previous finite element study, implied the force transmitted through the implant increases with the implant stiffness. The applied internal fixation method does not standardize mechanical loading over the defect to enable detecting small differences in bone regeneration performances of bone substituting biomaterials. In conclusion, the fixation method requires further optimization to reduce the effects of the operative procedure and make the mechanical loading more consistent and improve the overall sensitivity of this rat femur defect model.     

Design of intramedullary prostheses to prevent bone loss: predictions based on damage-stimulated remodelling

The adaptation of bone around intramedullary fixated prostheses, such as the femoral component of the hip joint or the radial component of the elbow joint, is well documented in follow-up studies. Bone adaptation takes the form of proximal bone atrophy accompanied, in some cases, by distal bone hypertrophy. A mechanistic model has been formulated to predict bone adaptation based on the concept that the continuous processes of damage and repair regulate bone adaptation. We apply the model to investigate the significance of two features of intramedullary prosthesis design on bone adaptation: prosthesis Young's modulus and the presence of a prosthesis collar. Results, as well as indicating some characteristics of accumulative-damage stimulated bone adaptation, predict that a low Young's modulus stem will very much reduce the extent of bone loss whereas the presence of a collar will have no significant effect. The results predict that a collarless low stiffness prosthesis is one possible approach for improving the secondary stability of intramedullary-fixated orthopaedic implants.     

Strain shielding in distal femur after patellofemoral arthroplasty under different activity conditions

Strain shielding, a mechanical effect occurring in structures combining stiff with more flexible materials, is considered to lead to a reduction of density in bone surrounding the implant. This effect can be related to the weakness of the implant fixation, which can promote implant loosening. Several studies describe a significant decrease in postoperative bone mineral density adjacent to joint implants, which can compromise their long-term fixation. The aim of the present study was to quantify the strain shielding effect on the distal femur after patellofemoral arthroplasty. For this purpose three activities of daily living were considered: level walking, stair climbing and deep bending at different angles of knee flexion. To determine the strain shielding effect, cortical bone strains were measured experimentally with triaxial strain gauges in synthetic femurs before and after patellofemoral arthroplasty for each of the different daily activities. The results showed that the patellofemoral arthroplasty in general reduced the strains in the medial and distal regions of the femur when deep bending activity occurred, consequently, strain shielding in these regions, with strain decreases of ?72.0% and ?67.5% were measured. On the other side, higher values of strain were found in the anterior region after patellofemoral replacement for this activity with an increase of +182.0%. The occurrence of strain shielding seems to be more significant when the angle of knee flexion and applied load increases. Strain shielding and over-loading may have relevant effects on bone remodeling surrounding the patellofemoral implant, suggesting a potential effect of later bone resorption in the medial and distal femur regions in case of regular deep bending activity.     

3 种不同弹性模量钛合金股骨假体在羊股骨置换模型中应力分布的三维 有限元分析

目的:通过三维有限元分析方法来观察并比较3种不同弹性模量钛合金股骨假体在羊股骨置换模型中von-Mises应力分布的情况。方法:采用64排螺旋CT对一健康成年羊的下肢股骨进行全长的CT扫描,扫描层厚为0.5 mm,扫描所得的数据存储为DICOM文件。将得到的DICOM文件导入到CT图像分析软件Mimics 10.0,然后利用Mimics 10.0软件来生成股骨的骨质点云数据,再将生成的骨质点云数据导入到Simpleware分析软件,通过机械加工反求中的复杂曲面造型技术建立起精确的三维实体模型。对三维实体模型进行网格划分,确定了髓腔的形状,并根据羊下肢股骨髓腔的形状设计了作者实验用的羊股骨假体模型,然后在ANSYS 12.1软件中进行网格划分。给予加载缓慢行走载荷以及扭转载荷,分析并比较羊股骨以及3种不同弹性模量钛合金股骨假体在股骨置换模型中von-Mises应力分布的情况。结果:在缓慢行走载荷以及扭转载荷条件下,3种不同弹性模量钛合金股骨假体von-Mises应力分布变化趋势一致,假体的柄颈结合部以及假体柄上1/3为应力集中区域。3种不同弹性模量的最大应力集中点均位于柄颈结合部,60 GPa弹性模量的股骨假体植入后假体的最大应力最小(37.8 MPa、29.1 MPa),股骨的最大应力最大(12.6 MPa、24.5 MPa);80 GPa的次之,假体的最大应力(38.4 MPa、33.4 MPa),股骨的最大应力(12.5 MPa、24.5 MPa);110 GPa的股骨假体植入后假体的最大应力最大(38.9 MPa、38.1 MPa),股骨的最大应力最小(12.3 MPa、24.5 MPa)。60 GPa弹性模量的股骨假体植入后的假体最大位移和相对位移均最小(缓慢行走载荷下假体最大位移为0.551 mm、相对位移为0.008 mm,扭转载荷下假体最大位移为0.730 mm、相对位移为0.011 mm)。结论:较低弹性模量的钛合金股骨假体(60 GPa)由于其弹性模量更接近于骨组织的弹性模量,股骨假体与股骨间的"应力遮挡"效应较小,更有利于应力在股骨假体及股骨间的传递,增加了股骨假体的早期稳定性,延长了其临床寿命。     

Computational simulation of simultaneous cortical and trabecular bone change in human proximal femur during bone remodeling

In this study, we developed a numerical framework that computationally determines simultaneous and interactive structural changes of cortical and trabecular bone types during bone remodeling, and we investigated the structural correlation between the two bone types in human proximal femur. We implemented a surface remodeling technique that performs bone remodeling in the exterior layer of the cortical bone while keeping its interior area unchanged. A micro-finite element (μFE) model was constructed that represents the entire cortical bone and full trabecular architecture in human proximal femur. This study simulated and compared the bone adaptation processes of two different structures: (1) femoral bone that has normal cortical bone shape and (2) perturbed femoral bone that has an artificial bone lump in the inferomedial cortex. Using the proposed numerical method in conjunction with design space optimization, we successfully obtained numerical results that resemble actual human proximal femur. The results revealed that actual cortical bone, as well as the trabecular bone, in human proximal femur has structurally optimal shapes, and it was also shown that a bone abnormality that has little contribution to bone structural integrity tends to disappear. This study also quantitatively determined the structural contribution of each bone: when the trabecular adaptation was complete, the trabecular bone supported 54% of the total load in the human proximal femur while the cortical bone carried 46%.     

Effects of Materials of Cementless Femoral Stem on the Functional Adaptation of Bone

The objective of this paper is to identify the effects of materials of cementless femoral stem on the functional adaptive behaviors of bone.The remodeling behaviors of a two-dimensional simplified model of cementless hip prosthesis with stiff stem,flexible 'iso-elastic' stem,one-dimensional Functionally Graded Material (FGM) stem and two-dimensional FGM stem for the period of four years after prosthesis replacement were quantified by incorporating the bone remodeling algorithm with finite element analysis.The distributions of bone density,von Mises stress,and interface shear stress were obtained.The results show that two-dimensional FGM stem may produce more mechanical stimuli and more uniform interface shear stress compared with the stems made of other materials,thus the host bone is well preserved.Accordingly,the two-dimensional FGM stem is an appropriate femoral implant from a biomechanical point of view.The numerical simulation in this paper can provide a quantitative computational paradigm for the changes of bone morphology caused by implants,which can help to improve the design of implant to reduce stress shielding and the risk of bone-prosthesis interface failure.     

Load-shift--numerical evaluation of a new design philosophy for uncemented hip prostheses

All hip replacement prostheses alter the load transfer from the hip joint into the femur by changing the mechanical loading of the proximal femur from an externally to an internally loaded system. This alteration of the load transfer causes stress shielding and might lead to severe bone loss, especially with uncemented prostheses. To minimize these effects, load transfer to the femur should occur as proximal as possible. In order to support sufficient primary stability, however, directly post operative (PO) distal stabilization is reasonable. Consequently, the prostheses have to alter its mechanical characteristics after implantation. This concept is referred to as load-shift concept. Primary stability during the early PO state is provided by a prosthesis shaft, which is widened at the tip by a biodegradable pin. After resorption of the pin load transfer occurs no longer distally. The objective of this study was the numerical evaluation of the load-shift concept. The analysis was performed with a finite element model. Three-dimensional non-linear dynamic gait analyses data were used to evaluate whether the load transfer during walking can be altered effectively by insertion and resorption of a distal pin. Directly PO 38% of the transverse load is transferred through the prosthesis shaft and micromotion of the proximal prostheses tip is below 55 microm. After resorption of the pin, no transverse loads are transferred through the prosthesis shaft. Therefore, the loading of the proximal bone tissue is far more pronounced than in the case of a standard prosthesis, demonstrating the feasibility of the load-shift concept. A balanced degradation of the pin simultaneously with the ingrowth of the prosthesis is expected to reduce hip replacement complications.     

Evaluation of bone remodeling around single dental implants of different lengths: a mechanobiological numerical simulation and validation using clinical data

Algorithmic models have been proposed to explain adaptive behavior of bone to loading; however, these models have not been applied to explain the biomechanics of short dental implants. Purpose of present study was to simulate bone remodeling around single implants of different lengths using mechanoregulatory tissue differentiation model derived from the Stanford theory, using finite elements analysis (FEA) and to validate the theoretical prediction with the clinical findings of crestal bone loss. Loading cycles were applied on 7-, 10-, or 13-mm-long dental implants to simulate daily mastication and bone remodeling was assessed by changes in the strain energy density of bone after a 3, 6, and 12 months of function. Moreover, clinical findings of marginal bone loss in 45 patients rehabilitated with same implant designs used in the simulation (n = 15) were computed to validate the theoretical results. FEA analysis showed that although the bone density values reduced over time in the cortical bone for all groups, bone remodeling was independent of implant length. Clinical data showed a similar pattern of bone resorption compared with the data generated from mathematical analyses, independent of implant length. The results of this study showed that the mechanoregulatory tissue model could be employed in monitoring the morphological changes in bone that is subjected to biomechanical loads. In addition, the implant length did not influence the bone remodeling around single dental implants during the first year of loading.     

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 \(\upmu \hbox {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.     

Stress analysis and failure prediction in the proximal femur before and after total hip replacement

An investigation was performed to determine the effects of the presence of two lengths of proximal Müller prosthesis on predicted failure loads, as compared to those for an intact femur. Three-dimensional stresses in a bone/cement/prosthesis system were determined using finite element methods, with both isotropic and transversely isotropic material properties used for the diaphyseal cortex. Significant increases in prosthesis stem stresses were found when the transversely isotropic material properties were employed in the diaphyseal cortex. This leads to the conclusion that accurate anisotropic material properties for bone are essential for precise stress determination and optimum design in prosthetic implants. Failure loads were also predicted for vertical compression and axial torque, similar to available experimental conditions, and were within the range of the experimental failure data found in the literature. The technique developed herein can be used to systematically assess existing as well as future implant designs, taking into account the complex three-dimensional interaction effects of the overall bone/cement/prosthesis system.     

The effect of muscle loading on the simulation of bone remodelling in the proximal femur

A large number of finite element analyses of the proximal femur rely on a simplified set of muscle and joint contact loads to represent the boundary conditions of the model. In the context of bone remodelling analysis around hip implants, muscle loading affects directly the spatial distribution of the remodelling signal. In the present study we performed a sensitivity analysis on the effect of different muscle loading configurations on the outcome of the bone remodelling simulation. An anatomical model of the femur with the implanted stem in place was constructed using the CT data of the Visible Human Project dataset of the National Institute of Health. The model was loaded with three muscle force configurations with increasing level of complexity. A strain adaptive remodelling rule was employed to simulate the post-operative bone changes around the implant stem and the results of the simulation were assessed quantitatively in terms of the bone mineral content changes in 18 periprosthetic regions of interest. The results showed considerable differences in the amount of bone loss predicted between the three cases. The simplified models generally predicted more pronounced bone loss. Although the overall remodelling patterns observed were similar, the bone conserving effect of additional muscle forces in the vicinity of their areas of attachment was clear. The results of this study suggest that the loading configuration of the FE model does play an important role in the outcome of the remodelling simulation.     

Predictions of bone remodeling around dental implant systems

This study presents the implementation of a mathematical bone remodeling algorithm to bone adaptation in the premolar area of the mandible around various dental implant systems, and thus sheds a new perspective to the complex interactions in dental implant mechanics. A two-dimensional, plane strain model of the bone was built from a CT-scan. The effect of implant contour on internal bone remodeling was investigated by considering four dental implant systems with contours similar to commercially available ones and another four with cylindrical and conical cross-sections. The remodeling algorithm predicts non-homogeneous density/elastic modulus distribution; and, implant contour has some effect on how this is distributed. Bone density is predicted to increase on the tips of the threads of the implants, but to decrease inside the grooves. Threadless implants favor to develop a softer bone around their periphery, compared to implant systems that have threads. The overall contour (dimensions and the shape) of an implant affect the bone density redistribution, but the differences between different implant systems are relatively small.