Computational study of Wolff's law with trabecular architecture in the human proximal femur using topology optimization

In the field of bone adaptation, it is believed that the morphology of bone is affected by its mechanical loads, and bone has self-optimizing capability; this phenomenon is well known as Wolff's law of the transformation of bone. In this paper, we simulated trabecular bone adaptation in the human proximal femur using topology optimization and quantitatively investigated the validity of Wolff's law. Topology optimization iteratively distributes material in a design domain producing optimal layout or configuration, and it has been widely and successfully used in many engineering fields. We used a two-dimensional micro-FE model with 50 microm pixel resolution to represent the full trabecular architecture in the proximal femur, and performed topology optimization to study the trabecular morphological changes under three loading cases in daily activities. The simulation results were compared to the actual trabecular architecture in previous experimental studies. We discovered that there are strong similarities in trabecular patterns between the computational results and observed data in the literature. The results showed that the strain energy distribution of the trabecular architecture became more uniform during the optimization; from the viewpoint of structural topology optimization, this bone morphology may be considered as an optimal structure. We also showed that the non-orthogonal intersections were constructed to support daily activity loadings in the sense of optimization, as opposed to Wolff's drawing.     

Computer simulation of trabecular remodeling in human proximal femur using large-scale voxel FE models: Approach to understanding Wolff's law

Ever since Julius Wolff proposed the law of bone transformation in the 19th century, it has been widely known that the trabecular structure of cancellous bone adapts functionally to the loading environment. To understand the mechanism of Wolff's law, a three-dimensional (3D) computer simulation of trabecular structural changes due to surface remodeling was performed for a human proximal femur. A large-scale voxel finite element model was constructed to simulate the structural changes of individual trabeculae over the entire cancellous region. As a simple remodeling model that considers bone cellular activities regulated by the local mechanical environment, nonuniformity of local stress was assumed to drive the trabecular surface remodeling to seek a uniform stress state. Simulation results demonstrated that cell-scale (～10 μm) remodeling in response to mechanical stimulation created complex 3D trabecular structures of the entire bone-scale (～10 cm), as illustrated in the reference of Wolff. The bone remodeling reproduced the characteristic anisotropic structure in the coronal cross section and the isotropic structures in other cross sections. The principal values and axes of a structure characterized by fabric ellipsoids corresponded to those of the apparent stress of the structure. The proposed large-scale computer simulation indicates that in a complex mechanical environment of a hierarchical bone structure of over 10⁴ length scale (from ～10 μm to ～10 cm), a simple remodeling at the cellular/trabecular levels creates a highly complex and functional trabecular structure, as characterized by bone density and orientation.     

Comparison of different hip prosthesis shapes considering micro-level bone remodeling and stress-shielding criteria using three-dimensional design space topology optimization

Since the late 1980s, computational analysis of total hip arthroplasty (THA) prosthesis components has been completed using macro-level bone remodeling algorithms. The utilization of macro-sized elements requires apparent bone densities to predict cancellous bone strength, thereby, preventing visualization and analysis of realistic trabecular architecture. In this study, we utilized a recently developed structural optimization algorithm, design space optimization (DSO), to perform a micro-level three-dimensional finite element bone remodeling simulation on the human proximal femur pre- and post-THA. The computational simulation facilitated direct performance comparison between two commercially available prosthetic implant stems from Zimmer Inc.: the Alloclassic and the Mayo conservative. The novel micro-level approach allowed the unique ability to visualize the trabecular bone adaption post-operation and to quantify the changes in bone mineral content by region. Stress-shielding and strain energy distribution were also quantified for the immediate post-operation and the stably fixated, post-remodeling conditions. Stress-shielding was highest in the proximal region and remained unchanged post-remodeling; conversely, the mid and distal portions show large increases in stress, suggesting a distal shift in the loadpath. The Mayo design conserves bone mass, while simultaneously reducing the incidence of stress-shielding compared to the Alloclassic, revealing a key benefit of the distinctive geometry. Several important factors for stable fixation, determined in clinical evaluations from the literature, were evident in both designs: high levels of proximal bone loss and distal bone densification. The results suggest this novel computational framework can be utilized for comparative hip prosthesis shape, uniquely considering the post-operation bone remodeling as a design criterion.     

Computational simulation of simultaneous cortical and trabecular bone change in human proximal femur during bone remodeling

In this study, we developed a numerical framework that computationally determines simultaneous and interactive structural changes of cortical and trabecular bone types during bone remodeling, and we investigated the structural correlation between the two bone types in human proximal femur. We implemented a surface remodeling technique that performs bone remodeling in the exterior layer of the cortical bone while keeping its interior area unchanged. A micro-finite element (μFE) model was constructed that represents the entire cortical bone and full trabecular architecture in human proximal femur. This study simulated and compared the bone adaptation processes of two different structures: (1) femoral bone that has normal cortical bone shape and (2) perturbed femoral bone that has an artificial bone lump in the inferomedial cortex. Using the proposed numerical method in conjunction with design space optimization, we successfully obtained numerical results that resemble actual human proximal femur. The results revealed that actual cortical bone, as well as the trabecular bone, in human proximal femur has structurally optimal shapes, and it was also shown that a bone abnormality that has little contribution to bone structural integrity tends to disappear. This study also quantitatively determined the structural contribution of each bone: when the trabecular adaptation was complete, the trabecular bone supported 54% of the total load in the human proximal femur while the cortical bone carried 46%.     

Functional adaptation of cancellous bone in human proximal femur predicted by trabecular surface remodeling simulation toward uniform stress state

Two-dimensional simulation of trabecular surface remodeling was conducted for a human proximal femur to investigate the structural change of cancellous bone toward a uniform stress state. Considering that a local mechanical stimulus plays an important role in cellular activities in bone remodeling, local stress nonuniformity was assumed to drive trabecular structural change to seek a uniform stress state. A large-scale pixel-based finite element model was used to simulate structural changes of individual trabeculae over the entire bone. As a result, the initial structure of trabeculae changed from isotropic to anisotropic due to trabecular microstructural changes caused by surface remodeling according to the mechanical environment in the proximal femur. Under a single-loading condition, it was shown that the apparent structural property evaluated by fabric ellipses corresponded to the apparent stress state in cancellous bone. As is observed in the actual bone, a distributed trabecular structure was obtained under a multiple-loading condition. Through these studies, it was concluded that trabecular surface remodeling toward a local uniform stress state at the trabecular level could naturally bring about functional adaptation phenomenon at the apparent tissue level. The proposed simulation model would be capable of providing insight into the hierarchical mechanism of trabecular surface remodeling at the microstructural level up to the apparent tissue level.     

Trabecular bone ontogeny in the human proximal femur

Ontogenetic changes in the human femur associated with the acquisition of bipedal locomotion, especially the development of the bicondylar angle, have been well documented. The purpose of this study is to quantify changes in the three-dimensional structure of trabecular bone in the human proximal femur in relation to changing functional and external loading patterns with age. High-resolution X-ray computed tomography scan data were collected for 15 juvenile femoral specimens ranging in age from prenatal to approximately nine years of age. Serial slices were collected for the entire proximal femur of each individual with voxel resolutions ranging from 0.017 to 0.046 mm depending on the size of the specimen. Spherical volumes of interest were defined within the proximal femur, and the bone volume fraction, trabecular thickness, trabecular number, and fabric anisotropy were calculated in three dimensions. Bone volume fraction, trabecular number, and degree of anisotropy decrease between the age of 6 months and 12 months, with the lowest values for these parameters occurring in individuals near 12 months of age. By age 2-3 years, the bone volume, thickness, and degree of anisotropy increase slightly, and regions in the femoral neck become more anisotropic corresponding to the thickening of the inferior cortical bone of the neck. These results suggest that trabecular structure in the proximal femur reflects the shift in external loading patterns associated with the initiation of unassisted walking in infants.     

Computational simulation of trabecular adaptation progress in human proximal femur during growth

There are a large number of clinical and experimental studies that analyzed trabecular architecture as a result of bone adaptation. However, only a limited amount of quantitative data is currently available on the progress of trabecular adaptation during growth. In this paper, we proposed a two-step numerical simulation method that predicts trabecular adaptation progress during growth using a recently developed topology optimization algorithm, design space optimization (DSO), under the hypothesis that the mechanisms of DSO are functionally equivalent to those of bone adaptation. We applied the proposed scheme to trabecular adaptation simulation in human proximal femur. For the simulation, the full trabecular architecture in human proximal femur was represented by a two-dimensional μFE model with 50 μm resolution. In Step 1, we determined a reference value that regulates trabecular adaptation in human proximal femur. In Step 2, we simulated trabecular adaptation in human proximal femur during growth with the reference value derived in Step 1. We analyzed the architectural and mechanical properties of trabecular patterns through iterations. From the comparison with experimental data in the literature, we showed that in the early growth stage trabecular adaptation was achieved mainly by increasing bone volume fraction (or trabecular thickness), while in the later stage of the development the trabecular architecture gained higher structural efficiency by increasing structural anisotropy with a relatively low level of bone volume fraction (or trabecular thickness). We demonstrated that the proposed numerical framework predicted the growing progress of trabecular bone that has a close correlation with experimental data.     

Three-dimensional topology optimization model to simulate the external shapes of bone

Elucidation of the mechanism by which the shape of bones is formed is essential for understanding vertebrate development. Bones support the body of vertebrates by withstanding external loads, such as those imposed by gravity and muscle tension. Many studies have reported that bone formation varies in response to external loads. An increased external load induces bone synthesis, whereas a decreased external load induces bone resorption. This relationship led to the hypothesis that bone shape adapts to external load. In fact, by simulating this relationship through topology optimization, the internal trabecular structure of bones can be successfully reproduced, thereby facilitating the study of bone diseases. In contrast, there have been few attempts to simulate the external structure of bones, which determines vertebrate morphology. However, the external shape of bones may be reproduced through topology optimization because cells of the same type form both the internal and external structures of bones. Here, we constructed a three-dimensional topology optimization model to attempt the reproduction of the external shape of teleost vertebrae. In teleosts, the internal structure of the vertebral bodies is invariable, exhibiting an hourglass shape, whereas the lateral structure supporting the internal structure differs among species. Based on the anatomical observations, we applied different external loads to the hourglass-shaped part. The simulations produced a variety of three-dimensional structures, some of which exhibited several structural features similar to those of actual teleost vertebrae. In addition, by adjusting the geometric parameters, such as the width of the hourglass shape, we reproduced the variation in the teleost vertebrae shapes. These results suggest that a simulation using topology optimization can successfully reproduce the external shapes of teleost vertebrae. By applying our topology optimization model to various bones of vertebrates, we can understand how the external shape of bones adapts to external loads.     

Compressive properties of trabecular bone in the distal femur

Early loosening and implant migration are two problems that lead to failures in cementless (press-fit) femoral knee components of total knee replacements. To begin to address these early failures, this study determined the anterior-posterior mechanical properties from four locations in the human distal femur. Thirty-three cylindrical specimens were removed perpendicular to the press-fit surface after the surgical cuts on 10 human cadaveric femurs (age 71.5+/-14.2 years) had been made. Compression testing was performed that utilized methods to reduce the effects of end-artifacts. The bone mineral apparent density (BMAD), apparent modulus of elasticity, yield and ultimate stress, and yield and ultimate strain were measured for 28 cylindrical specimens. The apparent modulus, yield and ultimate stress, and yield and ultimate strain each significantly differed (p<0.05) in the superior and inferior locations. Linear and power law relationships between superior and inferior mechanical properties and BMAD were determined. The inferior apparent modulus and stresses were higher than those in the superior locations. These results show that the press-fit fixation characteristics of the femoral knee component differ on the anterior shield and posterior condyles. This information will be useful in the assignment of mechanical properties in finite element models for further investigations of femoral knee components. The property-density relations also have applications for implant design and preoperative assessment of bone strength using clinically available tools.     

A method of quantification of stress shielding in the proximal femur using hierarchical computational modeling

Stress shielding is a biomechanical phenomenon causing adaptive changes in bone strength and stiffness around metallic implants, which potentially lead to implant loosening. Accordingly, there is a need for standard, objective engineering measures of the "stress shielding" performances of an implant that can be employed in the process of computer-aided implant design. To provide and test such measures, we developed hierarchical computational models of adaptation of the trabecular microarchitecture at different sites in the proximal femur, in response to insertion of orthopaedic screws and in response to hypothetical reductions in hip joint and gluteal muscle forces. By identifying similar bone adaptation outcomes from the two scenarios, we were able to quantify the stress shielding caused by screws in terms of analogous hypothetical reductions in hip joint and gluteal muscle forces. Specifically, we developed planar lattice models of trabecular microstructures at five regions of interest (ROI) in the proximal femur. The homeostatic and abnormal loading conditions for the lattices were determined from a finite element model of the femur at the continuum scale and fed to an iterative algorithm simulating the adaptation of each lattice to these loads. When screws were inserted to the femur model, maximal simulated bone loss (17% decrease in apparent density, 10% decrease in thickness of trabeculae) was at the greater trochanter and this effect was equivalent to the effect of 50% reduction in gluteal force and normal hip joint force. We conclude that stress shielding performances can be quantified for different screw designs using model-predicted hypothetical musculoskeletal load fractions that would cause a similar pattern and extent of bone loss to that caused by the implants.     

A method of quantification of stress shielding in the proximal femur using hierarchical computational modeling

Stress shielding is a biomechanical phenomenon causing adaptive changes in bone strength and stiffness around metallic implants, which potentially lead to implant loosening. Accordingly, there is a need for standard, objective engineering measures of the “stress shielding” performances of an implant that can be employed in the process of computer-aided implant design. To provide and test such measures, we developed hierarchical computational models of adaptation of the trabecular microarchitecture at different sites in the proximal femur, in response to insertion of orthopaedic screws and in response to hypothetical reductions in hip joint and gluteal muscle forces. By identifying similar bone adaptation outcomes from the two scenarios, we were able to quantify the stress shielding caused by screws in terms of analogous hypothetical reductions in hip joint and gluteal muscle forces. Specifically, we developed planar lattice models of trabecular microstructures at five regions of interest (ROI) in the proximal femur. The homeostatic and abnormal loading conditions for the lattices were determined from a finite element model of the femur at the continuum scale and fed to an iterative algorithm simulating the adaptation of each lattice to these loads. When screws were inserted to the femur model, maximal simulated bone loss (17% decrease in apparent density, 10% decrease in thickness of trabeculae) was at the greater trochanter and this effect was equivalent to the effect of 50% reduction in gluteal force and normal hip joint force. We conclude that stress shielding performances can be quantified for different screw designs using model-predicted hypothetical musculoskeletal load fractions that would cause a similar pattern and extent of bone loss to that caused by the implants.     

Mechanical properties of metaphyseal bone in the proximal femur

We used a three-point bending test to investigate the structural behavior of 123 rectangular flat plate specimens harvested from the metaphyseal shell of the cervical and intertrochanteric regions of five fresh/frozen human proximal femora. For comparison purposes, 36 specimens of similar geometry were also fabricated from bone of the femoral diaphysis. All specimens were oriented in either the local longitudinal or transverse directions. The mean longitudinal elastic modulus was 9650 +/- 2410 (SD) MPa and demonstrated a 24% decrease from that measured for the diaphysis (12500 +/- 2140 MPa) using the same testing technique. However, the transverse elastic moduli did not differ significantly between the proximal (5470 +/- 1720 MPa) and diaphyseal (5990 +/- 1520 MPa) specimens. The globally averaged values for the ultimate tensile strengths of the metaphyseal shell were 101 +/- 26 MPa in the longitudinal and 50 +/- 12 MPa in the transverse directions. These compared with diaphyseal values of 128 +/- 16 MPa and 47 +/- 12 MPa, respectively. While these differences were largely due to the reduced density of the proximal specimens, a slight decrease in transverse anisotropy for the proximal specimens was also noted by comparing the ratio of longitudinal to transverse moduli (1.76) and tensile strength (2.02) to the diaphyseal values (2.09 and 2.71, respectively). Use of these data should lead to improved performance of analytical models for the proximal femur, and thus help focus increased attention on the structural contribution of trabecular bone to the strength and rigidity of the proximal femur.     

Mathematical analysis of trabecular 'trajectories' in apparent trajectorial structures: the unfortunate historical emphasis on the human proximal femur

Wolff's "law" of the functional adaptation of bone is rooted in the trajectory hypothesis of cancellous bone architecture. Wolff often used the human proximal femur as an example of a trajectorial structure (i.e. arched trabecular patterns appear to be aligned along tension/compression stress trajectories). We examined two tenets of the trajectory hypothesis; namely, that the trabecular tracts from the tension- and compression-loaded sides of a bending environment will: (1) follow 'lines' (trajectories) of tension/compression stress that resemble an arch with its apex on a neutral axis, and (2) form orthogonal (90 degrees ) intersections. These predictions were analysed in proximal femora of chimpanzees and modern humans, and in calcanei of sheep and deer. Compared to complex loading of the human femoral neck, the chimpanzee femoral neck reputedly receives relatively simpler loading (i.e. temporally/spatially more consistent bending), and the artiodactyl calcaneus is even more simply loaded in bending. In order to directly consider Wolff's observations, measurements were also made on two-dimensional, cantilevered beams and curved beams, each with intersecting compression/tension stress trajectories. Results in the calcanei showed: (1) the same nonlinear equation best described the dorsal ("compression") and plantar ("tension") trabecular tracts, (2) these tracts could be exactly superimposed on the corresponding compression/tension stress trajectories of the cantilevered beams, and (3) trabecular tracts typically formed orthogonal intersections. In contrast, trabecular tracts in human and chimpanzee femoral necks were non-orthogonal (mean approximately 70 degrees ), with shapes differing from trabecular tracts in calcanei and stress trajectories in the beams. Although often being described by the same equations, the trajectories in the curved beams had lower r(2) values than calcaneal tracts. These results suggest that the trabecular patterns in the calcanei and stress trajectories in short beams are consistent with basic tenets of the trajectory hypothesis while those in human and chimpanzee femoral necks are not. Compared to calcanei, the more complexly loaded human and chimpanzee femoral necks probably receive more prevalent/predominant shear, which is best accommodated by non-orthogonal, asymmetric trabecular tracts. The asymmetrical trabecular patterns in the proximal femora may also reflect the different developmental 'fields' (trochanteric vs. neck/head) that formed these regions, of which there is no parallel in the calcanei.     

Gender-related changes in three-dimensional microstructure of trabecular bone at the human proximal tibia with aging

Despite increasing interest in age- and gender-related bone alterations, data on trabecular microstructure at the proximal tibia are scarce. The aim of this study was to identify trabecular microstructural change at the human proximal tibia with age and gender, using micro-computed tomography (micro-CT) and scanning electron microscopy (SEM). Fifty-six proximal tibias from 28 Japanese men and women (57-98 years of age) were used in this study. The subjects were chosen to give an even age and gender distribution. Both women and men were divided into three age groups, middle (57-68 years), old (72-82 years) and elderly (87-98 years) groups. The trabecular bone specimens from the medial compartment of the proximal tibial metaphysis were examined. Trabecular bone mineral density (BMD), bone volume fraction (BV/TV) and trabecular thickness (Tb.Th) decreased between the middle-aged and elderly groups similarly in women and men. However, trabecular number (Tb.N) decreased by 13% between the middle-aged and elderly groups in women and nearly double that in men. As compared with women, men had higher BV/TV and lower trabecular separation (Tb.Sp) in the old age and elderly groups, and higher Tb.N and connectivity density (Conn.D) in the elderly group. Increased trabecular resorbing surfaces, perforated or disconnected trabeculae and microcallus formations were observed with age. These findings indicate that both BMD and BV/TV decreased at the proximal tibia with age similarly for women and men, but significant differences between women and men were observed for some microstructural parameters. These findings illustrate potential mechanisms underlying osteoporotic proximal tibial fracture.     

Elastic modulus and hardness of cortical and trabecular bone lamellae measured by nanoindentation in the human femur.

The mechanical properties of bone tissue are determined by composition as well as structural, microstructural and nanostructural organization. The aim of this study was to quantify the elastic properties of bone at the lamellar level and compare these properties among osteonal, interstitial and trabecular microstructures from the diaphysis and the neck of the human femur. A nanoindentation technique with a custom irrigation system was used for simultaneously measuring force and displacement of a diamond tip pressed 500 nm into the moist bone tissue. An isotropic elastic modulus was calculated from the unloading curve with an assumed Poisson ratio of 0.3, while hardness was defined as the maximal force divided by the corresponding contact area. The elastic moduli ranged from 6.9 +/- 4.3 GPa in trabecular tissue from the femoral neck of a 74 yr old female up to 25.0 +/- 4.3 GPa in interstitial tissue from the diaphyseal cortex of a 69 yr old female. The mean elastic modulus was found to be significantly influenced by the type of lamella (p < 10(-6)) and by donor (p < 10(-6)). The interaction between the type of lamella and the donor was also highly significant (p < 10(-6)). Hardness followed a similar distribution as elastic modulus among types of lamellae and donor, but with lower statistical contrast. It is concluded that the nanostructure of bone tissue must differ substantially among lamellar types, anatomical sites and individuals and suggests that tissue heterogeneity is of potential importance in bone fragility and adaptation.     

Simulation of the human airways using virtual topology tools and meshing optimization

A method is proposed to improve the quality of the three-dimensional airway geometric models using a commercial software, checking the number of elements, meshing time, and aspect ratio and skewness parameters. The use of real and virtual topologies combined with patch-conforming and patch-independent meshing algorithms results in four different models being the best solution the combination of virtual topology and patch-independent algorithm, due to an excellent aspect ratio and skewness of the elements, and minimum meshing time. The result is a reduction in the computational time required for both meshing and simulation due to a smaller number of cells. The use of virtual topologies combined with patch-independent meshing algorithms could be extended in bioengineering because the geometries handling is similar to this case. The method is applied to a healthy person using their computed tomography images. The resulting numerical models are able to simulate correctly a forced spirometry.     

On the interrelationship of permeability and structural parameters of vertebral trabecular bone: a parametric computational study

To characterise the flow of a fluid through a uniform porous medium, the medium may be completely described by its permeability, a measure of flow resistance. Fluid flow in the intertrabecular spaces of cancellous bone has been recognised as an important factor in a number of physical phenomena. In order to investigate the interdependence of permeability, porosity and the structural parameters, we adapted a morphological model and systematically varied its structural parameters. By simulating a viscous Stokes flow regime, we were able to estimate the anisotropic permeability tensor and performed an extensive, stepwise multivariate regression analysis to establish empirical relationships between the morphological parameters and the permeability for the anatomical directions individually. The regression analysis indicated high values of determination coefficients [0.88 < R² < 0.89 (transversal directions) and R² = 0.60 (longitudinal direction), porosity-based prediction and R² = 0.98 for all directions and information presented to the regression model]. We conclude that a pooled set of structural parameters may explain up to 98% of the permeability variability, the regression model predicts permeability values that match experimental data, and a good prediction performance could be achieved by only incorporating the porosity and either the degree of anisotropy (0.89 < R² < 0.91) or the trabecular spacing predictor (0.96 < R² < 0.97). These conclusions imply that trabecular thickness and shape parameters only play a minor role in the determination of vertebral trabecular bone permeability. However, a major limitation of the model is that it reflects an idealisation of the real, regionally varying structure of trabecular bone. Therefore, the goodness-of-fit estimates we are presenting should be considered as an upper bound limitation regarding the prediction performance.