Trabecular bone adaptation with an orthotropic material model.

Most bone adaptation algorithms, that attempt to explain the connection between bone morphology and loads, assume that bone is effectively isotropic. An isotropic material model can explain the bone density distribution, but not the structure and pattern of trabecular bone, which clearly has a mechanical significance. In this paper, an orthotropic material model is utilized to predict the proximal femur trabecular structure. Two hypotheses are combined to determine the local orientation and material properties of each element in the model. First, it is suggested that trabecular directions, which correspond to the orthotropic material axes, are determined locally by the maximal principal stress directions due to the multiple load cases (MLC) the femur is subject to. The second hypothesis is that material properties in each material direction can be determined using directional stimuli, thus extending existing adaptation algorithms to include directionality. An algorithm is utilized, where each iteration comprises of two stages. First, material axes are rotated to the direction of the largest principal stress that occurs from a multiple load scheme applied to the proximal femur. Next, material properties are modified in each material direction, according to a directional stimulus. Results show that local material directions correspond with known trabecular patterns, reproducing all main groups of trabeculae very well. The local directional stiffnesses, degree of anisotropy and density distribution are shown to conform to real femur morphology.     

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

Prediction of shape and internal structure of the proximal femur using a modified level set method for structural topology optimisation

A computational framework was developed to simulate the bone remodelling process as a structural topology optimisation problem. The mathematical formulation of the Level Set technique was extended and then implemented into a coronal plane model of the proximal femur to simulate the remodelling of internal structure and external geometry of bone into the optimal state. Results indicated that the proposed approach could reasonably mimic the major geometrical and material features of the natural bone. Simulation of the internal bone remodelling on the typical gross shape of the proximal femur, resulted in a density distribution pattern with good consistency with that of the natural bone. When both external and internal bone remodelling were simulated simultaneously, the initial rectangular design domain with a regularly distributed mass reduced gradually to an optimal state with external shape and internal structure similar to those of the natural proximal femur.     

Prediction of shape and internal structure of the proximal femur using a modified level set method for structural topology optimisation

A computational framework was developed to simulate the bone remodelling process as a structural topology optimisation problem. The mathematical formulation of the Level Set technique was extended and then implemented into a coronal plane model of the proximal femur to simulate the remodelling of internal structure and external geometry of bone into the optimal state. Results indicated that the proposed approach could reasonably mimic the major geometrical and material features of the natural bone. Simulation of the internal bone remodelling on the typical gross shape of the proximal femur, resulted in a density distribution pattern with good consistency with that of the natural bone. When both external and internal bone remodelling were simulated simultaneously, the initial rectangular design domain with a regularly distributed mass reduced gradually to an optimal state with external shape and internal structure similar to those of the natural proximal femur.     

Tensile strength of bovine trabecular bone

Data on the tensile and compressive properties of trabecular bone are needed to define input parameters and failure criteria for modeling total joint replacements. To help resolve differences in reports comparing tensile and compressive properties of trabecular bone, we have developed new methods, based on porous foam technology, for tensile testing of fresh/frozen trabecular bone specimens. Using bovine trabecular bone from an isotropic region from the proximal humerus as a model material, we measured ultimate strengths in tension and compression for two groups of 24 specimens each. The average ultimate strength in tension was 7.6 +/- 2.2 (95% C.I.) MPa and in compression was 12.4 +/- 3.2 MPa. This difference was statistically significant (p = 0.013) and was not related to density differences between the test groups (p = 0.28). Strength was related by a power-law function of the local apparent density, but, even accounting for density influences, isotropic bovine trabecular bone exhibits significantly lower strengths in tension than in compression.     

The elastic and ultimate properties of compact bone tissue

The use of a tranversely isotropic model is tested for the elastic behavior of bovine and human bone and the five independent constants of this model are determined. The accuracy of the model is tested for eight cases by comparing the off-axis modulus predicted by a rotation of stiffness matrix with an experimentally determined off-axis modulus.

Ultimate properties are presented for bovine and human bone for tension, compression, and torsional loads. A Hankinson type failure criterion is proposed for off-axis ultimate stress and this predicted value compared with experimental values for nine cases.     

Compression or tension? The stress distribution in the proximal femur

Background  

Questions regarding the distribution of stress in the proximal human femur have never been adequately resolved. Traditionally, by considering the femur in isolation, it has been believed that the effect of body weight on the projecting neck and head places the superior aspect of the neck in tension. A minority view has proposed that this region is in compression because of muscular forces pulling the femur into the pelvis. Little has been done to study stress distributions in the proximal femur. We hypothesise that under physiological loading the majority of the proximal femur is in compression and that the internal trabecular structure functions as an arch, transferring compressive stresses to the femoral shaft.     

Nonlinear behavior of trabecular bone at small strains

Study of the behavior of trabecular bone at strains below 0.40 percent is of clinical and biomechanical importance. The goal of this work was to characterize, with respect to anatomic site, loading mode, and apparent density, the subtle concave downward stress-strain nonlinearity, that has been observed recently for trabecular bone at these strains. Using protocols designed to minimize end-artifacts, 155 cylindrical cores from human vertebrae, proximal tibiae, proximal femora, and bovine proximal tibiae were mechanically tested to yield at 0.50 percent strain per second in tension or compression. The nonlinearity was quantified by the reduction in tangent modulus at 0.20 percent and 0.40 percent strain as compared to the initial modulus. For the pooled data, the mean +/- SD percentage reduction in tangent modulus at 0.20 percent strain was 9.07+/- 3.24 percent in compression and 13.8 +/- 4.79 percent in tension. At 0.40 percent strain, these values were 23.5 +/- 5.71 and 35.7+/- 7.10 percent, respectively. The magnitude of the nonlineari't depended on both anatomic site (p < 0.001) and loading mode (p < 0.001), and in tension was positively correlated with density. Calculated values of elastic modulus and yield properties depended on the strain range chosen to define modulus via a linear curve fit (p < 0.005). Mean percent differences in 0.20 percent offset yield strains were as large as 10.65 percent for some human sites. These results establish that trabecular bone exhibits nonlinearity at low strains, and that this behavior can confound intersite comparisons of mechanical properties. A nonlinear characterization of the small strain behavior of trabecular bone was introduced to characterize the initial stress-strain behavior more thoroughly.     

Stress analysis and failure prediction in the proximal femur before and after total hip replacement

An investigation was performed to determine the effects of the presence of two lengths of proximal Müller prosthesis on predicted failure loads, as compared to those for an intact femur. Three-dimensional stresses in a bone/cement/prosthesis system were determined using finite element methods, with both isotropic and transversely isotropic material properties used for the diaphyseal cortex. Significant increases in prosthesis stem stresses were found when the transversely isotropic material properties were employed in the diaphyseal cortex. This leads to the conclusion that accurate anisotropic material properties for bone are essential for precise stress determination and optimum design in prosthetic implants. Failure loads were also predicted for vertical compression and axial torque, similar to available experimental conditions, and were within the range of the experimental failure data found in the literature. The technique developed herein can be used to systematically assess existing as well as future implant designs, taking into account the complex three-dimensional interaction effects of the overall bone/cement/prosthesis system.     

Elastic area compressibility modulus of red cell membrane.

Micropipette measurements of isotropic tension vs. area expansion in pre-swollen single human red cells gave a value of 288 +/- 50 SD dyn/cm for the elastic, area compressibility modulus of the total membrane at 25 degrees C. This elastic constant, characterizing the resistance to area expansion or compression, is about 4 X 10(4) times greater than the elastic modulus for shear rigidity; therefore, in situations where deformation of the membrane does not require large isotropic tensions (e.g., in passage through normal capillaries), the membrane can be treated by a simple constitutive relation for a two-dimensionally, incompressible material (i.e. fixed area). The tension was found to be linear and reversible for the range of area changes observed (within the experimental system resolution of 10%). The maximum fractional area expansion required to produce lysis was uniformly distributed between 2 and 4% with 3% average and 0.7% SD. By heating the cells to 50 degrees C, it appears that the structural matrix (responsible for the shear rigidity and most of the strength in isotropic tension) is disrupted and primarily the lipid bilayer resists lysis. Therefore, the relative contributions of the structural matrix and lipid bilayer to the elastic, area compressibility could be estimated. The maximum isotropic tension at 25 degrees C is 10-12 dyn/cm and at 50 degrees C is between 3 and 4 dyn/cm. From this data, the respective compressibilities are estimated at 193 dyn/cm and 95 dyn/cm for structural network and bilayer. The latter value correlates well with data on in vitro, monolayer surface pressure versus area curves at oil-water interfaces.     

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.     

Reliable simulations of the human proximal femur by high-order finite element analysis validated by experimental observations

Background: The mechanical response of patient-specific bone to various load conditions is of major clinical importance in orthopedics. Herein we enhance the methods presented in Yosibash et al. [2007. A CT-based high-order finite element analysis of the human proximal femur compared to in-vitro experiments. ASME Journal of Biomechanical Engineering 129(3), 297–309.] for the reliable simulations of the human proximal femur by high-order finite elements (FEs) and validate the simulations by experimental observations.

Method of approach: A fresh-frozen human femur was scanned by quantitative computed tomography (QCT) and thereafter loaded (in vitro experiments) by a quasi-static force of up to 1250 N. QCT scans were manipulated to generate a high-order FE bone model with distinct cortical and trabecular regions having inhomogeneous isotropic elastic properties with Young's modulus represented by continuous spatial functions. Sensitivity analyses were performed to quantify parameters that mostly influence the mechanical response. FE results were compared to displacements and strains measured in the experiments.

Results: Young moduli correlated to QCT Hounsfield Units by relations in Keyak and Falkinstein [2003. Comparison of in situ and in vitro CT scan-based finite element model predictions of proximal femoral fracture load. Medical Engineering and Physics 25, 781–787.] were found to provide predictions that match the experimental results closely. Excellent agreement was found for both the displacements and strains. The presented study demonstrates that reliable and validated high-order patient-specific FE simulations of human femurs based on QCT data are achievable for clinical computer-aided decision making.     

A B-spline based heterogeneous modeling and analysis of proximal femur with graded element

Bone is a complex biological tissue and natural heterogeneous object. The main objective of this study is to simulate quasi-static loading of bio-objects like human femur with B-spline based modeling and its 3D finite element analysis with graded element. B-spline surface representation method is extended to represent material composition to develop heterogeneous solid model of proximal femur. Lagrangian graded element is used to assign inhomogeneous isotropic elastic properties in finite element model to improve the performance. Convergence study is carried out with finite element model in single leg stance load condition. To test the feasibility of the model, sensitivity of simulation is investigated. To validate the model, numerical results are compared with those of an experimental work for the same specimen in simple stance load condition obtained from one of the reference paper. Good agreement is achieved for vertical displacement and strains in most of the locations.     

The elastic properties of trabecular and cortical bone tissues are similar: results from two microscopic measurement techniques

Acoustic microscopy (30-60 microm resolution) and nanoindentation (1-5 microm resolution) are techniques that can be used to evaluate the elastic properties of human bone at a microstructural level. The goals of the current study were (1) to measure and compare the Young's moduli of trabecular and cortical bone tissues from a common human donor, and (2) to compare the Young's moduli of bone tissue measured using acoustic microscopy to those measured using nanoindentation. The Young's modulus of cortical bone in the longitudinal direction was about 40% greater than (p<0.01) the Young's modulus in the transverse direction. The Young's modulus of trabecular bone tissue was slightly higher than the transverse Young's modulus of cortical bone, but substantially lower than the longitudinal Young's modulus of cortical bone. These findings were consistent for both measurement methods and suggest that elasticity of trabecular tissue is within the range of that of cortical bone tissue. The calculation of Young's modulus using nanoindentation assumes that the material is elastically isotropic. The current results, i.e., the average anisotropy ratio (E(L)/E(T)) for cortical bone determined by nanoindentation was similar to that determined by the acoustic microscope, suggest that this assumption does not limit nanoindentation as a technique for measurement of Young's modulus in anisotropic bone.     

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.     

The bone tissue of the canine mandible is elastically isotropic

This paper reports experimental measurements which show that canine mandibular bone tissue is elastically isotropic. Earlier work has established that human, canine and bovine cortical bone tissue of the femur, tibia and skull are elastically anisotropic and therefore the reported isotropy of mandibular tissue was unexpected. The isotropic elastic moduli of the canine mandible are represented by a Young's modulus of 7.5 GPa and a Poisson's ratio of 0.4. Earlier work gave the three orthotropic Young's moduli of the cortical one of the canine femur as 12.8 GPa, 15.6 GPa and 20.1 GPa. The experimental technique employed is elastic wave propagation at ultrasonic frequencies.     

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

Automatic determination of anatomical coordinate systems for three-dimensional bone models of the isolated human knee

The combination of three-dimensional (3-D) models with dual fluoroscopy is increasingly popular for evaluating joint function in vivo. Applying these modalities to study knee motion with high accuracy requires reliable anatomical coordinate systems (ACSs) for the femur and tibia. Therefore, a robust method for creating ACSs from 3-D models of the femur and tibia is required. We present and evaluate an automated method for constructing ACSs for the distal femur and proximal tibia based solely on 3-D bone models. The algorithm requires no observer interactions and uses model cross-sectional area, center of mass, principal axes of inertia, and cylindrical surface fitting to construct the ACSs. The algorithm was applied to the femur and tibia of 10 (unpaired) human cadaveric knees. Due to the automated nature of the algorithm, the within specimen variability is zero for a given bone model. The algorithm’s repeatability was evaluated by calculating variability in ACS location and orientation across specimens. Differences in ACS location and orientation between specimens were low (<1.5 mm and <2.5°). Variability arose primarily from natural anatomical and morphological differences between specimens. The presented algorithm provides an alternative method for automatically determining subject-specific ACSs from the distal femur and proximal tibia.     

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.     

Mapping anisotropy of the proximal femur for enhanced image based finite element analysis

Finite element (FE) models of bone derived from quantitative computed tomography (QCT) rely on realistic material properties to accurately predict bone strength. QCT cannot resolve bone microarchitecture, therefore QCT-based FE models lack the anisotropy apparent within the underlying bone tissue. This study proposes a method for mapping femoral anisotropy using high-resolution peripheral quantitative computed tomography (HR-pQCT) scans of human cadaver specimens. Femur HR-pQCT images were sub-divided into numerous overlapping cubic sub-volumes and the local anisotropy was quantified using a ‘direct-mechanics’ method. The resulting directionality reflected all the major stress lines visible within the trabecular lattice, and provided a realistic estimate of the alignment of Harvesian systems within the cortical compartment. QCT-based FE models of the proximal femur were constructed with isotropic and anisotropic material properties, with directionality interpolated from the map of anisotropy. Models were loaded in a sideways fall configuration and the resulting whole bone stiffness was compared to experimental stiffness and ultimate strength. Anisotropic models were consistently less stiff, but no statistically significant differences in correlation were observed between material models against experimental data. The mean difference in whole bone stiffness between model types was approximately 26%, suggesting that anisotropy can still effect considerable change in the mechanics of proximal femur models. The under prediction of whole bone stiffness in anisotropic models suggests that the orthotropic elastic constants require further investigation. The ability to map mechanical anisotropy from high-resolution images and interpolate information into clinical-resolution models will allow testing of new anisotropic material mapping strategies.