Abstract:Bone is an extraordinary biological material that continuously adapts its hierarchical microstructure to respond to static and dynamic loads for offering optimal mechanical features, in terms of stiffness and toughness, across different scales, from the sub-microscopic constituents within osteons-where the cyclic activity of osteoblasts, osteoclasts, and osteocytes redesigns shape and percentage of mineral crystals and collagen fibers-up to the macroscopic level, with growth and remodeling processes that modify the architecture of both compact and porous bone districts. Despite the intrinsic complexity of the bone mechanobiology, involving coupling phenomena of micro-damage, nutrients supply driven by fluid flowing throughout hierarchical networks, and cells turnover, successful models and numerical algorithms have been presented in the literature to predict, at the macroscale, how bone remodels under mechanical stimuli, a fundamental issue in many medical applications such as optimization of femur prostheses and diagnosis of the risk fracture. Within this framework, one of the most classical strategies employed in the studies is the so-called Stanford's law, which allows uploading the effect of the time-dependent load-induced stress stimulus into a biomechanical model to guess the bone structure evolution. In the present work, we generalize this approach by introducing the bone poroelasticity, thus incorporating in the model the role of the fluid content that, by driving nutrients and contributing to the removal of wastes of bone tissue cells, synergistically interacts with the classical stress fields to change homeostasis states, local saturation conditions, and reorients the bone density rate, in this way affecting growth and remodeling. Through two paradigmatic example applications, i.e. a cylindrical slice with internal prescribed displacements idealizing a tract of femoral diaphysis pushed out by the pressure exerted by a femur prosthesis and a bone element in a form of a bent beam, it is highlighted that the present model is capable to catch more realistically both the transition between spongy and cortical regions and the expected non-symmetrical evolution of bone tissue density in the medium-long term, unpredictable with the standard approach. A real study case of a femur is also considered at the end in order to show the effectiveness of the proposed remodeling algorithm.

Using Nonlinear Response to Estimate the Strength of an Elastic Network

A Model for Bone Strength and Osteoporotic Fractures

Fracture prediction for the proximal femur using finite element models: Part II--Nonlinear analysis

Network models for characterization of trabecular bone

Response of Biopolymer Networks Governed by the Physical Properties of Cross-Linking Molecules.

The sensitivity of nonlinear computational models of trabecular bone to tissue level constitutive model

A proposal for the combined analysis of bone quantity and quality of human cortical bone by quasi-brittle fracture mechanics

Neural network prediction of load from the morphology of trabecular bone

Inelastic Modelling of Bone Damage under Compressive Loading

Stochastic failure analysis of proximal femur using an isogeometric analysis based nonlocal gradient-enhanced damage model

Constitutive relationships for osteonal microcracking in human cortical bone using statistical mechanics

A new constitutive relation to describe the response of bones

Acoustic emission applied to stochastic modeling of microdamage in compact bone

Effect of including damage at the tissue level in the nonlinear homogenisation of trabecular bone

Patient-specific fracture risk assessment of vertebrae: A multiscale approach coupling X-ray physics and continuum micromechanics

Symmetry breaking and effects of nutrient walkway in time-dependent bone remodeling incorporating poroelasticity

Predicting the mechanical properties of spring networks

Estimation of the in-plane ultimate stress of lamellar tissue as a function of bone mineral density and osteocyte lacunae porosity

Evaluation of nonlocal approaches for modelling fracture near nonconvex boundaries

A method for patient-specific evaluation of vertebral cancellous bone strength: in vitro validation

Experimental validation of finite element analysis of human vertebral collapse under large compressive strains