On the Use of Bone Remodelling Models to Estimate the Density Distribution of Bones. Uniqueness of the Solution

Trabecula-level mechanoadaptation: Numerical analysis of morphological changes

BACKGROUND

Bone is a living material that, unlike man-made ones, demonstrates continuous adaptation of its structure and mechanical properties to resist the imposed mechanical loading. Adaptation in trabecular bone is characterised by improvement of its stiffness in the loading direction and respective realignment of trabecular load-bearing architecture. Considerable experimental and simulation evidence of trabecular bone adaptation to its mechanical environment at the tissue- and organ-levels was obtained, while little attention was given to the trabecula-level of this process. This study aims to describe and classify load-driven morphological changes at the level of individual trabeculae and to propose their drivers.

METHOD

For this purpose, a well-established mechanoregulation-based numerical model of bone adaptation was implemented in a user-defined subroutine that changed the structural and mechanical properties of trabeculae based on the magnitude of a mechanical stimulus. This subroutine was used in conjunction with finite-element models of variously shaped structures representing trabeculae loaded in compression or shear.

RESULTS

In all analysed cases, trabeculae underwent morphological evolution under applied compressive or shear loading. Among twelve cases analysed, six main mechanisms of morphological evolution were established: reorientation, splitting, merging, full resorption, thinning, and thickening. Moreover, all simulated cases presented the ability to reduce the mean value of von Mises stress while increasing their ability to resist compressive/shear loading during adaptation.

CONCLUSION

This study evaluated morphological and mechanical changes in trabeculae of different shapes in response to compressive or shear loadings and compared them based on the analysis of von Mises stress distribution as well as profiles of normal and shear stresses in the trabeculae at different stages of their adaptation.

Computational Analysis of the Influence of Menopause and Ageing on Bone Mineral Density, Exploring the Impact of Bone Turnover and Focal Bone Balance-A Study on Overload and Underload Scenarios

This study aims to investigate the impact of hormonal imbalances during menopause, compounded by the natural ageing process, on bone health. Specifically, it examines the effects of increased bone turnover and focal bone balance on bone mass. A three-dimensional computational bone remodeling model was employed to simulate the response of the femur to habitual loads over a 19-year period, spanning premenopause, menopause, and postmenopause. The model was calibrated using experimental bone mineral density data from the literature to ensure accurate simulations. The study reveals that individual alterations in bone turnover or focal bone balance do not fully account for the observed experimental outcomes. Instead, simultaneous changes in both factors provide a more comprehensive explanation, leading to increased porosity while maintaining the material-to-apparent density ratio. Additionally, different load scenarios were tested, demonstrating that reaching the clinical osteoporosis threshold is independent of the timing of load changes. However, underload scenarios resulted in the threshold being reached approximately 6 years earlier than overload scenarios. These findings hold significant implications for strategies aimed at delaying the onset of osteoporosis and minimizing fracture risks through targeted mechanical stimulation during the early stages of menopause.

Stress-shielding resistant design of custom pelvic prostheses using lattice-based topology optimization

Endoprosthetic reconstruction of the pelvic bone using 3D-printed, custom-made implants has delivered early load-bearing ability and good functional outcomes in the short term to individuals with pelvic sarcoma. However, excessive stress-shielding and subsequent resorption of peri‑prosthetic bone can imperil the long-term stability of such implants. To evaluate the stress-shielding performance of pelvic prostheses, we developed a sequential modeling scheme using subject-specific finite element models of the pelvic bone-implant complex and personalized neuromusculoskeletal models for pre- and post-surgery walking. A new topology optimization approach is introduced for the stress-shielding resistant (SSR) design of custom pelvic prostheses, which uses 3D-printable porous lattice structures. The SSR optimization was applied to a typical pelvic prosthesis to reconstruct a type II+III bone resection. The stress-shielding performance of the optimized implant based on the SSR approach was compared against the conventional optimization. The volume of the peri‑prosthetic bone predicted to undergo resorption post-surgery decreased from 44 to 18%. This improvement in stress-shielding resistance was achieved without compromising the structural integrity of the prosthesis. The SSR design approach has the potential to improve the long-term stability of custom-made pelvic prostheses.

Computational model of articular cartilage regeneration induced by scaffold implantation in vivo

Computational models allow to explain phenomena that cannot be observed through an animal model, such as the strain and stress states which can highly influence regeneration of the tissue. For this purpose, we have developed a simulation tool to determine the mechanical conditions provided by the polymeric scaffold. The computational model considered the articular cartilage, the subchondral bone, and the scaffold. All materials were modeled as poroelastic, and the cartilage had linear-elastic oriented collagen fibers. This model was able to explain the remodeling process that subchondral bone goes through, and how the scaffold allowed the conditions for cartilage regeneration. These results suggest that the use of scaffolds might lead the cartilaginous tissue growth in vivo by providing a better mechanical environment. Moreover, the developed computational model demonstrated to be useful as a tool prior experimental in vivo studies, by predicting the possible outcome of newly proposed treatments allowing to discard approaches that might not bring good results.

Mechanoregulated trabecular bone adaptation: Progress report on in silico approaches

Adaptation is the process by which bone responds to changes in loading environment and modulates its properties and spatial organization to meet the mechanical demands. Adaptation in trabecular bone is achieved through increase in bone mass and alignment of trabecular-bone morphology along the loading direction. This transformation of internal microstructure is governed by mechanical stimuli sensed by mechanosensory cells in the bone matrix. Realisation of adaptation in the form of local bone-resorption and -formation activities as a function of mechanical stimuli is still debated. In silico modelling is a useful tool for simulation of various scenarios that cannot be investigated in vivo and particularly well suited for prediction of trabecular bone adaptation. This progress report presents the recent advances in in silico modelling of mechanoregulated adaptation at the scale of trabecular bone tissue. Four well-established bone-adaptation models are reviewed in terms of their recent improvements and validation. They consider various mechanical factors: (i) strain energy density, (ii) strain and damage, (iii) stress nonuniformity and (iv) daily stress. Contradictions of these models are discussed and their ability to describe adequately a real-life mechanoregulation process in bone is compared.

A coupled computational framework for bone fracture healing and long-term remodelling: Investigating the role of internal fixation on bone fractures

In this study, a coupled computational modelling framework for bone fracture repair is presented that enables predictions of both healing and remodelling phases of the fracture region and is used to investigate the role of an internal fixation plate on the long-term healing performance of a fracture tibia under a range of different conditions. It was found that introduction of a titanium plate allowed the tibia to undergo successful healing at higher loading conditions and fracture gaps, compared with the non-plated versions. While these plated cases showed faster rates of repair in the healing phase, their performance was substantially different once they entered the remodelling phase, with substantial regions of stress shielding predicted. This framework is one of the few implementations of both fracture healing and remodelling phases of bone repair and includes several innovative approaches to smoothing, time-averaging and time incrementation in its implementation, thereby avoiding any unwanted abrupt changes between tissue phenotypes. This provides a better representation of tissue development in the fracture site when compared with fracture healing models alone and provides a suitable platform to investigate the long-term performance of orthopaedic fixation devices. This would enable the more effective design of permanent fixation devices and optimisation of the spatial and temporal performance of bioabsorbable implants.

Computational Analysis of Bone Remodeling in the Proximal Tibia Under Electrical Stimulation Considering the Piezoelectric Properties

The piezoelectricity of bone is known to play a crucial role in bone adaptation and remodeling. The application of an external stimulus such as mechanical strain or electric field has the potential to enhance bone formation and implant osseointegration. Therefore, in the present study, the objective is to investigate bone remodeling under electromechanical stimulation as a step towards establishing therapeutic strategies. For the first time, piezoelectric bone remodeling in the human proximal tibia under electro-mechanical loads was analyzed using the finite element method in an open-source framework. The predicted bone density distributions were qualitatively and quantitatively assessed by comparing with the computed tomography (CT) scan and the bone mineral density (BMD) calculated from the CT, respectively. The effect of model parameters such as uniform initial bone density and reference stimulus on the final density distribution was investigated. Results of the parametric study showed that for different values of initial bone density the model predicted similar but not identical final density distribution. It was also shown that higher reference stimulus value yielded lower average bone density at the final time. The present study demonstrates an increase in bone density as a result of electrical stimulation. Thus, to minimize bone loss, for example, due to physical impairment or osteoporosis, mechanical loads during daily physical activities could be partially replaced by therapeutic electrical stimulation.

Finite element analysis of bone remodelling with piezoelectric effects using an open-source framework

Bone tissue exhibits piezoelectric properties and thus is capable of transforming mechanical stress into electrical potential. Piezoelectricity has been shown to play a vital role in bone adaptation and remodelling processes. Therefore, to better understand the interplay between mechanical and electrical stimulation during these processes, strain-adaptive bone remodelling models without and with considering the piezoelectric effect were simulated using the Python-based open-source software framework. To discretise numerical attributes, the finite element method (FEM) was used for the spatial variables and an explicit Euler scheme for the temporal derivatives. The predicted bone apparent density distributions were qualitatively and quantitatively evaluated against the radiographic scan of a human proximal femur and the bone apparent density calculated using a bone mineral density (BMD) calibration phantom, respectively. Additionally, the effect of the initial bone density on the resulting predicted density distribution was investigated globally and locally. The simulation results showed that the electrically stimulated bone surface enhanced bone deposition and these are in good agreement with previous findings from the literature. Moreover, mechanical stimuli due to daily physical activities could be supported by therapeutic electrical stimulation to reduce bone loss in case of physical impairment or osteoporosis. The bone remodelling algorithm implemented using an open-source software framework facilitates easy accessibility and reproducibility of finite element analysis made.

Remodelling of trabecular bone in human distal tibia: A model based on an in-vivo HR-pQCT study

An abnormal remodelling process of bones can lead to various bone disorders, such as osteoporosis, making them prone to fracture. Simulations of load-induced remodelling of trabecular bone were used to investigate its response to mechanical signal. However, the role of mechanostat in trabecular-bone remodelling has not yet been investigated in simulations underpinned by a longitudinal in-vivo study in humans. In this work, a finite-element model based on a 6-month longitudinal in-vivo HR-pQCT study was developed and validated to investigate the effect of mechanical stimuli on bone remodelling. The simulated changes in microstructural parameters and density of trabecular bone were compared with respective experimental results. A maximum principal strain (MPS) and a maximum principal strain gradient (∇MPS) were used as mechanical signals to drive a five-stage mechanostat remodelling model, including additional over-strain and damage stages. It was found that the density distribution varied with the studied mechanical signals, along with decreasing with time levels of bone volume fraction BV/TV, trabecular thickness Tb.Th and bone surface area Tb.BS as well as increased trabecular separation Tb.Sp. Among these parameters, BV/TV and Tb.Th together with the bone-remodelling parameters from the MPS model demonstrated a significant correlation with the experimental data. The developed model provides a good foundation for further development and investigation of the relationships between mechanical loading and human-bone microarchitecture.

Analysis of the uniqueness and stability of solutions to problems regarding the bone-remodeling process

Simulation of the bone remodeling process is extremely important because it makes possible the structure forecast of one or several bones when anomalous situations, such as prosthesis installation, occur. Thus, it is necessary that the mathematical model to simulate the bone remodeling process be reliable; that is, the numerical solution must be stable regardless of initial density field for a phenomenological approach to model the process. For several models found in the literature, this characteristic of stability is not observed, largely due to the discontinuities present in the property values of the models (e.g., Young's modulus and Poisson's ratio). In addition, checkerboard formation and the lazy zone prevent the uniqueness of the solution. To correct these difficulties, this study proposes a set of modifications to guarantee the uniqueness and stability of the solutions, when a phenomenological approach is used. The proposed modifications are: (a) change the rate of remodeling curve in the lazy zone region and (b) create transition functions to guarantee the continuity of the expressions used to describe Young's modulus and Poisson's ratio. Moreover, the stress smoothing process controls the checkerboard formation. Numerical analysis is used to simulate the solution behavior from each proposed modification. The results show that, when all proposed modifications are applied to the three-dimensional models simulated here, it is possible to observe the tendency toward a unique solution.

Numerical Simulation of Mandible Bone Remodeling under Tooth Loading: A Parametric Study

Bone adapts to the change of mechanical stimulus by bone remodeling activities. A number of numerical algorithms have been developed to model the adaptive bone remodeling under mechanical loads for orthopedic and dental applications. This paper examines the effects of several model parameters on the computed apparent bone density in mandible under normal chewing and biting forces. The density change rate was based on the strain energy density per unit mass. The algorithms used in this study containing an equilibrium zone (lazy zone) and saturated values of density change rate provides certain stability to result in convergence without discontinuous checkerboard patterns. The parametric study shows that when different boundary conditions were applied, the bone density distributions at convergence were very different, except in the vicinity of the applied loads. Compared with the effects of boundary conditions, the models are less sensitive to the choice of initial density values. Several models starting from different initial density values resulted in similar but not exactly the same bone density distribution at convergence. The results also show that higher reference value of mechanical stimulus resulted in lower average bone density at convergence. Moreover, the width of equilibrium zone did not substantially affect the average density at convergence. However, with increasing width, the areas with the highest and the lowest bone density areas were all reduced. The limitations of the models and challenges for future work were discussed for the better agreement between the computed results and the in vivo data.

The influence of an extra-articular implant on bone remodelling of the knee joint

Bone remodelling is a crucial feature of maintaining healthy bones. The loading conditions on the bones are one of the key aspects which affect the bone remodelling cycle. Many implants, such as hip and knee implants, affect the natural loading conditions and hence influence bone remodelling. Theoretical and numerical methods, such as adaptive bone remodelling, can be used to investigate how an implant affects bone mineral density (BMD). This research aimed to study the influence of an extra-articular implant on bone remodelling of the knee joint using adaptive bone remodelling. Initially, a finite element (FE) model of the knee joint was created. A user-defined material subroutine was developed to generate a heterogeneous BMD distribution in the FE model. The heterogeneous density was then assigned to the knee model with the implant in order to investigate how the implant would affect BMD of the knee joint, five years postoperatively. It was observed that in the medial compartments of the femur and tibia, bone mineral density increased by approximately 3.4% and 4.1%, respectively, and the density for the fixation holes of both bones increased by around 2.2%. From these results, it is concluded that implanting of this load-sharing device does not result in significantly adverse BMD changes in the femur and tibia.

A comparative analysis of internal bone remodelling concepts in a novel implant for direct skeletal attachment of limb prosthesis evaluation: A finite element analysis

Nowadays, numerous internal bone remodelling concepts are under development, in order to estimate long-term functionality of implants by evaluating the intensity of stress-shielding effect. This effect is also analysed for the implants for direct skeletal attachment, considered as a better exoprosthesis fixation method than prosthetic sockets. Most of bone remodelling approaches are based on basic concepts, differing with certain assumptions, which may affect the accuracy of the results. This article compares commonly used internal bone remodelling concepts and evaluates the functionality of the proposed Limb Prosthesis Osseointegrated Fixation System for direct skeletal attachment of limb prosthesis in comparison with two currently available implants: the Intraosseous Transcutaneous Amputation Prosthesis and the Osseointegrated Prostheses for the Rehabilitation of Amputees. Three concepts were chosen: without and with lazy zone and with the use of quadratic formula which considers bone overloading. Therefore, three finite element models were created with identical femur implanted with each of analysed implants. The implants were loaded with loads that refer to two stages of gait cycle (heel strike and toe-off). The analysed concepts have given similar results, allowing to assume that each of them can be successfully used to estimate internal bone remodelling around analysed implants for direct skeletal attachment of limb prosthesis. The results also present higher functionality of the proposed implant for direct skeletal attachment of limb prosthesis due to a significant reduction in stress-shielding in the analysed areas around implant in comparison with the Intraosseous Transcutaneous Amputation Prosthesis and the Osseointegrated Prostheses for the Rehabilitation of Amputees. It suggests that the proposed design is a better alternative to the currently used solutions.