Time Dependent Behaviour of Trabecular Bone at Multiple Load Levels

How Artificial Intelligence and Machine Learning Is Assisting Us to Extract Meaning from Data on Bone Mechanics?

Dramatic advancements in interdisciplinary research with the fourth paradigm of science, especially the implementation of computer science, nourish the potential for artificial intelligence (AI), machine learning (ML), and artificial neural network (ANN) algorithms to be applied to studies concerning mechanics of bones. Despite recent enormous advancement in techniques, gaining deep knowledge to find correlations between bone shape, material, mechanical, and physical responses as well as properties is a daunting task. This is due to both complexity of the material itself and the convoluted shapes that this complex material forms. Moreover, many uncertainties and ambiguities exist concerning the use of traditional computational techniques that hinders gaining a full comprehension of this advanced biological material. This book chapter offers a review of literature on the use of AI, ML, and ANN in the study of bone mechanics research. A main question as to why to implement AI and ML in the mechanics of bones is fully addressed and explained. This chapter also introduces AI and ML and elaborates on the main features of ML algorithms such as learning paradigms, subtypes, main ideas with examples, performance metrics, training algorithms, and training datasets. As a frequently employed ML algorithm in bone mechanics, feedforward ANNs are discussed to make their taxonomy and working principles more readily comprehensible to researchers. A summary as well as detailed review of papers that employed ANNs to learn from collected data on bone mechanics are presented. Reviewing literature on the use of these data-driven tools is essential since their wider application has the potential to: improve clinical assessments enabling real-time simulations; avoid and/or minimize injuries; and, encourage early detection of such injuries in the first place.

Time-resolved in situ synchrotron-microCT: 4D deformation of bone and bone analogues using digital volume correlation

Digital volume correlation (DVC) in combination with high-resolution micro-computed tomography (microCT) imaging and in situ mechanical testing is gaining popularity for quantifying 3D full-field strains in bone and biomaterials. However, traditional in situ time-lapsed (i.e., interrupted) mechanical testing cannot fully capture the dynamic strain mechanisms in viscoelastic biological materials. The aim of this study was to investigate the time-resolved deformation of bone structures and analogues via continuous in situ synchrotron-radiation microCT (SR-microCT) compression and DVC to gain a better insight into their structure-function relationships. Fast SR-microCT imaging enabled the deformation behaviour to be captured with high temporal and spatial resolution. Time-resolved DVC highlighted the relationship between local strains and damage initiation and progression in the different biostructures undergoing plastic deformation, bending and/or buckling of their main microstructural elements. The results showed that SR-microCT continuous mechanical testing complemented and enhanced the information obtained from time-lapsed testing, which may underestimate the 3D strain magnitudes as a result of the stress relaxation occurring in between steps before image acquisition in porous biomaterials. Altogether, the findings of this study highlight the importance of time-resolved in situ experiments to fully characterise the time-dependent mechanical behaviour of biological tissues and biomaterials and to further explore their micromechanics under physiologically relevant conditions. STATEMENT OF SIGNIFICANCE: Time-resolved synchrotron X-ray tomography in combination with in situ mechanical testing provided the first four-dimensional analysis of the mechanical deformation of bone and bone analogues. To unravel the interplay of damage initiation and progression with local deformation, digital volume correlation was used to map the local strain field while microstructural changes were tracked with high temporal and spatial resolution. The results highlighted the importance of fast imaging and time-resolved in situ experiments to capture the real deformation of complex porous materials to fully characterize the local strain-damage relationship. The findings are notably improving the understanding of time-dependent mechanical behaviour of bone tissue, with the potential to be extend to highly viscoelastic biomaterials and soft tissues.

A predictive model for creep deformation following vertebral compression fractures

Many vertebral compression fractures continue to collapse over time, resulting in spinal deformity and chronic back pain. Currently, there is no adequate screening strategy to identify patients at risk of progressive vertebral collapse. This study developed a mathematical model to describe the quantitative relationship between initial bone damage and progressive ("creep") deformation in human vertebrae. The model uses creep rate before damage, and the degree of vertebral bone damage, to predict creep rate of a fractured vertebra following bone damage. Mechanical testing data were obtained from 27 vertebral trabeculae samples, and 38 motion segments, from 26 human spines. These were analysed to evaluate bone damage intensity, and creep rates before and after damage, in order to estimate the model parameter, p, which represents how bone damage affects the change of creep rate after damage. Results of the model showed that p was 1.38 (R² = 0.72, p < 0.001) for vertebral trabeculae, and 1.48 for motion segments (R² = 0.22, p = 0.003). These values were not significantly different from each other (P > 0.05). Further analyses revealed that p was not significantly influenced by cortical bone damage, endplate damage, disc degeneration, vertebral size, or vertebral areal bone mineral density (aBMD) (P > 0.05). The key determinant of creep deformation following vertebral compression fracture was the degree of trabecular bone damage. The proposed model could be used to identify the measures of bone damage on routine MR images that are associated with creep deformation so that a screening tool can be developed to predict progressive vertebral collapse following compression fracture.

A finite element analysis of tibial tritanium cones without stems in varying bone defects

BACKGROUND

In the UK around 10% of hip and knee arthroplasties are revision operations. At revision total knee arthroplasty (rTKA), bone loss management is critical to achieving a stable bone-implant construct. Though tritanium cones have been used to manage bone defects in rTKA, their biomechanical performance with varying defects remains unknown.

METHODS

Uncontained tibial bone defects at four anatomic locations, with varying depths and widths (Type T2A and T2B) were investigated computationally in a composite tibia which was subjected to four loading scenarios. The ability of the tritanium cone to replace the tibial bone defect was examined using the outcome measures of bone strain distribution and interface micromotions.

RESULTS

It was found that anterior and lateral defects do not significantly alter the strain distribution compared with intact bone. For medial defects, strain distribution is sensitive to defect width; while strain distributions for posterior defects are associated with defect width and depth. In general, micromotions at the bone-implant interface are small and are primarily influenced by defect depth.

CONCLUSIONS

Our models show that the cone is an acceptable choice for bone defect management in rTKA. Since all observed micromotions were small, successful osteointegration would be expected in all types of uncontained defects considered in this study. Tritanium cones safely accommodate uncontained tibial defects up to 10 mm deep and extending up to 9 mm from the centre of the cone. Medial and posteriorly based defects managed with symmetric cones display the greatest bone strains and asymmetric cones may be useful in this context.

The influence of long-term shoulder loading on sagittal spino-pelvic morphology: a population-based retrospective study of Chinese farmers from radiology

Background

To investigate associations between long-term shoulder loading and sagittal spino-pelvic morphology in Chinese farmers from radiology evidences.

Methods

We retrospectively analyzed 463 back pain patients who attended outpatient and inpatient departments of two hospitals from January 2016 to December 2018, and who had long, standing lateral X-rays according to inclusion and exclusion criteria. One hundred eighty-four of them were farmers with a long history of heavy shoulder loading for over 20 years in their young age, while others were office workers with no reported long-term shoulder loading history. The following parameters were measured by three researchers independently and then analyzed statistically: thoracic kyphosis (TK), lumbar lordosis (LL), thoracolumbar kyphosis (TLK), T9 sagittal offset (T9SO), T1 sagittal offset (T1SO), sacral slope (SS), pelvic incidence (PI), pelvic tilt (PT), C7 tilt (C7T), spino-pelvic angle (SSA), and sagittal vertical axis (SVA).

Results

The “Loading group” included 86 males and 98 females with average age 73.3 (SD 8.3) years, whereas the “Non-loading group” included 126 males and 153 females with average age 63.7 (SD 14.1) years. Age was significantly higher in the loading group (p < 0.001), but gender, height, weight, BMI, and BMD were not significantly different (p > 0.05). The following spino-pelvic parameters were significantly greater (p < 0.05) in the loading group: TK (mean 39.1° vs 32.8°), TLK (25.8° vs 10.1°), and T9SO (12.2° vs 10.1°). Other values were not significantly different between the two groups (p > 0.05).

Conclusion

Long-term shoulder loading in youth is a risk factor for pathological thoracic kyphosis especially in the lower thoracic spinal segments when farmers getting older.

Metaphyseal cones in revision total knee arthroplasty: The role of stems

Aims

Metaphyseal tritanium cones can be used to manage the tibial bone loss commonly encountered at revision total knee arthroplasty (rTKA). Tibial stems provide additional fixation and are generally used in combination with cones. The aim of this study was to examine the role of the stems in the overall stability of tibial implants when metaphyseal cones are used for rTKA.

Methods

This computational study investigates whether stems are required to augment metaphyseal cones at rTKA. Three cemented stem scenarios (no stem, 50 mm stem, and 100 mm stem) were investigated with 10 mm-deep uncontained posterior and medial tibial defects using four loading scenarios designed to mimic activities of daily living.

Results

Small micromotions (mean < 12 µm) were found to occur at the bone-implant interface for all loading cases with or without a stem. Stem inclusion was associated with lower micromotion, however these reductions were too small to have any clinical significance. Peak interface micromotion, even when the cone is used without a stem, was too small to effect osseointegration. The maximum difference occurred with stair descent loading. Stress concentrations in the bone occurred around the inferior aspect of each implant, with the largest occurring at the end of the long stem; these may lead to end-of-stem pain. Stem use is also found to result in stress shielding in the bone along the stem.

Conclusion

When a metaphyseal cone is used at rTKA to manage uncontained posterior or medial defects of up to 10 mm depth, stem use may not be necessary.

Cite this article:Bone Joint Res. 2020;9(4):162–172.

Time-dependent behaviour of demineralised trabecular bone - Experimental investigation and development of a constitutive model

Trabecular bone is a cellular composite material comprising primarily of mineral and organic phases and its mechanical response to loads is time-dependent. The contribution of the organic phase to the time-dependent behaviour of bone is not yet understood. We investigated the time-dependent response of demineralised trabecular bone through tensile multiple-load-creep-unload-recovery experiments. We found that demineralised trabecular bone's time-dependent response is nonlinearly related to the applied stress levels - it stiffens with increased stress levels. Our results also indicated that the time-dependent behaviour is associated with the original bone volume ratio (BV/TV). Irrecoverable strain exists, even at the low strain levels, but are not associated with BV/TV. Furthermore, we found that the nonlinear viscoelastic model can accurately predict the time-dependent behaviour of the trabecular bone's organic phase, which can be incorporated together with the properties of mineral to generate a composite model of bone. This study will help to provide a better understanding of this natural composite material.

Effect of loading frequency on deformations at the bone-implant interface

This study considers the time-dependent behaviour of bone in the context of loosening of metal implants, which is one of the typical complications of joint replacement and fracture-fixation surgeries. We employed viscoelastic properties developed from our previous experimental studies for trabecular bone in a representative bone-implant construct, which was subjected to cyclic loading at varying loading frequencies. We found that the separation between the bone and the implant is a function of loading frequency and increases with number of loading cycles applied. Our analysis shows that at the start of cyclic loading, a higher frequency results in a lower displacement response of bone at the bone-implant interface; however, after the bone-implant system has been subjected to a large number of cycles (>500 cycles in this study), higher interfacial displacements are observed at higher loading frequencies. In other words, higher loading frequencies will not result in bone-implant separation if limited number of cycles are applied. In all cases, interfacial displacements increase as bone volume ratio decreases. This simple approach can be used to evaluate the mechanical environment in bone-implant systems due to cyclic loading which commonly used time-independent models that are unable to simulate. The approach can also be used to evaluate implant loosening due to cyclic loading.

The risk of loosening of extramedullary fracture fixation devices

Extramedullary devices that use screws, pins or wires are used extensively to treat fractures in normal and diseased bone. A common failure mode is implant loosening at the bone-screw/pin/wire interface before fracture healing occurs. This review first considers the fundamental mechanics of the bone-fixator construct with focus on interfacial strains that result in loosening. It then evaluates the time-independent and time-dependent material models of bone that have been used to simulate and predict loosening. It is shown that the recently developed time-dependent models are capable of predicting loosening due to cyclic loads in bone of varying quality.

Time-dependent behaviour of bone accentuates loosening in the fixation of fractures using bone-screw systems

Aims

Loosening is a well-known complication in the fixation of fractures using devices such as locking plates or unilateral fixators. It is believed that high strains in the bone at the bone-screw interface can initiate loosening, which can result in infection, and further loosening. Here, we present a new theory of loosening of implants. The time-dependent response of bone subjected to loads results in interfacial deformations in the bone which accumulate with cyclical loading and thus accentuates loosening.

Methods

We used an 'ideal' bone-screw system, in which the screw is subjected to cyclical lateral loads and trabecular bone is modelled as non-linear viscoelastic and non-linear viscoelastic-viscoplastic material, based on recent experiments, which we conducted.

Results

We found that the interfacial deformation in the bone increases with the number of cycles, and the use of a non-linear viscoelastic-viscoplastic model results in larger deformations, some of which are irrecoverable. There is an apparent trend in which interfacial deformations increase with increasing porosity of bone.

Conclusion

The developed time-dependent model of the mechanical behaviour of bone permits prediction of loosening due to cyclical loads, which has not been possible previously. Application of this model shows that implant loosening will be accentuated by cyclical loading due to physiological activities, and the risks of loosening are greater in osteoporotic patients.Cite this article: S. Xie, K. Manda, P. Pankaj. Time-dependent behaviour of bone accentuates loosening in the fixation of fractures using bone-screw systems. Bone Joint Res 2018;7:580-586. DOI: 10.1302/2046-3758.710.BJR-2018-0085.R1.

Mechanical Stimuli in the Local In Vivo Environment in Bone: Computational Approaches Linking Organ-Scale Loads to Cellular Signals

PURPOSE OF REVIEW

Connecting organ-scale loads to cellular signals in their local in vivo environment is a current challenge in the field of bone (re)modelling. Understanding this critical missing link would greatly improve our ability to anticipate mechanotransduction during different modes of stimuli and the resultant cellular responses. This review characterises computational approaches that could enable coupling links across the multiple scales of bone.

RECENT FINDINGS

Current approaches using strain and fluid shear stress concepts have begun to link organ-scale loads to cellular signals; however, these approaches fail to capture localised micro-structural heterogeneities. Furthermore, models that incorporate downstream communication from osteocytes to osteoclasts, bone-lining cells and osteoblasts, will help improve the understanding of (re)modelling activities. Incorporating this potentially key information in the local in vivo environment will aid in developing multiscale models of mechanotransduction that can predict or help describe resultant biological events related to bone (re)modelling. Progress towards multiscale determination of the cell mechanical environment from organ-scale loads remains elusive. Construction of organ-, tissue- and cell-scale computational models that include localised environmental variation, strain amplification and intercellular communication mechanisms will ultimately help couple the hierarchal levels of bone.

From Tension to Compression: Asymmetric Mechanical Behaviour of Trabecular Bone's Organic Phase

Trabecular bone is a cellular composite material comprising primarily of mineral and organic phases with their content ratio known to change with age. Therefore, the contribution of bone constituents on bone's mechanical behaviour, in tension and compression, at varying load levels and with changing porosity (which increases with age) is of great interest, but remains unknown. We investigated the mechanical response of demineralised bone by subjecting a set of bone samples to fully reversed cyclic tension-compression loads with varying magnitudes. We show that the tension to compression response of the organic phase of trabecular bone is asymmetric; it stiffens in tension and undergoes stiffness reduction in compression. Our results indicate that demineralised trabecular bone struts experience inelastic buckling under compression which causes irreversible damage, while irreversible strains due to microcracking are less visible in tension. We also identified that the values of this asymmetric mechanical response is associated to the original bone volume ratio (BV/TV).

The relationship of whole human vertebral body creep to geometric, microstructural, and material properties

Creep, the time dependent deformation of a structure under load, is an important viscoelastic property of bone and may play a role in the development of permanent deformity of the vertebrae in vivo leading to clinically observable spinal fractures. To date, creep properties and their relationship to geometric, microstructural, and material properties have not been described in isolated human vertebral bodies. In this study, a range of image-based measures of vertebral bone geometry, bone mass, microarchitecture and mineralization were examined in multiple regression models in an effort to understand their contribution to creep behavior. Several variables, such as measures of mineralization heterogeneity, average bone density, and connectivity density persistently appeared as significant effects in multiple regression models (adjusted r²: 0.17-0.56). Although further work is needed to identify additional tissue properties to fully describe the portion of variability not explained by these models, these data are expected to help understand mechanisms underlying creep and improve prediction of vertebral deformities that eventually progress to a clinically observable fracture.

Quantification of interfacial motions following primary and revision total knee arthroplasty: A verification study versus experimental data

Motion at the bone-implant interface, following primary or revision knee arthroplasty, can be detrimental to the long-term survival of the implant. This study employs experimentally verified computational models of the distal femur to characterize the relative motion at the bone-implant interface for three different implant types; a posterior stabilizing implant (PS), a total stabilizing implant (TS) with short stem (12 mm × 50 mm), and a total stabilizing implant (TS) with long offset stem (19 mm × 150 mm with a 4 mm lateral offset). Relative motion was investigated for both cemented and uncemented interface conditions. Monitoring relative motion about a single reference point, though useful for discerning global differences between implant types, was found to not be representative of the true pattern and distribution of motions which occur at the interface. The contribution of elastic deformation to apparent reference point motion varied based on implant type, with the PS and TSSS implanted femurs experiencing larger deformations (43 and 39 μm, respectively) than the TSLS implanted femur (22 μm). Furthermore, the pattern of applied loading was observed to greatly influence location and magnitude of peak motions, as well as the surface area under increased motion. Interestingly, the influence was not uniform across all implant types, with motions at the interface of long stemmed prosthesis found to be less susceptible to changes in pattern of loading. These findings have important implications for the optimization and testing of orthopedic implants in vitro and in silico. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:387-396, 2018.