Variation of trabecular architecture in proximal femur of postmenopausal women

This investigation of microstructure in the human proximal femur probes the relationship between the parameters of the FRAX index of fracture risk and the parameters of bone microstructure. The specificity of fracture sites at the proximal femur raises the question of whether trabecular parameters are site-specific during post-menopause, before occurrence of fragility fracture. The donated proximal femurs of sixteen post-menopausal women in the sixth and seventh decades of life, free of metabolic pathologies and therapeutic interventions that could have altered the bone tissue, constituted the material of the study. We assessed bone mineral density of the proximal femurs by dual energy X-ray absorptiometry and then sectioned the femurs through the center of the femoral head and along the femoral neck axis. For each proximal femur, morphometry of trabeculae was conducted on the plane of the section divided into conventional regions and sub-regions consistent with the previously identified trabecular families that provide regions of relatively homogeneous microstructure. Mean trabecular width and percent bone area were calculated at such sites. Our findings indicate that each of mean trabecular width and percent bone area vary within each proximal femur independently from each other, with dependence on site. Both trabecular parameters show significant differences between pairs of sites. We speculate that a high FRAX index at the hip corresponds to a reduced percent bone area among sites that gives a more homogeneous and less site-specific quality to the proximal femur. This phenomenon may render the local tissue less able to carry out the expected mechanical function.     

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.     

Multi-organ hydatidosis with extensive involvement of the hemi-pelvis and ipsilateral femur

Bone involvement of hydatid disease is uncommon but when encountered, it presents few unique pathological features. The pattern of tissue involvement is largely different from that of visceral hydatid cyst. We describe the case of a 47 year-old man from northern India, a case of systemic hydatidosis including the liver and the right lung, presenting with an abscess like lesion in the left gluteal region with pathological fractures of the left femur. Radiographs and CT-scan images showed extensive invasion of the left hemi-pelvis and left proximal femur. Debridement of the honeycombed ilium yielded hydatid fluid, numerous small cysts and necrotic material. Multiple large devitalized and sequestrated bone pieces were recovered from the bone cavity of the affected ilium. A histopathological study of the bone sequestrums revealed the unique pattern of bone invasion by the characteristic laminated multi-layered cyst walls into areas of least resistance. Bone sequestration has not often been described or demonstrated elaborately in published studies of the past. The bone defects formed after debridement of the ilium and proximal femur were filled with bone cement along with augmentation of the femur using intra-medullary nail. The surgical technique adopted in our case although was not expected to be curative owing to the multi-system disease; it did result in significant functional improvement in the patient.     

Optimal Mass Distribution Prediction for Human Proximal Femur with Bi-modulus Property

Simulation of the mass distribution in a human proximal femur is important to provide a reasonable therapy scheme for a patient with osteoporosis. An algorithm is developed for prediction of optimal mass distribution in a human proximal femur under a given loading environment. In this algorithm, the bone material is assumed to be bi-modulus, i.e., the tension modulus is not identical to the compression modulus in the same direction. With this bi-modulus bone material, a topology optimization method, i.e., modified SIMP approach, is employed to determine the optimal mass distribution in a proximal femur. The effects of the difference between two moduli on the final material distribution are numerically investigated. Numerical results obtained show that the mass distribution in bi-modular bone materials is different from that in traditional isotropic material. As the tension modulus is less than the compression modulus for bone tissues, the amount of mass required to support tension loads is greater than that required by isotropic material for the same daily activities including one-leg stance, abduction and adduction.     

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.     

Sensitivity analysis of a cemented hip stem to implant position and cement mantle thickness

Patient-specific finite element models of the implanted proximal femur can be built from pre-operative computed tomography scans and post-operative X-rays. However, estimating three-dimensional positioning from two-dimensional radiographs introduces uncertainty in the implant position. Further, accurately measuring the thin cement mantle and the degree of cement–bone interdigitation from imaging data is challenging. To quantify the effect of these uncertainties in stem position and cement thickness, a sensitivity study was performed. A design-of-experiment study was implemented, simulating both gait and stair ascent. Cement mantle stresses and bone–implant interface strains were monitored. The results show that small variations in alignment affect the implant biomechanics, especially around the most proximal and most distal ends of the stem. The results suggest that implant position is more influential than cement thickness. Rotation around the medial–lateral axis is the dominant factor in the proximal zones and stem translations are the dominant factors around the distal tip.     

Design considerations for ceramic resurfaced femoral head: effect of interface characteristics on failure mechanisms

Ceramic hip resurfacing may offer improved wear resistance compared to metallic components. The study is aimed at investigating the effects of stiffer ceramic components on the stress/strain-related failure mechanisms in the resurfaced femur, using three-dimensional finite element models of intact and resurfaced femurs with varying stem–bone interface conditions. Tensile stresses in the cement varied between 1 and 5 MPa. Postoperatively, 20–85% strain shielding was observed inside the resurfaced head. The variability in stem–bone interface condition strongly influenced the stresses and strains generated within the resurfaced femoral head. For full stem–bone contact, high tensile (151–158 MPa) stresses were generated at the cup–stem junction, indicating risk of fracture. Moreover, there was risk of femoral neck fracture due to elevated bone strains (0.60–0.80% strain) in the proximal femoral neck region. Stresses in the ceramic component are reduced if a frictionless gap condition exists at the stem–bone interface. High stresses, coupled with increased strain shielding in the ceramic resurfaced femur, appear to be major concerns regarding its use as an alternative material.     

A finite element analysis of hollow stemmed hip prostheses as a means of reducing stress shielding of the femur.

Stress shielding of the femur is known to be a principal factor in aseptic loosening of hip replacements. This paper considers the use of a hollow stemmed hip implant for reducing the effects of stress shielding, while maintaining acceptably low levels of stress in the cement. Using finite element modelling, the stresses in the proximal femur using different shapes of hollow stem were compared with those produced using comparable sizes of solid stem with different values of elastic modulus. A reduction in stress shielding could be achieved with a hollow stem. A cylindrical hollow stem design was then optimised in order to control the maximum allowable stress in the cement, the minimum allowable stresses in the bone, and a combination of the two. The resulting stems achieved an increase in proximal bone stress of about 15% for the first case and 32% for a model using high strength cement, compared with solid stems of the same nominal outside diameter. The gains of these theoretically optimised designs dropped off rapidly further down the stem. Linearly tapered hollow stems reached a 22% gain, which could be a good compromise between acceptable cement stresses and ease of manufacture.     

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.     

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.     

Critical evaluation of known bone material properties to realize anisotropic FE-simulation of the proximal femur

PURPOSE: In a meta-analysis of the literature we evaluated the present knowledge of the material properties of cortical and cancellous bone to answer the question whether the available data are sufficient to realize anisotropic finite element (FE)-models of the proximal femur. MATERIAL AND METHOD: All studies that met the following criteria were analyzed: Young's modulus, tensile, compressive and torsional strengths, Poisson's ratio, the shear modulus and the viscoelastic properties had to be determined experimentally. The experiments had to be carried out in a moist environment and at room temperature with freshly removed and untreated human cadaverous femurs. All material properties had to be determined in defined load directions (axial, transverse) and should have been correlated to apparent density (g/cm(3)), reflecting the individually variable and age-dependent changes of bone material properties. RESULTS: Differences in Young's modulus of cortical [cancellous] bone at a rate of between 33% (58%) (at low apparent density) and 62% (80%) (at high apparent density), are higher in the axial than in the transverse load direction. Similar results have been seen for the compressive strength of femoral bone. For the tensile and torsional strengths, Poisson's ratio and the shear modulus, only ultimate values have been found without a correlation to apparent density. For the viscoelastic behaviour of bone only data of cortical bone and in axial load direction have been described up to now. CONCLUSIONS: Anisotropic FE-models of the femur could be realized for most part with the summarized material properties of bone if characterized by apparent density and load directions. Because several mechanical properties have not been correlated to these main criteria, further experimental investigations will be necessary in future.     

A model of bone adaptation as an optimization process

The internal bone adaptation of the proximal femur is considered. A three-dimensional finite element model of the proximal femur is used. The bone remodeling in this work is numerically described by an evolutionary remodeling scheme with anisotropic material parameters and time-dependent loading. The memory of past loading is included in the model to account for the delay in the bone response from the load changes. The remodeling rate equation is derived from the structural optimization task of maximizing the stiffness in each time step. Additional information can be extracted from the optimization process; the remodeling equilibrium parameter where no apposition or resorption takes place, is defined as the element optimality conditions and the optimal design is used as an initial design for the onset of the remodeling simulation. Two examples of bone adaptation resulting from load changes are given, and the irreversible nature of the model is illustrated.     

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

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

The influence of cement mantle thickness and stem geometry on fatigue damage in two different cemented hip femoral prostheses

Experimental models can be used for pre-clinical testing of cemented and other type of hip replacements. Total hip replacement (THR) failure scenarios include, among others, cement damage accumulation and the assessment of accurate stress and strain magnitudes at the cement mantle interfaces (stem-cement and cement-bone) can be used to predict mechanical failure. The aseptic loosening scenario in cemented hip replacements is currently not fully understood, and methods of evaluating medical devices must be developed to improve clinical performance. Different results and conclusions concerning the cement micro-cracking mechanism have been reported.The aim of this study was to verify the in vitro behavior of two cemented femoral stems with respect to fatigue crack formation. Fatigue crack damage was assessed at the medial, lateral, anterior and posterior sides of the Lubinus SPII and Charnley stems. All stems were loaded and tested in stair climbing fatigue loading during one million cycles at 2 Hz. After the experiments each implanted synthetic femur was sectioned and analyzed. We observed more damage (cracks per area) for the Lubinus SPII stem, mainly on the proximal part of the cement mantle. The micro-cracking formation initiated in the stem–cement interface and grew towards the direction of cortical bone of the femur.Overall, the cement–bone interface seems to be crucial for the success of the hip replacement. The Charnley stem provoked more damage on the cement–bone interface. A failure index (maximum length of crack/maximum thickness of cement) considered was higher for the cement–stem interface of the Lubinus SPII stem. For a cement mantle thickness higher than 5 mm, cracking initiated at the cement–bone interface and depended on the opening canal process (reaming procedure and instrumentation). The analysis also showed that fatigue-induced damage on the cement mantle, increasing proximally, and depended on the axial position of the stem. The cement thickness is an important factor for the success of THR and this study evidenced that cement thickness higher than 2 mm apparently does not affect the mechanical behavior of the cement mantel and induce more crack formation on the cement–bone interface.     

The mechanical properties of trabecular bone: Dependence on anatomic location and function

In 1961, Evans and King documented the mechanical properties of trabecular bone from multiple locations in the proximal human femur. Since this time, many investigators have cataloged the distribution of trabecular bone material properties from multiple locations within the human skeleton to include femur, tibia, humerus, radius, vertebral bodies, and iliac crest. The results of these studies have revealed tremendous variations in material properties and anisotropy. These variations have been attributed to functional remodeling as dictated by Wolff's Law. Both linear and power functions have been found to explain the relationship between trabecular bone density and material properties. Recent studies have re-emphasized the need to accurately quantify trabecular bone architecture proposing several algorithms capable of determining the anisotropy, connectivity and morphology of the bone. These past studies, as well as continuing work, have significantly increased the accuracy of analytical and experimental models investigating bone, and bone/implant interfaces as well as enhanced our perspective towards understanding the factors which may influence bone formation or resorption.     

Finite element analysis of a cemented ceramic femoral component for the assembly situation in total knee arthroplasty]

The femoral components of the total knee replacements are generally made of metal. In contrast, ceramic femoral components promise improved tribological and allergological properties. However, ceramic components present a risk of failure as a result of stress peaks. Stress peaks can be minimised through adequate implant design, proper material composition and optimum force transmission between bone and implant. Thus, the quality of the implant fixation is a crucial factor. The objective of the present study was to analyse the influence of the cement layer thickness on stress states in the ceramic femoral component and in the femur. Two- and three- dimensional finite element analyses of an artificial knee joint with cement layers of different thickness and with an unbalanced cement layer thickness between the ceramic femoral component and the femur were performed. Higher stress regions occurred in the area of force transmission and in the median plane. The maximum calculated stresses were below the accepted tensile strength. Stresses were found to be lower for cement layer thickness of <2.0 mm.     

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.     

Masquelet technique: The effect of altering implant material and topography on membrane matrix composition,mechanical and barrier properties in a rat defect model

The Masquelet technique is a surgical procedure to regenerate segmental bone defects. The two-phase treatment relies on the production of a vascularized foreign-body membrane to support bone grafts over three times larger than the traditional maximum. Historically, the procedure has always utilized a bone cement spacer to evoke membrane production. However, membrane formation can easily be effected by implant surface properties such as material and topology. This study sought to determine if the membrane’s mechanical or barrier properties are affected by changing the spacer material to titanium or roughening the surface finish. Ten-week-old, male Sprague Dawley rats were given an externally stabilized, 6 mm femur defect which was filled with a pre-made spacer of bone cement (PMMA) or titanium (TI) with a smooth (∼1 μm) or roughened (∼8 μm) finish. After 4 weeks of implantation, the membranes were harvested, and the matrix composition, tensile mechanics, shrinkage, and barrier function was assessed. Roughening the spacers resulted in significantly more compliant membranes. TI spacers created membranes that inhibited solute transport more. There were no differences between groups in collagen or elastin distribution. This suggests that different membrane characteristics can be created by altering the spacer surface properties. Surgeons may unknowingly effecting membrane formation via bone cement preparation techniques.