Abstract:The microstructure of trabecular bone is known to adapt its morphology in response to mechanical loads for achieving a biomechanical homeostasis. Based on this form–function relationship, previous investigators either simulated the remodeling of bone to predict the resulting density and architecture for a specific loading or retraced physiological loading conditions from local density and architecture. The latter inverse approach includes quantifying bone morphology using computed tomography and calculating the relative importance of selected load cases by minimizing the fluctuation of a tissue loading level metric. Along this concept, the present study aims at identifying an optimal, personalized, multiaxial load case at the distal section of the human radius using in vivo HR-pQCT-based isotropic, homogenized finite element (hFE) analysis. The dataset consisted of HR-pQCT reconstructions of the 20 mm most distal section of 21 human fresh-frozen radii. We simulated six different unit canonical load cases (FX palmar–dorsal force, FY ulnar–radial force, FZ distal–proximal force, MX moment about palmar–dorsal, MY moment about ulnar–radial, MZ moment about distal–proximal) using a simplified and efficient hFE method based on a single isotropic bone phase. Once we used a homogeneous mean density (shape model) and once the original heterogeneous density distribution (shape + density model). Using an analytical formulation, we minimized the deviation of the resulting strain tensors ε ( x ) to a hydrostatic compressive reference strain ε 0 , once for the 6 degrees of freedom (DOF) optimal (OPT) load case and for all individual 1 DOF load cases (FX, FY, FZ, MX, MY, MZ). All seven load cases were then extended in the nonlinear regime using the scaled displacements of the linear load cases as loading boundary conditions (MAX). We then compared the load cases and models for their objective function (OF) values, the stored energies and their ultimate strength using a specific torsor norm. Both shape and shape + density linear-optimized OPT models were dominated by a positive force in the z -direction (FZ). Transversal force DOFs were close to zero and mean moment DOFs were different depending on the model type. The inclusion of density distribution increased the influence and changed direction of MX and MY, while MZ was small in both models. The OPT load case had 12–15% lower objective function (OF) values than the FZ load case, depending on the model. Stored energies at the optimum were consistently 142–178% higher for the OPT load case than for the FZ load case. Differences in the nonlinear response maximum torsor norm ‖t‖ were heterogeneous, but consistently higher for OPT_MAX than FZ_MAX. We presented the proof of concept of an optimization procedure to estimate patient-specific loading conditions for hFE methods. In contrast to similar models, we included canonical load cases in all six DOFs and used a strain metric that favors hydrostatic compression. Based on a biomechanical analysis of the distal joint surfaces at the radius, the estimated load directions are plausible. For our dataset, the resulting OPT load case is close to the standard axial compression boundary conditions, usually used in HR-pQCT-based FE analysis today. But even using the present simplified hFE model, the optimized linear six DOF load case achieves a more homogeneous tissue loading and can absorb more than twice the energy than the standard uniaxial load case. The ultimate strength calculated with a torsor norm was consistently higher for the 6-DOF nonlinear model (OPT_MAX) than for the 1-DOF nonlinear uniaxial model (FZ_MAX). Defining patient-specific boundary conditions may decrease angulation errors during CT measurements and improve repeatability as well as reproducibility of bone stiffness and strength estimated by HR-pQCT-based hFE analysis. These results encourage the extension of the present method to anisotropic hFE models and their application to repeatability data sets to test the hypothesis of reduced angulation errors during measurement.

Exploring conditions that make cortical bone geometry optimal for physiological loading

Three-dimensional micro-level computational study of Wolff's law via trabecular bone remodeling in the human proximal femur using design space topology optimization

Topology Optimization of Continuum Structures Based on a New Bionics Method

Optimization based synthesis of pelvic structure for loads in running gait cycle

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

Design of complex bone internal structure using topology optimization with perimeter control

Do bone ingrowth processes produce a globally optimized structure?

Personalized loading conditions for homogenized finite element analysis of the distal sections of the radius

Biomimetic design of implants for long bone critical-sized defects

Skeletal structural adaptations to mechanical usage (SATMU): 1. Redefining Wolff's Law: The bone modeling problem

Computer modelling of bone’s adaptation: the role of normal strain, shear strain and fluid flow

Computationally derived endosteal strain and strain gradients correlate with increased bone formation in an axially loaded murine tibia model

Influence of interstitial fluid pressure, porosity, loading magnitude, and anisotropy in cortical bone adaptation

A Mathematical Model of Cortical Bone Remodeling at Cellular Level under Mechanical Stimulus

Effects of mechanical forces on maintenance and adaptation of form in trabecular bone

Cortical Thickness Adaptive Response to Mechanical Loading Depends on Periosteal Position and Varies Linearly With Loading Magnitude

Differences in the mechanical behavior of cortical bone between compression and tension when subjected to progressive loading

Role of Dynamic Loading on Early Stage of Bone Fracture Healing

The Constructure Alteration Simulation of A Type of Self-growth Material of Whistar Rat

Spatiotemporal Characterization of Microstructure Morphology, Mechanical Properties and Bone Remodeling of Rat Tibia Under Uniaxial Compressive Overload Loading

Application research of a structural topology optimization method based on a bionic principle