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.     

Predicting bone remodeling around tissue- and bone-level dental implants used in reduced bone width

The objective of this study was to predict time-dependent bone remodeling around tissue- and bone-level dental implants used in patients with reduced bone width. The remodeling of bone around titanium tissue-level, and titanium and titanium–zirconium alloy bone-level implants was studied under 100 N oblique load for one month by implementing the Stanford theory into three-dimensional finite element models. Maximum principal stress, minimum principal stress, and strain energy density in peri-implant bone and displacement in x- and y- axes of the implant were evaluated. Maximum and minimum principal stresses around tissue-level implant were higher than bone-level implants and both bone-level implants experienced comparable stresses. Total strain energy density in bone around titanium implants slightly decreased during the first two weeks of loading followed by a recovery, and the titanium–zirconium implant showed minor changes in the axial plane. Total strain energy density changes in the loading and contralateral sides were higher in tissue-level implant than other implants in the cortical bone at the horizontal plane. The displacement values of the implants were almost constant over time. Tissue-level implants were associated with higher stresses than bone-level implants. The time-dependent biomechanical outcome of titanium–zirconium alloy bone-level implant was comparable to the titanium implant.     

Numerical analysis of an osseointegrated prosthesis fixation with reduced bone failure risk and periprosthetic bone loss

Currently available implants for direct attachment of prosthesis to the skeletal system after transfemoral amputation (OPRA system, Integrum AB, Sweden and ISP Endo/Exo prosthesis, ESKA Implants AG, Germany) show many advantages over the conventional socket fixation. However, restraining biomechanical issues such as considerable bone loss around the stem and peri-prosthetic bone fractures are present. To overcome these limiting issues a new concept of the direct intramedullary fixation was developed. We hypothesize that the new design will reduce the peri-prosthetic bone failure risk and adverse bone remodeling by restoring the natural load transfer in the femur. Generic CT-based finite element models of an intact femur and amputated bones implanted with 3 analyzed implants were created and loaded with a normal walking and a forward fall load. The strain adaptive bone remodeling theory was used to predict long-term bone changes around the implants and the periprosthetic bone failure risk was evaluated by the von Mises stress criterion. The results show that the new design provides close to physiological distribution of stresses in the bone and lower bone failure risk for the normal walking as compared to the OPRA and the ISP implants. The bone remodeling simulations did not reveal any overall bone loss around the new design, as opposed to the OPRA and the ISP implants, which induce considerable bone loss in the distal end of the femur. This positive outcome shows that the presented concept has a potential to considerably improve safety of the rehabilitation with the direct fixation implants.     

Trabecular surface remodeling simulation for cancellous bone using microstructural voxel finite element models

A computational simulation method for three-dimensional trabecular surface remodeling was proposed, using voxel finite element models of cancellous bone, and was applied to the experimental data. In the simulation, the trabecular microstructure was modeled based on digital images, and its morphological changes due to surface movement at the trabecular level were directly expressed by removing/adding the voxel elements from/to the trabecular surface. A remodeling simulation at the single trabecular level under uniaxial compressive loading demonstrated smooth morphological changes even though the trabeculae were modeled with discrete voxel elements. Moreover, the trabecular axis rotated toward the loading direction with increasing stiffness, simulating functional adaptation to the applied load. In the remodeling simulation at the trabecular structural level, a cancellous bone cube was modeled using a digital image obtained by microcomputed tomography (microCT), and was uniaxially compressed. As a result, the apparent stiffness against the applied load increased by remodeling, in which the trabeculae reoriented to the loading direction. In addition, changes in the structural indices of the trabecular architecture coincided qualitatively with previously published experimental observations. Through these studies, it was demonstrated that the newly proposed voxel simulation technique enables us to simulate the trabecular surface remodeling and to compare the results obtained using this technique with the in vivo experimental data in the investigation of the adaptive bone remodeling phenomenon.     

Finite element analysis of peri-implant bone volume affected by stresses around Morse taper implants: effects of implant positioning to the bone crest

Abstract

Objectives: The purpose of the present study was to evaluate the distribution and magnitude of stresses through the bone tissue surrounding Morse taper dental implants at different positioning relative to the bone crest. Materials and Methods: A mandibular bone model was obtained from a computed tomography scan. A three-dimensional (3D) model of Morse taper implant-abutment systems placed at the bone crest (equicrestal) and 2?mm bellow the bone crest (subcrestal) were assessed by finite element analysis (FEA). FEA was carried out on axial and oblique (45°) loading at 150 N relatively to the central axis of the implant. The von Mises stresses were analysed considering magnitude and volume of affected peri-implant bone. Results: On vertical loading, maximum von Mises stresses were recorded at 6-7?MPa for trabecular bone while values ranging from 73 up to 118?MPa were recorded for cortical bone. On oblique loading at the equiquestral or subcrestal positioning, the maximum von Mises stresses ranged from 15 to 21?MPa for trabecular bone while values at 150?MPa were recorded for the cortical bone. On vertical loading, >99.9vol.% cortical bone volume was subjected to a maximum of 2?MPa while von Mises stress values at 15?MPa were recorded for trabecular bone. On oblique loading, >99.9vol.% trabecular bone volume was subjected to maximum stress values at 5?MPa, while von Mises stress values at 35?MPa were recorded for >99.4vol.% cortical bone. Conclusions: Bone volume-based stress analysis revealed that most of the bone volume (>99% by vol) was subjected to significantly lower stress values around Morse taper implants placed at equicrestal or subcrestal positioning. Such analysis is commentary to the ordinary biomechanical assessment of dental implants concerning the stress distribution through peri-implant sites.     

Predicting changes in mechanical properties of trabecular bone by adaptive remodeling

Because changes in the mechanical properties of bone are closely related to trabecular bone remodeling, methods that consider the temporal morphological changes induced by adaptive remodeling of trabecular bone are needed to estimate long-term fracture risk and bone quality in osteoporosis. We simulated bone remodeling using simplified and pig trabecular bone models and estimated the morphology of healthy and osteoporotic cases. We then displayed the fracture risk of the remodeled models based on a cumulative histogram from high stress. The histogram showed more elements had higher stresses in the osteoporosis model, indicating that the osteoporosis model had a greater risk.     

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.     

Incorporating tissue anisotropy and heterogeneity in finite element models of trabecular bone altered predicted local stress distributions

Trabecular bone is composed of organized mineralized collagen fibrils, which results in heterogeneous and anisotropic mechanical properties at the tissue level. Recently, biomechanical models computing stresses and strains in trabecular bone have indicated a significant effect of tissue heterogeneity on predicted stresses and strains. However, the effect of the tissue-level mechanical anisotropy on the trabecular bone biomechanical response is unknown. Here, a computational method was established to automatically impose physiologically relevant orientation inherent in trabecular bone tissue on a trabecular bone microscale finite element model. Spatially varying tissue-level anisotropic elastic properties were then applied according to the bone mineral density and the local tissue orientation. The model was used to test the hypothesis that anisotropy in both homogeneous and heterogeneous models alters the predicted distribution of stress invariants. Linear elastic finite element computations were performed on a 3 mm cube model isolated from a microcomputed tomography scan of human trabecular bone from the distal femur. Hydrostatic stress and von Mises equivalent stress were recorded at every element, and the distributions of these values were analyzed. Anisotropy reduced the range of hydrostatic stress in both tension and compression more strongly than the associated increase in von Mises equivalent stress. The effect of anisotropy was independent of the spatial redistribution high compressive stresses due to tissue elastic heterogeneity. Tissue anisotropy and heterogeneity are likely important mechanisms to protect bone from failure and should be included for stress analyses in trabecular bone.     

Changes in the Fabric and Compliance Tensors of Cancellous Bone due to Trabecular Surface Remodeling,Predicted by a Digital Image-based Model

Fabric and compliance tensors of a cube of cancellous bone with a complicated three-dimensional trabecular structure were obtained for trabecular surface remodeling by using a digital image-based model combined with a large-scale finite element method. Using mean intercept length and a homogenization method, the fabric and compliance tensors were determined for the trabecular structure obtained in the computer remodeling simulation. The tensorial quantities obtained indicated that anisotropic structural changes occur in cancellous bone adapting to the compressive loading condition. There were good correlations between the fabric tensor, bone volume fraction, and compliance tensor in the remodeling process. The result demonstrates that changes in the structural and mechanical properties of cancellous bone are essentially anisotropic and should be expressed by tensorial quantities.     

Changes in the fabric and compliance tensors of cancellous bone due to trabecular surface remodeling, predicted by a digital image-based model

Fabric and compliance tensors of a cube of cancellous bone with a complicated three-dimensional trabecular structure were obtained for trabecular surface remodeling by using a digital image-based model combined with a large-scale finite element method. Using mean intercept length and a homogenization method, the fabric and compliance tensors were determined for the trabecular structure obtained in the computer remodeling simulation. The tensorial quantities obtained indicated that anisotropic structural changes occur in cancellous bone adapting to the compressive loading condition. There were good correlations between the fabric tensor, bone volume fraction, and compliance tensor in the remodeling process. The result demonstrates that changes in the structural and mechanical properties of cancellous bone are essentially anisotropic and should be expressed by tensorial quantities.     

The discrete nature of trabecular bone microarchitecture affects implant stability

Small endosseous implants, such as screws, are important components of modern orthopedics and dentistry. Hence they have to reliably fulfill a variety of requirements, which makes the development of such implants challenging. Finite element analysis is a widely used computational tool used to analyze and optimize implant stability in bone. For these purposes, bone is generally modeled as a continuum material. However, bone failure and bone adaptation processes are occurring at the discrete level of individual trabeculae; hence the assessment of stresses and strains at this level is relevant. Therefore, the aim of the present study was to investigate how peri-implant strain distribution and load transfer between implant and bone are affected by the continuum assumption. We performed a computational study in which cancellous screws were inserted in continuum and discrete models of trabecular bone; axial loading was simulated. We found strong differences in bone-implant stiffness between the discrete and continuum bone model. They depended on bone density and applied boundary conditions. Furthermore, load transfer from the screw to the surrounding bone differed strongly between the continuum and discrete models, especially for low-density bone. Based on our findings we conclude that continuum bone models are of limited use for finite element analysis of peri-implant mechanical loading in trabecular bone when a precise quantification of peri-implant stresses and strains is required. Therefore, for the assessment and improvement of trabecular bone implants, finite element models which accurately represent trabecular microarchitecture should be used.     

Implementation of contact definitions calculated by FEA to describe the healing process of basal implants

Aims: Bone structure around basal implants shows a dual healing mode: direct contact areas manifest primary osteonal remodeling, in the void osteotomy-induced spaces, the repair begins with woven bone formation. This woven bone is later converted into osteonal bone. The purpose of this study was to develop a model to accurately represent the interface between bone and basal implant throughout the healing process. The model was applied to the biological scenario of changing load distribution in a basal implant system over time. Methods: Computations were made through finite element analysis using multiple models with changing boneimplant contact definitions which reflected the dynamic nature of the interface throughout the bony healing process. Five stages of bony healing were calculated taking into account the changes in mineral content of bone in the vicinity of the load transmitting implant surfaces. Results: As the bony integration of basal implants proceeds during healing, peak stresses within the metal structure shift geographically. While bony repair may still weaken osteonal bone, woven bone has already matured. This leads to changes in the load distribution between and within the direct contact areas, and bone areas which make later contact with implant. Conclusions: This study shows that basal implants undergo an intrinsic shift of maximum stress regions during osseointegration. Fatigue testing methods in the case of basal implants must therefore take into account this gradual shift from early healing phase until full osseointegration is achieved.     

Numerical simulation of the remodelling process of trabecular architecture around dental implants

Dental implants may alter the mechanical environment in the jawbone, thereby causing remodelling and adaptation of the surrounding trabecular bone tissues. To improve the efficacy of dental implant systems, it is necessary to consider the effect of bone remodelling on the performance of the prosthetic systems. In this study, finite element simulations were implemented to predict the evolution of microarchitecture around four implant systems using a previously developed model that combines both adaptive and microdamage-based mechano-sensory mechanisms in bone remodelling process. Changes in the trabecular architecture around dental implants were mainly focused. The simulation results indicate that the orientational and ladder-like architecture around the implants predicted herein is in good agreement with those observed in animal experiments and clinical observations. The proposed algorithms were shown to be effective in simulating the remodelling process of trabecular architecture around dental implant systems. In addition, the architectural features around four typical dental implant systems in alveolar bone were evaluated comparatively.     

A Femur-Implant Model for the Prediction of Bone Remodeling Behavior Induced by Cementless Stem

Bone remodeling simulation is an effective tool for the prediction of long-term effect of implant on the bone tissue, as well as the selection of an appropriate implant in terms of architecture and material. In this paper, a finite element model of proximal femur was developed to simulate the structures of internal trabecular and cortical bones by incorporating quantitative bone functional adaptation theory with finite element analysis. Cementless stems made of titanium, two types of Functionally Graded Material (FGM) and flexible ‘iso-elastic’ material as comparison were implanted in the structure of proximal femur respectively to simulate the bone remodeling behaviors of host bone. The distributions of bone density, von Mises stress, and interface shear stress were obtained. All the prosthetic stems had effects on the bone remodeling behaviors of proximal femur, but the degrees of stress shielding were different. The amount of bone loss caused by titanium implant was in agreement with the clinical observation. The FGM stems caused less bone loss than that of the titanium stem, in which FGM I stem (titanium richer at the top to more HAP/Col towards the bottom) could relieve stress shielding effectively, and the interface shear stresses were more evenly distributed in the model with FGM I stem in comparison with those in the models with FGM II (titanium and bioglass) and titanium stems. The numerical simulations in the present study provided theoretical basis for FGM as an appropriate material of femoral implant from a biomechanical point of view. The next steps are to fabricate FGM stem and to conduct animal experiments to investigate the effects of FGM stem on the remodeling behaviors using animal model.     

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.     

Influence of the stiffness of bone defect implants on the mechanical conditions at the interface--a finite element analysis with contact

The study focused on the influence of the implant material stiffness on stress distribution and micromotion at the interface of bone defect implants. We hypothesized that a low-stiffness implant with a modulus closer to that of the surrounding trabecular bone would yield a more homogeneous stress distribution and less micromotion at the interface with the bony bed. To prove this hypothesis we generated a three-dimensional, non-linear, anisotropic finite element (FE) model. The FE model corresponded to a previously developed animal model in sheep. A prismatic implant filled a standardized defect in the load-bearing area of the trabecular bone beneath the tibial plateau. The interface was described by face-to-face contact elements, which allow press fits, friction, sliding, and gapping. We assumed a physiological load condition and calculated contact pressures, shear stresses, and shear movements at the interface for two implants of different stiffness (titanium: E=110GPa; composite: E=2.2GPa). The FE model showed that the stress distribution was more homogeneous for the low-stiffness implant. The maximum pressure for the composite implant (2.1 MPa) was lower than for the titanium implant (5.6 MPa). Contrary to our hypothesis, we found more micromotion for the composite (up to 6 microm) than for the titanium implant (up to 4.5 microm). However, for both implants peak stresses and micromotion were in a range that predicts adequate conditions for the osseointegration. This was confirmed by the histological results from the animal studies.     

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.     

Mineral heterogeneity affects predictions of intratrabecular stress and strain

Knowledge of the influence of mineral variations (i.e., mineral heterogeneity) on biomechanical bone behavior at the trabecular level is limited. The aim of this study is to investigate how this material property affects the intratrabecular distributions of stress and strain in human adult trabecular bone. Two different sets of finite element (FE) models of trabecular samples were constructed; tissue stiffness was either scaled to the local degree of mineralization of bone as measured with microCT (heterogeneous) or tissue stiffness was assumed to be homogeneous. The influence of intratrabecular mineral heterogeneity was analyzed by comparing both models. Interesting effects were seen regarding intratrabecular stress and strain distributions. In the homogeneous model, the highest stresses were found at the surface with a significant decrease towards the core. Higher superficial stresses could indicate a higher predicted fracture risk in the trabeculae. In the heterogeneous model this pattern was different. A significant increase in stress with increasing distance from the trabecular surface was found followed by a significant decrease towards the core. This suggests trabecular bending during a compression. In both models a decrease in strain values from surface to core was predicted, which is consistent with trabecular bending. When mineral heterogeneity was taken into account, the predicted intratrabecular patterns of stress and strain are more consistent with the expected biomechanical behavior as based on mineral variations in trabeculae. Our findings indicate that mineral heterogeneity should not be neglected when performing biomechanical studies on topics such as the (long-term or dose dependent) effects of antiresorptive treatments.     

Trabecular bone microdamage and microstructural stresses under uniaxial compression

The balance between local remodeling and accumulation of trabecular bone microdamage is believed to play an important role in the maintenance of skeletal integrity. However, the local mechanical parameters associated with microdamage initiation are not well understood. Using histological damage labeling, micro-CT imaging, and image-based finite element analysis, regions of trabecular bone microdamage were detected and registered to estimated microstructural von Mises effective stresses and strains, maximum principal stresses and strains, and strain energy density (SED). Bovine tibial trabecular bone cores underwent a stepwise uniaxial compression routine in which specimens were micro-CT imaged following each compression step. The results indicate that the mode of trabecular failure observed by micro-CT imaging agreed well with the polarity and distribution of stresses within an individual trabecula. Analysis of on-axis subsections within specimens provided significant positive relationships between microdamage and each estimated tissue stress, strain and SED parameter. In a more localized analysis, individual microdamaged and undamaged trabeculae were extracted from specimens loaded within the elastic region and to the apparent yield point. As expected, damaged trabeculae in both groups possessed significantly higher local stresses and strains than undamaged trabeculae. The results also indicated that microdamage initiation occurred prior to apparent yield at local principal stresses in the range of 88-121 MPa for compression and 35-43 MPa for tension and local principal strains of 0.46-0.63% in compression and 0.18-0.24% in tension. These data provide an important step towards understanding factors contributing to microdamage initiation and establishing local failure criteria for normal and diseased trabecular bone.     

Tissue stresses and strain in trabeculae of a canine proximal femur can be quantified from computer reconstructions

A quantitative assessment of bone tissue stresses and strains is essential for the understanding of failure mechanisms associated with osteoporosis, osteoarthritis, loosening of implants and cell- mediated adaptive bone-remodeling processes. According to Wolff's trajectorial hypothesis, the trabecular architecture is such that minimal tissue stresses are paired with minimal weight. This paradigm at least suggests that, normally, stresses and strains should be distributed rather evenly over the trabecular architecture. Although bone stresses at the apparent level were determined with finite element analysis (FEA), by assuming it to be continuous, there is no data available on trabecular tissue stresses or strains of bones in situ under physiological loading conditions. The objectives of this project were to supply reasonable estimates of these quantities for the canine femur, to compare trabecular-tissue to apparent stresses, and to test Wolff's hypothesis in a quantitative sense. For that purpose, the newly developed method of large-scale micro-FEA was applied in conjunction with micro-CT structural measurements. A three-dimensional high-resolution computer reconstruction of a proximal canine femur was made using a micro-CT scanner. This was converted to a large-scale FE-model with 7.6 million elements, adequately refined to represent individual trabeculae. Using a special-purpose FE-solver, analyses were conducted for three different orthogonal hip-joint loading cases, one of which represented the stance-phase of walking. By superimposing the results, the tissue stress and strain distributions could also be calculated for other force directions. Further analyses of results were concentrated on a trabecular volume of interest (VOI) located in the center of the head. For the stance phase of walking an average tissue principal strain in the VOI of 279 strain was found, with a standard deviation of 212 microstrain. The standard deviation depended not only on the hip-force magnitude, but also on its direction. In more than 95% of the tissue volume the principal stresses and strains were in a range from zero to three times the averages, for all hip-force directions. This indicates that no single load creates even stress or strain distributions in the trabecular architecture. Nevertheless, excessive values occurred at few locations only, and the maximum tissue stress was approximately half the value reported for the tissue fatigue strength. These results thus indicate that trabecular bone tissue has a safety factor of approximately two for hip-joint loads that occur during normal activities.