An enhanced version of a bone-remodelling model based on the continuum damage mechanics theory

The purpose of this work was to propose an enhancement of Doblaré and García's internal bone remodelling model based on the continuum damage mechanics (CDM) theory. In their paper, they stated that the evolution of the internal variables of the bone microstructure, and its incidence on the modification of the elastic constitutive parameters, may be formulated following the principles of CDM, although no actual damage was considered. The resorption and apposition criteria (similar to the damage criterion) were expressed in terms of a mechanical stimulus. However, the resorption criterion is lacking a dimensional consistency with the remodelling rate. We propose here an enhancement to this resorption criterion, insuring the dimensional consistency while retaining the physical properties of the original remodelling model. We then analyse the change in the resorption criterion hypersurface in the stress space for a two-dimensional (2D) analysis. We finally apply the new formulation to analyse the structural evolution of a 2D femur. This analysis gives results consistent with the original model but with a faster and more stable convergence rate.     

Improving the damage accumulation in a biomechanical bone remodelling model

We extend, reformulate and analyse a phenomenological model for bone remodelling. The original macrobiomechanical model (MBM), proposed by Hazelwood et al. [J Biomech 2001; 34:299–308], couples a population equation for the cellular activities of the basic multicellular units (BMUs) in the bone and a rate equation to account for microdamage and repair. We propose to account for bone failure under severe overstressing by incorporating a Paris-like power-law damage accumulation term. The extended model agrees with the Hazelwood et al. predictions when the bone is under-stressed, and allows for suitably loaded bones to fail, in agreement with other MBM and experimental data regarding damage by fatigue. We numerically solve the extended model using a convergent algorithm and show that for unchanging loads, the stationary solution captures fully the model behaviour. We compute and analyse the stationary solutions. Our analysis helps guide additional extensions to this and other BMU activity based models.     

Informing phenomenological structural bone remodelling with a mechanistic poroelastic model

Studies suggest that fluid motion in the extracellular space may be involved in the cellular mechanosensitivity at play in the bone tissue adaptation process. Previously, the authors developed a mesoscale predictive structural model of the femur using truss elements to represent trabecular bone, relying on a phenomenological strain-based bone adaptation algorithm. In order to introduce a response to bending and shear, the authors considered the use of beam elements, requiring a new formulation of the bone adaptation drivers. The primary goal of the study presented here was to isolate phenomenological drivers based on the results of a mechanistic approach to be used with a beam element representation of trabecular bone in mesoscale structural modelling. A single-beam model and a microscale poroelastic model of a single trabecula were developed. A mechanistic iterative adaptation algorithm was implemented based on fluid motion velocity through the bone matrix pores to predict the remodelled geometries of the poroelastic trabecula under 42 different loading scenarios. Regression analyses were used to correlate the changes in poroelastic trabecula thickness and orientation to the initial strain outputs of the beam model. Linear (\(R^2>0.998\)) and third-order polynomial (\(R^2 >0.98\)) relationships were found between change in cross section and axial strain at the central axis, and between beam reorientation and ratio of bending strain to axial strain, respectively. Implementing these relationships into the phenomenological predictive algorithm for the mesoscale structural femur has the potential to produce a model combining biofidelic structure and mechanical behaviour with computational efficiency.     

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.     

A theoretical model of the effect of continuum damage on a bone adaptation model

Throughout life, bone is continuously turning over by the well-regulated processes of bone formation and resorption. Everyday activities damage bone, and this damage is normally repaired in a continuous remodelling process. When an imbalance in this remodelling process occurs, bones may become more susceptible to fracture. This paper is devoted to a theoretical modelling of the competition between damage and internal remodelling in bones. The general theory of adaptive damaged-elastic materials proposed here as a model for the physiological process of damaged-bone remodelling follows the general framework of continuum thermodynamics where new damaged-bone remodelling law and associated thermodynamical restrictions are stated, and specialized to the case of small strain in isothermal processes. An attempt is also made to derive: (a) the damage force (adaptive damage energy release rate ) which controls the microcracks propagation and arrest, and (b) the damage rule by introducing damage thresholds and loading/unloading conditions.     

An anisotropic internal-external bone adaptation model based on a combination of CAO and continuum damage mechanics technologies

In this work, a complete internal-external bone-remodelling scheme is presented and implemented into a finite element code. This model uses a combination of an anisotropic internal remodelling model based on a new Continuum "Damage-Repair" theory and an external adaptation approach that follows the idea, early introduced by Mattheck et al., to simulate the growth behaviour of biological systems, known as CAO method. This combined scheme qualitatively resembles most of the main features of the bone adaptive behaviour, like the bone mass distribution (heterogeneity and porosity), the directional internal structure (anisotropy), the alignment of the microstructure with the constitutive principal directions and these with those of the stress tensor when permanently loaded by a unique stress state (Wolff's law). It is also thermodynamically consistent, fulfilling a principle of minimum mechanical dissipation. Finally, the comparison between the predicted results and the ones obtained by different experimental tests allows us to conclude that this model is able of reproducing qualitatively the global behaviour of bone tissue when subjected to external mechanical loads.     

Allosteric control model of bone remodelling containing periodical modes

To help to understand the modelling process that occurs when a scaffold is implanted it is vital to understand the rather complex bone remodelling process prevalent in native bone. We have formulated a mathematical model that predicts osteoactivity both in scaffolds, as well as in bone in vivo and could set a basis for the more detailed allosteric models. The model is extended towards a bio-cybernetic vision of basic multicellular unit (BMU) action, when some of the regulation loops have been modified to reflect the allosteric control mechanisms, developed by Michaels-Menten, Hill, Koshland-Nemethy-Filmer, Monod-Wyman-Changeux. By implementation of this approach a four-dimensional system was obtained that shows steady cyclic behaviour using a wide range of constants with clear biological meaning. We have observed that a local steady state appears as a limiting cycle in multi-dimensional phase space and this is discussed in this paper. Physiological interpretation of this limiting four-dimension cycle possibly related to a conservative-like value has been proposed. Analysis and simulation of the model has shown an analogy between this conservative value, as a kind of substrate-energy regenerative potential of the bone remodelling system with a molecular nature, and to the classical physical value--energy. This dynamic recovery potential is directed against both mechanical and biomechanical damage to the bone. Furthermore, the current model has credibility when compared to the normal bone remodelling process. In the framework of widely recognised Hill mechanisms of allosteric regulation the cyclic attractor, described formerly for a pure cellular model, prevails for different forms of feedback control. This result indicates the viability of the proposed existence of a conservative value (analogous to energy) that characterises the recovery potential of the bone remodelling cycle. Linear stability analysis has been performed in order to determine the robustness of the basic state, however, additional work is required to study a wider range of constants.     

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 bone remodelling model including the directional activity of BMUs

Bone is able to adapt itself to the mechanical and biological environment by changing its porosity and/or orientation of its internal microstructure in a process known as bone remodelling. As a consequence, a change of bone mechanical properties is produced leading to an optimum structure, able to bear the external loads with the minimum weight. This adaptation is carried out by a temporal association of cells known as BMUs (basic multicellular units) that resorb old bone and sometimes produce new organic extracellular matrix (osteoid) that is later mineralized. This involves changes in porosity, damage level (density of microcracks accumulated by cyclic loads) and mineral content. All of these features were taken into account in a previous model, but the whole process and therefore the resulting bone constitutive behaviour was considered isotropic. The model proposed herein, recognizing that bone is actually anisotropic, tries to explain how BMUs modify the anisotropy by changing their progressing direction. We check the potential of the model to predict the alignment of the bone microstructure with the external loads in different situations. Then, the model is also applied to obtain the anisotropy and mechanical properties of the human proximal femur under physiological loads with initial conditions corresponding to a heterogeneous, but otherwise isotropic bone.     

A mechanostatistical approach to cortical bone remodelling: an equine model

In this study, the development of a mechanostatistical model of three-dimensional cortical bone remodelling informed with in vivo equine data is presented. The equine model was chosen as it is highly translational to the human condition due to similar Haversian systems, availability of in vivo bone strain and biomarker data, and furthermore, equine models are recommended by the US Federal Drugs Administration for comparative joint research. The model was derived from micro-computed tomography imaged specimens taken from the equine third metacarpal bone, and the Frost-based ‘mechanostat’ was informed from both in vivo strain gauges and biomarkers to estimate bone growth rates. The model also described the well-known ‘cutting cone’ phenomena where Haversian canals tunnel and replace bone. In order to make this model useful in practice, a partial least squares regression (PLSR) surrogate model was derived based on training data from finite element simulations with different loads. The PLSR model was able to predict microstructure and homogenised Young’s modulus with errors less than 2.2 % and \(0.6\,\% \), respectively.     

Effect of material damping on bone remodelling

This paper considers the effect of internal material damping on the stresses, strains, and surface and internal remodelling behaviour in a section of axisymmetrical bone with a force-fitted axially oriented medullary pin. The bone response to several loading situations is modelled using visco-elastic equations. An approximate method is developed to analyse the proposed mathematical model. By considering a numerical example, the effect of material damping on the remodelling stresses is quantified.     

A novel mathematical model of bone remodelling cycles for trabecular bone at the cellular level

After an initial phase of growth and development, bone undergoes a continuous cycle of repair, renewal and optimisation by a process called remodelling. This paper describes a novel mathematical model of the trabecular bone remodelling cycle. It is essentially formulated to simulate a remodelling event at a fixed position in the bone, integrating bone removal by osteoclasts and formation by osteoblasts. The model is developed to construct the variation in bone thickness at a particular point during the remodelling event, derived from standard bone histomorphometric analyses. The novelties of the approach are the adoption of a predator-prey model to describe the dynamic interaction between osteoclasts and osteoblasts, using a genetic algorithm-based solution; quantitative reconstruction of the bone remodelling cycle; and the introduction of a feedback mechanism in the bone formation activity to co-regulate bone thickness. The application of the model is first demonstrated by using experimental data recorded for normal (healthy) bone remodelling to predict the temporal variation in the number of osteoblasts and osteoclasts. The simulated histomorphometric data and remodelling cycle characteristics compare well with the specified input data. Sensitivity studies then reveal how variations in the model's parameters affect its output; it is hoped that these parameters can be linked to specific biochemical factors in the future. Two sample pathological conditions, hypothyroidism and primary hyperparathyroidism, are examined to demonstrate how the model could be applied more broadly, and, for the first time, the osteoblast and osteoclast populations are predicted for these conditions. Further data are required to fully validate the model's predictive capacity, but this work shows it has potential, especially in the modelling of pathological conditions and the optimisation of the treatment of those conditions.     

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.     

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.     

Inverse remodelling algorithm identifies habitual manual activities of primates based on metacarpal bone architecture

Previously, a micro-finite element (micro-FE)-based inverse remodelling method was presented in the literature that reconstructs the loading history of a bone based on its architecture alone. Despite promising preliminary results, it remains unclear whether this method is sensitive enough to detect differences of bone loading related to pathologies or habitual activities. The goal of this study was to test the sensitivity of the inverse remodelling method by predicting joint loading histories of metacarpal bones of species with similar anatomy but clearly distinct habitual hand use. Three groups of habitual hand use were defined using the most representative primate species: manipulation (human), suspensory locomotion (orangutan), and knuckle-walking locomotion (bonobo, chimpanzee, gorilla). Nine to ten micro-computed tomography scans of each species (\(n=48\) in total) were used to create micro-FE models of the metacarpal head region. The most probable joint loading history was predicted by optimally scaling six load cases representing joint postures ranging from \(-\,75^{\circ }\) (extension) to \(+\,75^{\circ }\) (flexion). Predicted mean joint load directions were significantly different between knuckle-walking and non-knuckle-walking groups (\(p<0.05\)) and in line with expected primary hand postures. Mean joint load magnitudes tended to be larger in species using their hands for locomotion compared to species using them for manipulation. In conclusion, this study shows that the micro-FE-based inverse remodelling method is sensitive enough to detect differences of joint loading related to habitual manual activities of primates and might, therefore, be useful for palaeoanthropologists to reconstruct the behaviour of extinct species and for biomedical applications such as detecting pathological joint loading.

     

Optimized intermolecular potential for nitriles based on Anisotropic United Atoms model

An extension of the anisotropic united atoms intermolecular potential model is proposed for nitriles. The electrostatic part of the intermolecular potential is calculated using atomic charges obtained by a simple Mulliken population analysis. The repulsion-dispersion interaction parameters for methyl and methylene groups are taken from transferable AUA4 literature parameters [Ungerer et al., J. Chem. Phys., 2000, 112, 5499]. Non-bonding Lennard-Jones intermolecular potential parameters are regressed for the carbon and nitrogen atoms of the nitrile group (–C≡N) from experimental vapor-liquid equilibrium data of acetonitrile. Gibbs Ensemble Monte Carlo simulations and experimental data agreement is very good for acetonitrile, and better than previous molecular potential proposed by Hloucha et al. [J. Chem. Phys., 2000, 113, 5401]. The transferability of the resulting potential is then successfully tested, without any further readjustment, to predict vapor-liquid phase equilibrium of propionitrile and n-butyronitrile. Figure Saturated vapour pressure of nitriles calculated in this work by molecular simulation compared to experimental data: a) for acetonitrile and b) for both propionitrile and butyronitrile