A Bone-remodelling Scheme Based on Principal Strains Applied to a Tooth During Translation

This paper investigates the role of principal strains within the periodontal ligament (PDL) during bone remodelling in orthodontics and particularly in the case of bodily motion (pure translation). Using analytical formulas of stress and strains within the PDL for the particular case of a paraboloidal central incisor during translation, the strains are directly related to the motion of the interface between the alveolar bone and the PDL, called bone surface. It is shown that both normal and shear strains within the PDL are of the same importance for bone surface motion. Moreover, both “mean average” and “geometrical average” of principal strains within the PDL play a significant role in the bone remodelling process, as they contribute with the same proportionality. In summary, the proposed formulas differ than previous ones that had been successfully applied to describe remodelling within long bones. The proposed theory is also sustained by a linear finite element analysis.     

Numerical study of tension and strain distribution around rat molars]

A knowledge of the mechanical processes triggered in the bone and periodontal ligament (PDL) by orthodontic forces applied to a tooth is of decisive importance for an understanding of the subsequent remodelling around the tooth. To investigate these mechanical relationships, three-dimensional finite element (FE) models of the first lower molar in the rat were established. On the basis of digitized serial histological sections, these FE models were generated semi-automatically. Using various simplified geometrical variations, an appropriate FE model for the analysis of the stress and strain distributions was established. The numerical analyses were carried out under a mesially directed force of 0.1 N. Stress distributions in the bone and PDL showed a similar pattern, while strains in the bone were lower than in the PDL by a factor of 10-5. The data confirm the assumption that strain patterns in the PDL may be the key stimulus of bone remodelling.     

The importance of the elastic and plastic components of strain in tensile and compressive fatigue of human cortical bone in relation to orthopaedic biomechanics

The longevity, success, or failure of an orthopaedic implant is dependent on its osseointegration especially within the initial six months of the initial surgery. The development of strains plays a crucial role in both bone modelling and remodelling. For remodelling, in particular, strains of substantial values are required to activate the osteoblastic and osteoclastic activity for the osseointegration of the implant. Bone, however, is subject to "damage" when strain levels exceed a certain threshold level. Damage is manifested in the form of microcracks; it is linked to increased elastic strain amplitudes and is accompanied by the development of "plastic" (irrecoverable, residual) strains. Such strains increase the likelihood for the implant to subside or loosen. The present study examines the rates (per cycle) by which these two components of strain (elastic and "plastic") develop during fatigue cycling in two loading modes, tension and compression. The results of this study show that these strain rates depend on the applied stress in both loading modes. It also shows that elastic and plastic strain rates can be linked to each other through simple power law relationships so that one can calculate or predict the latter from the former and vice versa. We anticipate that such basic bone biomechanics data would be of great benefit to both clinicians and bioengineers working in the field of FEA modelling applications and orthopaedic implant surgery.     

Flexible multibody simulation approach in the analysis of tibial strain during walking

Strains within the bone tissue play a major role in bone (re)modeling. These small strains can be assessed using experimental strain gage measurements, which are challenging and invasive. Further, the strain measurements are, in practise, limited to certain regions of superficial bones only, such as the anterior surface of the tibia. In this study, tibial strains occurring during walking were estimated using a numerical approach based on flexible multibody dynamics. In the introduced approach, a lower body musculoskeletal model was developed by employing motion capture data obtained from walking at a constant velocity. The motion capture data were used in inverse dynamics simulation to teach the muscles in the model to replicate the motion in forward dynamics simulation. The maximum and minimum tibial principal strains predicted by the model were 490 and -588 microstrain, respectively, which are in line with literature values from in vivo measurements. In conclusion, the non-invasive flexible multibody simulation approach may be used as a surrogate for experimental bone strain measurements and thus be of use in detailed strain estimations of bones in different applications.     

A visco-elastic model for the prediction of orthodontic tooth movement

This study presents a biomechanical model of orthodontic tooth movement. Although such models have already been presented in the literature, most of them incorporate computationally expensive finite elements (FE) methods to determine the strain distribution in the periodontal ligament (PDL). In contrast, the biomechanical model presented in this work avoids the use of FE methods. The elastic deformation of the PDL is modelled using an analytical approach, which does not require setting up a 3D model of the tooth. The duration of the lag phase is estimated using the calculated hydrostatic stresses, and bone remodelling is predicted by modelling the alveolar bone as a viscous material. To evaluate the model, some typically used motion patterns were simulated and a sensitivity analysis was carried out on the parameters. Results show that despite some shortcomings, the model is able to describe commonly used motion patterns in orthodontic tooth movement, in both single- and multi-rooted teeth.     

Experimental model of tooth mobility in the human "in vivo"]

The aim of the present study was to investigate experimentally the mechanical properties of tooth deflection under external loading. These properties have a significant impact on tooth movement during orthodontic treatment. The stresses and strains caused by tooth movement influence bone remodelling, which is the basis of orthodontic treatment. The movement of a tooth as a direct reaction to the forces acting on it is termed "initial" movement. It is nonlinear and has a clearly time-dependent component. While the initial tooth movement represents the totality of the reaction mechanisms of all the tissues of the tooth unit, it is determined primarily by the mechanical properties of the periodontal ligament (PDL). The PDL is the softest tissue of the tooth unit and is therefore subject to the largest deformations when forces act on the crown of the tooth. The objective of orthodontic treatment is to achieve as precise and rapid tooth movement as possible, without provoking such undesired effects as bone and root resorption. To enable the implementation of an optimal orthodontic force system that meets these requirements, a thorough knowledge of the biomechanics of tooth movement is a must.     

The Biomechanical Function of Periodontal Ligament Fibres in Orthodontic Tooth Movement

Orthodontic tooth movement occurs as a result of resorption and formation of the alveolar bone due to an applied load, but the stimulus responsible for triggering orthodontic tooth movement remains the subject of debate. It has been suggested that the periodontal ligament (PDL) plays a key role. However, the mechanical function of the PDL in orthodontic tooth movement is not well understood as most mechanical models of the PDL to date have ignored the fibrous structure of the PDL. In this study we use finite element (FE) analysis to investigate the strains in the alveolar bone due to occlusal and orthodontic loads when PDL is modelled as a fibrous structure as compared to modelling PDL as a layer of solid material. The results show that the tension-only nature of the fibres essentially suspends the tooth in the tooth socket and their inclusion in FE models makes a significant difference to both the magnitude and distribution of strains produced in the surrounding bone. The results indicate that the PDL fibres have a very important role in load transfer between the teeth and alveolar bone and should be considered in FE studies investigating the biomechanics of orthodontic tooth movement.     

Anisotropic bone remodelling model based on a continuum damage-repair theory.

The purpose of this paper is to present and discuss a new model for internal bone remodelling. We state 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 exact principles of Continuum Damage Mechanics (CDM), although no actual damage is considered. A remodelling tensor, analogous to the standard damage tensor, is proposed which completely characterises the state of the homogenised bone microstructure and, therefore, its stiffness. This tensor is defined in terms of the apparent density and the "fabric tensor" associated with porosity and directionality of the trabeculae, respectively. Contrary to standard damage mechanics, its variation may be negative to allow for material "repair". The different elements that compose the theory are then established, such as the mechanical stimulus which, as is common in CDM, was chosen as the tensor that is thermodynamically associated to the remodelling one. The resorption and apposition criteria (similar to the damage criterion) were also expressed in terms of this mechanical stimulus and, finally, the evolution law, to define the change rate of the remodelling tensor. An associated flow rule is proposed that fulfils a principle of minimum mechanical dissipation for convex remodelling criteria (the case in question here). Some other important experimental features are deduced naturally from the model, such as the coincidence of the principal directions of the fabric tensor with those of the elasticity tensor, or the fact that the principal directions of the fabric tensor tend to align with the principal directions of stress. This model is applied to the study of the remodelling evolution of the proximal extremity of the femur, obtaining results very similar to experimental data.     

Response of the tooth-periodontal ligament-bone complex to load: A microCT study of the minipig molar

One strategy evolved by teeth to avoid irreversible damage is to move and deform under the loads incurred during mastication. A key component in this regard is the periodontal ligament (PDL). The role of the bone underlying the PDL is less well defined. We study the interplay between the PDL and the underlying alveolar bone when loaded in the minipig. Using an Instron loading device we confirmed that the force-displacement curves of the molars and premolars of relatively fresh minipig intact mandibles are similar to those obtained for humans and other animals. We then used this information to obtain 3D images of the teeth before and after loading the tooth in a microCT such that the load applied is in the third linear part of the force displacement curve. We observed that at many locations there is a complimentary topography of the cementum and alveolar bone surface, strongly suggesting an active interplay between the tooth and the bone during mastication. We also observed that the loaded tooth does not come into direct contact with the underlying bone surface. A highly compressed layer of PDL is present between the tooth and the bone. The structure of the bone in the upper furcation region has a unique appearance with little obvious microstructure, abundant pores that have a large size range and at many locations the bone at the PDL interface has a needle-like shape. We conclude that there is a close interaction between the tooth, the PDL and the underlying alveolar bone during mastication. The highly compressed PDL layer that separates the tooth from the bone may fulfill a key shock absorbing function.     

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.     

The effects of the periodontal ligament on mandibular stiffness: a study combining finite element analysis and geometric morphometrics

It is generally accepted that the periodontal ligament (PDL) plays a crucial role in transferring occlusal forces from the teeth to the alveolar bone. Studies using finite element analysis (FEA) have helped to better understand this role and show that the stresses and strains in the alveolar bone are influenced by whether and how PDL is included in FE models. However, when the overall distribution of stresses and strains in crania and mandibles are of interest, PDL is often not included in FE models, although little is known about how this affects the results. Here we study the effect of representing PDL as a layer of solid material with isotropic homogeneous properties in an FE model of a human mandible using a novel application of geometric morphometrics. The results show that the modelling of the PDL affects the deformation and thus strain magnitudes not only of the alveolar bone around the biting tooth, but that the whole mandible deforms differently under load. As a result, the strain in the mandibular corpus is significantly increased when PDL is included, while the strain in the bone beneath the biting tooth is reduced. These results indicate the importance of the PDL in FE studies. Thus we recommend that the PDL should be included in FE models of the masticatory apparatus, with tests to assess the sensitivity of the results to changes in the Young's modulus of the PDL material.     

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).     

Kinematics of cytoplasmic deformation in neutrophils during active motion

A procedure is proposed to measure the cytoplasmic deformation in active motile neutrophils in the form of cytoplasmic strains and strain rates. Three neighboring microspheres in a local region of the cytoplasm serve as markers for local motion. Their positions are tracked by means of a high resolution light microscope and serve to compute nonlinear measures of strains and strain rates together with the principal strains and principal directions. Active neutrophils exhibit large cytoplasmic strains both during periodic pseudopod projections and during continuous locomotion in a polarized shape. The cytoplasmic motion is often synchronized with the whole cell deformation. The local cytoplasmic strains exceed the strains estimated for the whole cell and are not reversible except in some cases of single pseudopod projections. Large strains are observed both in attached and freely suspended cells. Strain rates are relatively constant but show an increase during the pseudopod retraction phase. Local cytoplasmic strains in neutrophils are inhomogeneous and reach large values during passage of the contraction rings. Neutrophils rendered passive by treatment with cytochalasin or EDTA show a random motion of microspheres with much smaller displacements. These observations suggest that the cytoplasm of active neutrophil exhibits large cytoplasmic strains and strain rates in the absence of an external stress resulting in a high degree of intracellular mixing. The proposed technique may be applied to a wide range of problems in cell biology.     

Determination of dynamically adapting anisotropic material properties of bone under cyclic loading

Because bone tissue adapts to loading conditions, finite element simulations of remodelling bone require a precise prediction of dynamically changing anisotropic elastic parameters. We present a phenomenological theory that refers to the tissue in terms of the tendency of the structure to align with principal stress directions. We describe the material parameters of remodelling bone. This work follows findings by the same research group and independently by Danilov (1971) in the field of plasticity, where the dependencies of the components of the stiffness tensor in terms of time are based on Hill's anisotropy. We modify such an approach in this novel theory that addresses bone tissue that can regenerate. The computational assumption of the theory is that bone trabeculae have the tendency to orient along one of the principal stress directions but during remodelling the principal stresses change continuously and the resulting orientation of the trabeculae can differ from the principal stress direction at any given time. The novelty of this work consists in the limited number of parameters needed to compute the twenty-one anisotropic material parameters at any given location in the bone tissue. In addition to the theory, we present here two cases of simplified geometry, loading and boundary conditions to show the effect of (1) time on the material properties; and (2) change of loading conditions on the anisotropic parameters. The long term goal is to experimentally verify that the predictions generated by theory provide a reliable simulation of cancellous bone properties.     

Trabecular bone fracture healing simulation with finite element analysis and fuzzy logic

Trabecular bone fractures heal through intramembraneous ossification. This process differs from diaphyseal fracture healing in that the trabecular marrow provides a rich vascular supply to the healing bone, there is very little callus formation, woven bone forms directly without a cartilage intermediary, and the woven bone is remodelled to form trabecular bone. Previous studies have used numerical methods to simulate diaphyseal fracture healing or bone remodelling, however not trabecular fracture healing, which involves both tissue differentiation and trabecular formation. The objective of this study was to determine if intramembraneous bone formation and remodelling during trabecular bone fracture healing could be simulated using the same mechanobiological principles as those proposed for diaphyseal fracture healing. Using finite element analysis and the fuzzy logic for diaphyseal healing, the model simulated formation of woven bone in the fracture gap and subsequent remodelling of the bone to form trabecular bone. We also demonstrated that the trabecular structure is dependent on the applied loading conditions. A single model that can simulate bone healing and remodelling may prove to be a useful tool in predicting musculoskeletal tissue differentiation in different vascular and mechanical environments.     

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.     

Functional strain in bone tissue as an objective, and controlling stimulus for adaptive bone remodelling

The skeleton consists of a series of elements with a variety of functions. In locations where shape or protection are of prime importance the bone's architecture is achieved during growth under direct genetic control. In locations where resistance to repetitive loading is important only the general form of the bone will be achieved as a result of growth alone, the remaining characteristics result from functional adaptation. This mechanism ensures that bone architecture is modelled and remodelled until prevailing strains match those genetically prescribed for that location. For this match to be established, and subsequently maintained, bone cells must be able to 'assess' feedback derived directly or indirectly from the functional strains produced within the tissue. These strains are therefore the objective of functionally adaptive remodelling, and the stimulus for its control. Evans was the first person to refer to the recording of functional strains from gauges attached to bone in vivo. This technique has allowed quantitative investigations on bone's normal functional strain environment, and its adaptive response to changes in its state of strain. Recent investigations have extended to the immediate effects of dynamic strains on the structure of the bone matrix, and the biochemical behaviour of the resident bone cells. Such studies should reveal the mechanism by which strains within the matrix are transduced into the biochemical signals by which adaptive remodelling is controlled.     

Bone remodelling algorithms incorporating both strain and microdamage stimuli

Biomechanical theories to predict bone remodelling have used either mechanical strain or microdamage as the stimulus driving cellular responses. Even though experimental data have implicated both stimuli in bone cell regulation, a mechano-regulatory system incorporating both stimuli has not yet been proposed. In this paper, we test the hypothesis that bone remodelling may be regulated by signals due to both strain and microdamage. Four mechano-regulation algorithms are studied where the stimulus is: strain, damage, combined strain/damage, and either strain or damage with damage-adaptive remodelling prioritised when damage is above a critical level. Each algorithm is implemented with both bone lining cell (surface) sensors and osteocyte cell (internal) sensors. Each algorithm is applied to prediction of a bone multicellular unit (BMU) remodelling on the surface of a bone trabecula. It is predicted that a regulatory system capable of responding to changes in either strain or microdamage but which prioritises removal of damaged bone when damage is above a critical level, is the only one that provides a plausible prediction of BMU behaviour. A mechanism for this may be that, below a certain damage threshold, osteocyte processes can sense changes in strain and fluid flow but above the threshold damage interferes with the signalling mechanism, or causes osteocyte apoptosis so that a remodelling response occurs to remove the dead osteocytes.     

Homogenization theory and digital imaging: A basis for studying the mechanics and design principles of bone tissue

Bone tissue is a complex multilevel composite which has the ability to sense ad respond to its mechanical environment. It is believed that bone cells called osteocytes within the bone matrix sense the mechanical environment and determine whether structural alterations are needed. At present it is not known, however, how loads are transferred from the whole bone level to cells. A computational procedure combining representative volume element (RVE) based homogenization theory with digital imaging is proposed to estimate strains at various levels of bone structure. Bone tissue structural organization and RVE based analysis are briefly reviewed. The digital image based computational procedure was applied to estimate strains in individual trabeculae (first-level microstructure). Homogenization analysis of an idealized model was used to estimate strains at one level of bone structure around osteocyte lacunae (second-level trabecular microstructure). The results showed that strain at one level of bone structure is amplified to a broad range at the next microstructural level. In one case, a zeor-level tensile principal strain of 495 muE engendered strains ranging between -1000 and 7000 muE in individual trabeculae (first-level microstructure). Subsequently, a first-level tensile principal strains of 1325 muE within an inidividual trabecula engendered strains ranging between 782 and 2530 muE around osteocyte lacunae. Lacunar orientation was found to influence strains around osteocyte lacunae much more than lacunar ellipticity. In conclusion, the computational procedure combining homogenization theory with digital imaging can proveide estimates of cell level strains within whole bones. Such results may be used to bridge experimental studies of bone adaptation at the whole bone and cell culture level. (c) 1994 John Wiley & Sons, Inc.     

A musculoskeletal model of the equine forelimb for determining surface stresses and strains in the humerus-part II. Experimental testing and model validation

The first objective of this study was to experimentally determine surface bone strain magnitudes and directions at the donor site for bone grafts, the site predisposed to stress fracture, the medial and cranial aspects of the transverse cross section corresponding to the stress fracture site, and the middle of the diaphysis of the humerus of a simplified in vitro laboratory preparation. The second objective was to determine whether computing strains solely in the direction of the longitudinal axis of the humerus in the mathematical model was inherently limited by comparing the strains measured along the longitudinal axis of the bone to the principal strain magnitudes and directions. The final objective was to determine whether the mathematical model formulated in Part I [Pollock et al., 2008, ASME J. Biomech. Eng., 130, p. 041006] is valid for determining the bone surface strains at the various locations on the humerus where experimentally measured longitudinal strains are comparable to principal strains. Triple rosette strain gauges were applied at four locations circumferentially on each of two cross sections of interest using a simplified in vitro laboratory preparation. The muscles included the biceps brachii muscle in addition to loaded shoulder muscles that were predicted active by the mathematical model. Strains from the middle grid of each rosette, aligned along the longitudinal axis of the humerus, were compared with calculated principal strain magnitudes and directions. The results indicated that calculating strains solely in the direction of the longitudinal axis is appropriate at six of eight locations. At the cranial and medial aspects of the middle of the diaphysis, the average minimum principal strain was not comparable to the average experimental longitudinal strain. Further analysis at the remaining six locations indicated that the mathematical model formulated in Part I predicts strains within +/-2 standard deviations of experimental strains at four of these locations and predicts negligible strains at the remaining two locations, which is consistent with experimental strains. Experimentally determined longitudinal strains at the middle of the diaphysis of the humerus indicate that tensile strains occur at the cranial aspect and compressive strains occur at the caudal aspect while the horse is standing, which is useful for fracture fixation.