Activity intensity,assistive devices and joint replacement influence predicted remodelling in the proximal femur

Bone morphology and density changes are commonly observed following joint replacement, may contribute to the risks of implant loosening and periprosthetic fracture and reduce the available bone stock for revision surgery. This study was presented in the ‘Bone and Cartilage Mechanobiology across the scales’ WCCM symposium to review the development of remodelling prediction methods and to demonstrate simulation of adaptive bone remodelling around hip replacement femoral components, incorporating intrinsic (prosthesis) and extrinsic (activity and loading) factors. An iterative bone remodelling process was applied to finite element models of a femur implanted with a cementless total hip replacement (THR) and a hip resurfacing implant. Previously developed for a cemented THR implant, this modified process enabled the influence of pre- to post-operative changes in patient activity and joint loading to be evaluated. A control algorithm used identical pre- and post-operative conditions, and the predicted extents and temporal trends of remodelling were measured by generating virtual X-rays and DXA scans. The modified process improved qualitative and quantitative remodelling predictions for both the cementless THR and resurfacing implants, but demonstrated the sensitivity to DXA scan region definition and appropriate implant–bone position and sizing. Predicted remodelling in the intact femur in response to changed activity and loading demonstrated that in this simplified model, although the influence of the extrinsic effects were important, the mechanics of implantation were dominant. This study supports the application of predictive bone remodelling as one element in the range of physical and computational studies, which should be conducted in the preclinical evaluation of new prostheses.     

Long-term study of bone remodelling after femoral stem: a comparison between dexa and finite element simulation

A hip replacement with a cemented or cementless femoral stem produces an effect on the bone called adaptive remodelling, attributable to mechanical and biological factors. The objective of all of cementless prostheses designs has been to achieve a perfect transfer of loads in order to avoid stress-shielding, which produces an osteopenia. In order to quantify this, the long term and mass-produced study with dual energy X-ray absorptiometry (DEXA) is necessary. Finite element (FE) simulation makes possible the explanation of the biomechanical changes which are produced in the femur after stem implantation. The good correlation obtained between the results of the FE simulation and the densitometric study allow, on one hand, to explain from the point of view of biomechanical performance the changes observed in bone density in the long-term, where it is clear that these are due to a different transfer of load in the implanted model compared to the healthy femur; on the other hand, it validates the simulation model, in a way that it can be used in different conditions and at different time periods, to carry out a sufficiently precise prediction of the evolution of the bone density from the biomechanical behaviour in the interaction between the prosthesis and femur.     

Numerical simulation of the remodelling process of trabecular architecture around dental implants

Dental implants may alter the mechanical environment in the jawbone, thereby causing remodelling and adaptation of the surrounding trabecular bone tissues. To improve the efficacy of dental implant systems, it is necessary to consider the effect of bone remodelling on the performance of the prosthetic systems. In this study, finite element simulations were implemented to predict the evolution of microarchitecture around four implant systems using a previously developed model that combines both adaptive and microdamage-based mechano-sensory mechanisms in bone remodelling process. Changes in the trabecular architecture around dental implants were mainly focused. The simulation results indicate that the orientational and ladder-like architecture around the implants predicted herein is in good agreement with those observed in animal experiments and clinical observations. The proposed algorithms were shown to be effective in simulating the remodelling process of trabecular architecture around dental implant systems. In addition, the architectural features around four typical dental implant systems in alveolar bone were evaluated comparatively.     

A method of quantification of stress shielding in the proximal femur using hierarchical computational modeling

Stress shielding is a biomechanical phenomenon causing adaptive changes in bone strength and stiffness around metallic implants, which potentially lead to implant loosening. Accordingly, there is a need for standard, objective engineering measures of the “stress shielding” performances of an implant that can be employed in the process of computer-aided implant design. To provide and test such measures, we developed hierarchical computational models of adaptation of the trabecular microarchitecture at different sites in the proximal femur, in response to insertion of orthopaedic screws and in response to hypothetical reductions in hip joint and gluteal muscle forces. By identifying similar bone adaptation outcomes from the two scenarios, we were able to quantify the stress shielding caused by screws in terms of analogous hypothetical reductions in hip joint and gluteal muscle forces. Specifically, we developed planar lattice models of trabecular microstructures at five regions of interest (ROI) in the proximal femur. The homeostatic and abnormal loading conditions for the lattices were determined from a finite element model of the femur at the continuum scale and fed to an iterative algorithm simulating the adaptation of each lattice to these loads. When screws were inserted to the femur model, maximal simulated bone loss (17% decrease in apparent density, 10% decrease in thickness of trabeculae) was at the greater trochanter and this effect was equivalent to the effect of 50% reduction in gluteal force and normal hip joint force. We conclude that stress shielding performances can be quantified for different screw designs using model-predicted hypothetical musculoskeletal load fractions that would cause a similar pattern and extent of bone loss to that caused by the implants.     

Relation between subject-specific hip joint loading, stress distribution in the proximal femur and bone mineral density changes after total hip replacement

In the prediction of bone remodelling processes after total hip replacement (THR), modelling of the subject-specific geometry is now state-of-the-art. In this study, we demonstrate that inclusion of subject-specific loading conditions drastically influences the calculated stress distribution, and hence influences the correlation between calculated stress distributions and changes in bone mineral density (BMD) after THR.For two patients who received cementless THR, personalized finite element (FE) models of the proximal femur were generated representing the pre- and post-operative geometry. FE analyses were performed by imposing subject-specific three-dimensional hip joint contact forces as well as muscle forces calculated based on gait analysis data. Average values of the von Mises stress were calculated for relevant zones of the proximal femur. Subsequently, the load cases were interchanged and the effect on the stress distribution was evaluated. Finally, the subject-specific stress distribution was correlated to the changes in BMD at 3 and 6 months after THR.We found subject-specific differences in the stress distribution induced by specific loading conditions, as interchanging of the loading also interchanged the patterns of the stress distribution. The correlation between the calculated stress distribution and the changes in BMD were affected by the two-dimensional nature of the BMD measurement.Our results confirm the hypothesis that inclusion of subject-specific hip contact forces and muscle forces drastically influences the stress distribution in the proximal femur. In addition to patient-specific geometry, inclusion of patient-specific loading is, therefore, essential to obtain accurate input for the analysis of stress distribution after THR.     

A method of quantification of stress shielding in the proximal femur using hierarchical computational modeling

Stress shielding is a biomechanical phenomenon causing adaptive changes in bone strength and stiffness around metallic implants, which potentially lead to implant loosening. Accordingly, there is a need for standard, objective engineering measures of the "stress shielding" performances of an implant that can be employed in the process of computer-aided implant design. To provide and test such measures, we developed hierarchical computational models of adaptation of the trabecular microarchitecture at different sites in the proximal femur, in response to insertion of orthopaedic screws and in response to hypothetical reductions in hip joint and gluteal muscle forces. By identifying similar bone adaptation outcomes from the two scenarios, we were able to quantify the stress shielding caused by screws in terms of analogous hypothetical reductions in hip joint and gluteal muscle forces. Specifically, we developed planar lattice models of trabecular microstructures at five regions of interest (ROI) in the proximal femur. The homeostatic and abnormal loading conditions for the lattices were determined from a finite element model of the femur at the continuum scale and fed to an iterative algorithm simulating the adaptation of each lattice to these loads. When screws were inserted to the femur model, maximal simulated bone loss (17% decrease in apparent density, 10% decrease in thickness of trabeculae) was at the greater trochanter and this effect was equivalent to the effect of 50% reduction in gluteal force and normal hip joint force. We conclude that stress shielding performances can be quantified for different screw designs using model-predicted hypothetical musculoskeletal load fractions that would cause a similar pattern and extent of bone loss to that caused by the implants.     

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.

     

Concept and development of an orthotropic FE model of the proximal femur

PURPOSE: In contrast to many isotropic finite-element (FE) models of the femur in literature, it was the object of our study to develop an orthotropic FE "model femur" to realistically simulate three-dimensional bone remodelling. METHODS: The three-dimensional geometry of the proximal femur was reconstructed by CT scans of a pair of cadaveric femurs at equal distances of 2mm. These three-dimensional CT models were implemented into an FE simulation tool. Well-known "density-determined" bony material properties (Young's modulus; Poisson's ratio; ultimate strength in pressure, tension and torsion; shear modulus) were assigned to each FE of the same "CT-density-characterized" volumetric group.In order to fix the principal directions of stiffness in FE areas with the same "density characterization", the cadaveric femurs were cut in 2mm slices in frontal (left femur) and sagittal plane (right femur). Each femoral slice was scanned into a computer-based image processing system. On these images, the principal directions of stiffness of cancellous and cortical bone were determined manually using the orientation of the trabecular structures and the Haversian system. Finally, these geometric data were matched with the "CT-density characterized" three-dimensional femur model. In addition, the time and density-dependent adaptive behaviour of bone remodelling was taken into account by implementation of Carter's criterion. RESULTS: In the constructed "model femur", each FE is characterized by the principal directions of the stiffness and the "CT-density-determined" material properties of cortical and cancellous bone. Thus, on the basis of anatomic data a three-dimensional FE simulation reference model of the proximal femur was realized considering orthotropic conditions of bone behaviour. CONCLUSIONS: With the orthotropic "model femur", the fundamental basis has been formed to realize realistic simulations of the dynamical processes of bone remodelling under different loading conditions or operative procedures (osteotomies, total hip replacements, etc).     

A multibody knee model with discrete cartilage prediction of tibio-femoral contact mechanics

Combining musculoskeletal simulations with anatomical joint models capable of predicting cartilage contact mechanics would provide a valuable tool for studying the relationships between muscle force and cartilage loading. As a step towards producing multibody musculoskeletal models that include representation of cartilage tissue mechanics, this research developed a subject-specific multibody knee model that represented the tibia plateau cartilage as discrete rigid bodies that interacted with the femur through deformable contacts. Parameters for the compliant contact law were derived using three methods: (1) simplified Hertzian contact theory, (2) simplified elastic foundation contact theory and (3) parameter optimisation from a finite element (FE) solution. The contact parameters and contact friction were evaluated during a simulated walk in a virtual dynamic knee simulator, and the resulting kinematics were compared with measured in vitro kinematics. The effects on predicted contact pressures and cartilage–bone interface shear forces during the simulated walk were also evaluated. The compliant contact stiffness parameters had a statistically significant effect on predicted contact pressures as well as all tibio-femoral motions except flexion–extension. The contact friction was not statistically significant to contact pressures, but was statistically significant to medial–lateral translation and all rotations except flexion–extension. The magnitude of kinematic differences between model formulations was relatively small, but contact pressure predictions were sensitive to model formulation. The developed multibody knee model was computationally efficient and had a computation time 283 times faster than a FE simulation using the same geometries and boundary conditions.     

Numerical analysis of an osseointegrated prosthesis fixation with reduced bone failure risk and periprosthetic bone loss

Currently available implants for direct attachment of prosthesis to the skeletal system after transfemoral amputation (OPRA system, Integrum AB, Sweden and ISP Endo/Exo prosthesis, ESKA Implants AG, Germany) show many advantages over the conventional socket fixation. However, restraining biomechanical issues such as considerable bone loss around the stem and peri-prosthetic bone fractures are present. To overcome these limiting issues a new concept of the direct intramedullary fixation was developed. We hypothesize that the new design will reduce the peri-prosthetic bone failure risk and adverse bone remodeling by restoring the natural load transfer in the femur. Generic CT-based finite element models of an intact femur and amputated bones implanted with 3 analyzed implants were created and loaded with a normal walking and a forward fall load. The strain adaptive bone remodeling theory was used to predict long-term bone changes around the implants and the periprosthetic bone failure risk was evaluated by the von Mises stress criterion. The results show that the new design provides close to physiological distribution of stresses in the bone and lower bone failure risk for the normal walking as compared to the OPRA and the ISP implants. The bone remodeling simulations did not reveal any overall bone loss around the new design, as opposed to the OPRA and the ISP implants, which induce considerable bone loss in the distal end of the femur. This positive outcome shows that the presented concept has a potential to considerably improve safety of the rehabilitation with the direct fixation implants.     

A contact model with ingrowth control for bone remodelling around cementless stems.

This work presents a computational model for bone remodelling around cementless stems. The problem is formulated as a material optimisation problem considering the bone and stem surfaces to be in contact. To emphasise the behaviour of the bone/stem interface, the computer model detects the existence of bone ingrowth during the remodelling; consequently, the contact conditions are changed for a better interface simulation. The trabecular bone is modelled as a strictly orthotropic material with equivalent properties computed by homogenisation. The distribution of bone relative density is obtained by the minimisation of a function that considers both the bone structural stiffness and the biological cost associated with metabolic maintenance of bone tissue. The situation of multiple load conditions is considered. The remodelling law, obtained from the necessary conditions for an optimum, is derived analytically from the optimisation problem and solved numerically using a suitable finite element mesh. The formulation is applied to an implanted femur. Results of bone density and ingrowth distribution are obtained for different coating conditions. Bone ingrowth does not occur over the entire coated surfaces. Indeed, we observed regions where separation or high relative displacement occurs that preclude bone ingrowth attachment. This prediction of the model is consistent with clinical observations of bone ingrowth. Thus, this model, which detect bone ingrowth and allow modification of the interface conditions, are useful for analysis of existing stems as well as design optimisation of coating extent and location on such stems.     

A Comparative Analysis of Different Treatments for Distal Femur Fractures using the Finite Element Method

The main objective of this work is the evaluation, by means of the finite element method (FEM) of the mechanical stability and long-term microstructural modifications in bone induced to three different kinds of fractures of the distal femur by three types of implants: the Condyle Plate, the less invasive stabilization system plate (LISS) and the distal femur nail (DFN). The displacement and the stress distributions both in bone and implants and the internal bone remodelling process after fracture and fixation are obtained and analysed by computational simulation. The main conclusions of this work are that distal femoral fractures can be treated correctly with the Condyle Plate, the LISS plate and the DFN. The stresses both in LISS and DFN implant are high especially around the screws. When respect to remodelling, the LISS produces an important resorption in the fractured region, while the other two implants do not strongly modify bone tissue microstructure.     

A comparative analysis of different treatments for distal femur fractures using the finite element method

The main objective of this work is the evaluation, by means of the finite element method (FEM) of the mechanical stability and long-term microstructural modifications in bone induced to three different kinds of fractures of the distal femur by three types of implants: the Condyle Plate, the less invasive stabilization system plate (LISS) and the distal femur nail (DFN). The displacement and the stress distributions both in bone and implants and the internal bone remodelling process after fracture and fixation are obtained and analysed by computational simulation. The main conclusions of this work are that distal femoral fractures can be treated correctly with the Condyle Plate, the LISS plate and the DFN. The stresses both in LISS and DFN implant are high especially around the screws. When respect to remodelling, the LISS produces an important resorption in the fractured region, while the other two implants do not strongly modify bone tissue microstructure.     

Bone remodelling around the tibia due to total ankle replacement: effects of implant material and implant–bone interfacial conditions

Abstract

One of the major causes of implant loosening is due to excessive bone resorption surrounding the implant due to bone remodelling. The objective of the study is to investigate the effects of implant material and implant–bone interface conditions on bone remodelling around tibia bone due to total ankle replacement. Finite element models of intact and implanted ankles were developed using CT scan data sets. Bone remodelling algorithm was used in combination with FE analysis to predict the bone density changes around the ankle joint. Dorsiflexion, neutral, and plantar flexion positions were considered, along with muscle force and ligaments. Implant–bone interfacial conditions were assumed as debonded and bonded to represent non-osseointegration and fully osseointegration at the porous coated surface of the implant. To investigate the effect of implant material, three finite element models having different material combinations of the implant were developed. For model 1, tibial and talar components were made of Co–Cr–Mo, and meniscal bearing was made of UHMWPE. For model 2, tibial and talar components were made of ceramic and meniscal bearing was made of UHMWPE. For model 3, tibial and talar components were made of ceramic and meniscal bearing was made of CFR-PEEK. Changes in implant material showed no significant changes in bone density due to bone remodelling. Therefore, ceramic appears to be a viable alternative to metal and CFR-PEEK can be used in place of UHMWPE. This study also indicates that proper bonding between implant and bone is essential for long-term survival of the prosthetic components.     

A Femur-Implant Model for the Prediction of Bone Remodeling Behavior Induced by Cementless Stem

Bone remodeling simulation is an effective tool for the prediction of long-term effect of implant on the bone tissue, as well as the selection of an appropriate implant in terms of architecture and material. In this paper, a finite element model of proximal femur was developed to simulate the structures of internal trabecular and cortical bones by incorporating quantitative bone functional adaptation theory with finite element analysis. Cementless stems made of titanium, two types of Functionally Graded Material (FGM) and flexible ‘iso-elastic’ material as comparison were implanted in the structure of proximal femur respectively to simulate the bone remodeling behaviors of host bone. The distributions of bone density, von Mises stress, and interface shear stress were obtained. All the prosthetic stems had effects on the bone remodeling behaviors of proximal femur, but the degrees of stress shielding were different. The amount of bone loss caused by titanium implant was in agreement with the clinical observation. The FGM stems caused less bone loss than that of the titanium stem, in which FGM I stem (titanium richer at the top to more HAP/Col towards the bottom) could relieve stress shielding effectively, and the interface shear stresses were more evenly distributed in the model with FGM I stem in comparison with those in the models with FGM II (titanium and bioglass) and titanium stems. The numerical simulations in the present study provided theoretical basis for FGM as an appropriate material of femoral implant from a biomechanical point of view. The next steps are to fabricate FGM stem and to conduct animal experiments to investigate the effects of FGM stem on the remodeling behaviors using animal model.     

Accuracy of finite element predictions in sideways load configurations for the proximal human femur

Subject-specific finite element models have been used to predict stress-state and fracture risk in individual patients. While many studies analysed quasi-axial loading configurations, only few works simulated sideways load configurations, such as those arising in a fall. The majority among these latter directly predicted bone strength, without assessing elastic strain prediction accuracy. The aim of the present work was to evaluate if a subject-specific finite element modelling technique from CT data that accurately predicted strains in quasi-axial loading configurations is suitable to accurately predict strains also when applying low magnitude loads in sideways configurations. To this aim, a combined numerical-experimental study was performed to compare finite element predicted strains with strain-gauge measurements from three cadaver proximal femurs instrumented with sixteen strain rosettes and tested non-destructively under twelve loading configurations, spanning a wide cone (0-30° for both adduction and internal rotation angles) of sideways fall scenarios. The results of the present study evidenced a satisfactory agreement between experimentally measured and predicted strains (R(2) greater than 0.9, RMSE% lower than 10%) and displacements. The achieved strain prediction accuracy is comparable to those obtained in state of the art studies in quasi-axial loading configurations. Still, the presence of the highest strain prediction errors (around 30%) in the lateral neck aspect would deserve attention in future studies targeting bone failure.     

Consideration of multiple load cases is critical in modelling orthotropic bone adaptation in the femur

Functional adaptation of the femur has been investigated in several studies by embedding bone remodelling algorithms in finite element (FE) models, with simplifications often made to the representation of bone’s material symmetry and mechanical environment. An orthotropic strain-driven adaptation algorithm is proposed in order to predict the femur’s volumetric material property distribution and directionality of its internal structures within a continuum. The algorithm was applied to a FE model of the femur, with muscles, ligaments and joints included explicitly. Multiple load cases representing distinct frames of two activities of daily living (walking and stair climbing) were considered. It is hypothesised that low shear moduli occur in areas of bone that are simply loaded and high shear moduli in areas subjected to complex loading conditions. In addition, it is investigated whether material properties of different femoral regions are stimulated by different activities. The loading and boundary conditions were considered to provide a physiological mechanical environment. The resulting volumetric material property distribution and directionalities agreed with ex vivo imaging data for the whole femur. Regions where non-orthogonal trabecular crossing has been documented coincided with higher values of predicted shear moduli. The topological influence of the different activities modelled was analysed. The influence of stair climbing on the properties of the femoral neck region is highlighted. It is recommended that multiple load cases should be considered when modelling bone adaptation. The orthotropic model of the complete femur is released with this study.     

Application of an anisotropic bone-remodelling model based on a damage-repair theory to the analysis of the proximal femur before and after total hip replacement

In this work, a new model for internal anisotropic bone remodelling is applied to the study of the remodelling behaviour of the proximal femur before and after total hip replacement (THR). This model considers bone remodelling under the scope of a general damage-repair theory following the principles of continuum damage mechanics. A "damage-repair" tensor is defined in terms of the apparent density and Cowin's "fabric tensor", respectively, associated with porosity and directionality of the trabeculae. The different elements of a thermodynamically consistent damage theory are established, including resorption and apposition criteria, evolution law and rate of remodelling. All of these elements were introduced and discussed in detail in a previous paper (García, J. M., Martinez, M. A., Doblaré, M., 2001. An anisotrophic internal-external bone adaptation model based on a combination of CAO and continuum damage mechanics technologies. Computer Methods in Biomechanics and Biomedical Engineering 4(4), 355-378.), including the definition of the proposed mechanical stimulus and the qualitative properties of the model. In this paper, the fundamentals of the proposed model are briefly reviewed and the computational aspects of its implementation are discussed. This model is then applied to the analysis of the remodelling behaviour of the intact femur obtaining densities and mass principal values and directions very close to the experimental data. The second application involved the proximal femoral extremity after THR and the inclusion of an Exeter prosthesis. As a result of the simulation process, some well-known features previously detected in medical clinics were recovered, such as the stress yielding effect in the proximal part of the implant or the enlargement of the cortical layer at the distal part of the implant. With respect to the anisotropic properties, bone microstructure and local stiffness are known to tend to align with the stress principal directions. This experimental fact is mathematically proved in the framework of this remodelling model and clearly shown in the results corresponding to the intact femur. After THR the degree of anisotropy decreases tending, specifically in the proximal femur, to a more isotropic behaviour.     

Patient-specific bone modelling and remodelling simulation of hypoparathyroidism based on human iliac crest biopsies

We previously developed a load-adaptive bone modelling and remodelling simulation model that can predict changes in the bone micro-architecture as a result of changes in mechanical loading or cell activity. In combination with a novel algorithm to estimate loading conditions, this offers the possibility for patient-specific predictions of bone modelling and remodelling. Based on such models, the underlying mechanisms of bone diseases and/or the effects of certain drugs and their influence on the bone micro-architecture can be investigated. In the present study we test the ability of this approach to predict changes in bone micro-architecture during hypoparathyroidism (HypoPT), as an illustrative example. We hypothesize that, apart from reducing bone turnover, HypoPT must also lead to increased osteocyte mechanosensitivity in order to explain the changes in bone mass seen in patients. Healthy human iliac crest biopsies were used as the starting point for the simulations that mimic HypoPT conditions and the resultant micro-architectures were compared to age-matched clinical HypoPT biopsies. Simulation results were in good agreement with the clinical data when osteocyte mechanosensitivity was increased by 40%. In conclusion, the results confirm our hypothesis, and also demonstrate that patient-specific bone modelling and remodelling simulations are feasible.     

Numerical simulation of load-induced bone structural remodelling using stress-limit criterion

A simple and efficient numerical method for predicting the remodelling of adaptive materials and structures under applied loading was presented and implemented within a finite element framework. The model uses the trajectorial architecture theory of optimisation to predict the remodelling of material microstructure and structural organisation under mechanical loading. We used the proposed model to calculate the density distribution of proximal femur in the frontal plane. The loading considered was the hip joint contact forces and muscular forces at the attachment sites of the muscles to the bone. These forces were estimated from a separate finite element calculation using a heterogeneous three-dimensional model of the proximal femur. The density distributions obtained by this procedure has a qualitative similarity with in vivo observations. Solutions displayed the characteristic high-density channels that are evident in the Dual X-ray Absorptiometry scan. There is also evidence of the intramedullary canal, as well as low-density regions in the femoral neck. Several parametric studies were carried out to highlight the advantages of the proposed method, which includes fast convergence and low-computational cost. The potential applications of the proposed method in predicting bone structural remodelling in cancer are also briefly discussed.