Femoral strain changes after total hip arthroplasty--patient-specific finite element analyses 12 years after operation

Stress shielding effect after total hip arthroplasty varies between combinations of stem design and stiffness-a comparing biomechanical finite element analysis

PURPOSE

Total hip arthroplasty (THA) has become a highly frequent orthopaedic procedure. Multiple approaches have been made to design the femoral component for THA with a mechanical behaviour as close as possible to a natural femur. The aim of this study was to compare different combinations of design and biomechanical properties of THA prostheses and their impact on stress shielding of the periprosthetic bone.

METHODS

Virtual implantation of different stem designs (straight standard stem, straight short stem, anatomical short stem) by finite element analysis based on in vivo data from computer tomography was performed. For each stem, three grades of stiffness were generated, followed by a strain analysis.

RESULTS

Reduction of stem stiffness led to less stress shielding. Implantation of an anatomical short-stem prosthesis with low stiffness provided the most physiological strain-loading effect (p < 0.001).

CONCLUSION

A combination of a short and an anatomically designed stem with a low stiffness might provide a more physiological strain transfer during THA. Biomechanical properties of the femoral component for THA should be considered as a multifactorial function of dimensions, design, and stiffness.

Effect of anatomical variability on stress-shielding induced by short calcar-guided stems: Automated finite element analysis of 90 femora

Short stem hip implants are becoming increasingly popular since they preserve bone stock and presumably reduce stress-shielding. However, concerns remain whether they are suitable for a wide range of patients with varying anatomy. The aim of this study was to investigate how femoral anatomy influences stress-shielding induced by a short calcar-guided stem across a set of 90 CT-based femur models. A computational tool was developed that automatically selected the optimal size and position of the stem. Finite element models of the intact and implanted femurs were constructed and subjected to walking loads. Stress-shielding was evaluated in relevant volumes of interest of the proximal femur. After a detailed anatomical analysis, linear regression was performed to find potential correlations between anatomy and stress-shielding. Stress-shielding was found to be highest in the proximal regions on the medial and posterior side. A highly significant negative relationship was observed between stress-shielding and bone density; a strong positive relationship was observed with stem size and the valgus orientation of the stem with respect to the femur. The results reveal how anatomy influences stress-shielding, and they highlight the importance of evaluating new implant designs across a large population taking into account the anatomical variability. The study demonstrates that such large population studies can be conducted in an efficient way using an automated workflow. © 2019 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 9999:1-8, 2019.

Virtual trial to evaluate the robustness of cementless femoral stems to patient and surgical variation

Primary stability is essential for the success of cementless femoral stems. In this study, patient specific finite element (FE) models were used to assess changes in primary stability due to variability in patient anatomy, bone properties and stem alignment for two commonly used cementless femoral stems, Corail® and Summit® (DePuy Synthes, Warsaw, USA). Computed-tomography images of the femur were obtained for 8 males and 8 females. An automated algorithm was used to determine the stem position and size which minimized the endo-cortical space, and then span the plausible surgical envelope of implant positions constrained by the endo-cortical boundary. A total of 1952 models were generated and ran, each with a unique alignment scenario. Peak hip contact and muscle forces for stair climbing were scaled to the donor's body weight and applied to the model. The primary stability was assessed by comparing the implant micromotion and peri-prosthetic strains to thresholds (150 μm and 7000 µε, respectively) above which fibrous tissue differentiation and bone damage are expected to prevail. Despite the wide range of implant positions included, FE prediction were mostly below the thresholds (medians: Corail®: 20-74 µm and 1150-2884 µε, Summit®: 25-111 µm and 860-3010 µε), but sensitivity of micromotion and interfacial strains varied across femora, with the majority being sensitive (p < 0.0029) to average bone mineral density, cranio-caudal angle, post-implantation anteversion angle and lateral offset of the femur. The results confirm the relationship between implant position and primary stability was highly dependent on the patient and the stem design used.

Finite element analysis of the impact of screw insertion distal to the trochanter minor on the risk of iatrogenic subtrochanteric fracture

Iatrogenic subtrochanteric fractures are rarely encountered after cannulated screw fixation of femoral neck fractures; however, when they do occur, there can be several complications. Many orthopedic surgeons have concerns about the potential for iatrogenic subtrochanteric fractures after screw fixation distal to the trochanter minor; therefore, some surgeons are typically reluctant to perform this procedure. This study focused on the risk of an iatrogenic subtrochanteric fracture after treating femoral neck fractures with cannulated screws. The main purpose of the study was to understand iatrogenic subtrochanteric fractures and evaluate the effects on these fractures of an inverted triangular–shaped configuration for placement of the cannulated screws. A femur bone with an femoral neck fracture at a 40° incline to the horizontal plane was generated along with a representation of a three-dimensional finite element model, and three inverted triangular–shaped configurations for placement of the cannulated screws were investigated using finite element analyses. Statistical results indicated that the occurrence risk of ISF increases when the screw is located distal to the trochanter minor. Moreover, the risk of occurrence of intertrochanteric fracture increases when the screw is located medial to the trochanter minor because of local concentrated stress on the surface of the screw canals. To avoid the vulnerability of the subtrochanteric region, it was found that proximal placement of the screws using the inverted triangular–shaped configuration could yield better results. In addition, the results of this study provide suggestions on improved screw configurations.

Bone preserving level of osteotomy in short-stem total hip arthroplasty does not influence stress shielding dimensions - a comparing finite elements analysis

Background

The main objective of every new development in total hip arthroplasty (THA) is the longest possible survival of the implant. Periprosthetic stress shielding is a scientifically proven phenomenon which leads to inadvertent bone loss. So far, many studies have analysed whether implanting different hip stem prostheses result in significant preservation of bone stock. The aim of this preclinical study was to investigate design-depended differences of the stress shielding effect after implantation of a selection of short-stem THA-prostheses that are currently available.

Methods

Based on computerised tomography (CT), a finite elements (FE) model was generated and a virtual THA was performed with different stem designs of the implant. Stems were chosen by osteotomy level at the femoral neck (collum, partial collum, trochanter sparing, trochanter harming). Analyses were performed with previously validated FE models to identify changes in the strain energy density (SED).

Results

In the trochanteric region, only the collum-type stem demonstrated a biomechanical behaviour similar to the native femur. In contrast, no difference in biomechanical behaviour was found between partial collum, trochanter harming and trochanter sparing models. All of the short stem-prostheses showed lower stress-shielding than a standard stem.

Conclusion

Based on the results of this study, we cannot confirm that the design of current short stem THA-implants leads to a different stress shielding effect with regard to the level of osteotomy. Somehow unexpected, we found a bone stock protection in metadiaphyseal bone by simulating a more distal approach for osteotomy. Further clinical and biomechanical research including long-term results is needed to understand the influence of short-stem THA on bone remodelling and to find the optimal stem-design for a reduction of the stress shielding effect.

Electronic supplementary material

The online version of this article (doi:10.1186/s12891-017-1702-2) contains supplementary material, which is available to authorized users.

In vitro implant-bone interface pressure measurements for a cementless femoral implant. A preliminary study

PURPOSE

Implants endurance as well as a good clinical tolerance depends on the recovery of a physiological stress distribution within bone after implantation. The purpose of the present work was to develop an alternative technique using Force Sensing Resistors (FSR) to gather in vitro pressure values at the implant-bone interface for a cementless implant.

METHOD

Eight cementless femoral stems were instrumented with six calibrated FSR bonded on each facet and then implanted in eight cadaver femurs. Compression tests were performed until failure and FSR pressure values were recorded.

RESULTS

The average failure load was 4241 N. The maximum contact pressure measured with the FSR averaged 1.965 MPa.

CONCLUSION

FSR reached many of the requirements for an ideal implant-bone interfacial sensor. This experimentation provided in vitro quantitative data on contact pressure at the implant-bone interface, which could help understanding stress shielding phenomenon and developing relevant numerical model.

Long-term implant fixation and stress-shielding in total hip replacement

Implant fixation implies a strong and durable mechanical bond between the prosthetic component and host skeleton. Assuming the short-term impediments to implant fixation are successfully addressed and that longer-term issues such as late infection and mechanical failure of the components are avoided, the biological response of the host tissue to the presence of the implant is critical to long-term success. In particular, maintenance of adequate peri-prosthetic bone stock is a key factor. Two major causes of bone loss in the supporting bone are adverse bone remodeling in response to debris shed from the implant and stress-shielding. Here, I review some of the major lessons learned from studying stress-shielding-induced bone loss. It is well known that stress-shielding can be manipulated by altering implant design, but less well appreciated that the development of bone anabolic agents may make it possible to reduce the severity of stress-shielding and the associated bone loss by augmenting the host skeleton through the use of locally or systemically delivered agents. In most cases, mechanical, material and biological factors do not act in isolation, emphasizing that it is often not possible to optimize all boundary conditions.

Four decades of finite element analysis of orthopaedic devices: where are we now and what are the opportunities?

Finite element has been used for more than four decades to study and evaluate the mechanical behaviour total joint replacements. In Huiskes seminal paper "Failed innovation in total hip replacement: diagnosis and proposals for a cure", finite element modelling was one of the potential cures to avoid poorly performing designs reaching the market place. The size and sophistication of models has increased significantly since that paper and a range of techniques are available from predicting the initial mechanical environment through to advanced adaptive simulations including bone adaptation, tissue differentiation, damage accumulation and wear. However, are we any closer to FE becoming an effective screening tool for new devices? This review contains a critical analysis of currently available finite element modelling techniques including (i) development of the basic model, the application of appropriate material properties, loading and boundary conditions, (ii) describing the initial mechanical environment of the bone-implant system, (iii) capturing the time dependent behaviour in adaptive simulations, (iv) the design and implementation of computer based experiments and (v) determining suitable performance metrics. The development of the underlying tools and techniques appears to have plateaued and further advances appear to be limited either by a lack of data to populate the models or the need to better understand the fundamentals of the mechanical and biological processes. There has been progress in the design of computer based experiments. Historically, FE has been used in a similar way to in vitro tests, by running only a limited set of analyses, typically of a single bone segment or joint under idealised conditions. The power of finite element is the ability to run multiple simulations and explore the performance of a device under a variety of conditions. There has been increasing usage of design of experiments, probabilistic techniques and more recently population based modelling to account for patient and surgical variability. In order to have effective screening methods, we need to continue to develop these approaches to examine the behaviour and performance of total joint replacements and benchmark them for devices with known clinical performance. Finite element will increasingly be used in the design, development and pre-clinical testing of total joint replacements. However, simulations must include holistic, closely corroborated, multi-domain analyses which account for real world variability.

Three-dimensional finite element analysis and comparison of a new intramedullary fixation with interlocking intramedullary nail

This study was set to introduce a new intramedullary fixation, explore its biomechanical properties, and provide guidance for further biomechanical experiments. With the help of CT scans and finite element modeling software, finite element model was established for a new intramedullary fixation and intramedullary nailing of femoral shaft fractures in a volunteer adult. By finite element analysis software ANSYS 10.0, we conducted 235-2,100 N axial load, 200-1,000 N bending loads and 2-15 Nm torsional loading, respectively, and analyzed maximum stress distribution, size, and displacement of the fracture fragments of the femur and intramedullary nail. During the loading process, the maximum stress of our new intramedullary fixation were within the normal range, and the displacement of the fracture fragments was less than 1 mm. Our new intramedullary fixation exhibited mechanical reliability and unique advantages of anti-rotation, which provides effective supports during fracture recovery.

Biomechanical assessment of composite versus metallic intramedullary nailing system in femoral shaft fractures: A finite element study

BACKGROUND

Intramedullary nails are the primary choice for treating long bone fractures. However, complications following nail surgery including non-union, delayed union, and fracture of the bone or the implant still exist. Reducing nail stiffness while still maintaining sufficient stability seems to be the ideal solution to overcome the abovementioned complications.

METHODS

In this study, a new hybrid concept for nails made of carbon fibers/flax/epoxy was developed in order to reduce stress shielding. The mechanical performance of this new implant in terms of fracture stability and load sharing was assessed using a comprehensive non-linear FE model. This model considers several mechanical factors in nine fracture configurations at immediately post-operative, and in the healed bone stages.

RESULTS

Post-operative results showed that the hybrid composite nail increases the average normal force at the fracture site by 319.23N (P<0.05), and the mean stress in the vicinity of fracture by 2.11MPa (P<0.05) at 45% gait cycle, while only 0.33mm and 0.39mm (P<0.05) increases in the fracture opening and the fragments' shear movement were observed. The healed bone results revealed that implantation of the titanium nail caused 20.2% reduction in bone stiffness, while the composite nail lowered the stiffness by 11.8% as compared to an intact femur.

INTERPRETATION

Our results suggest that the composite nail can provide a preferred mechanical environment for healing, particularly in transverse shaft fractures. This may help bioengineers better understand the biomechanics of fracture healing, and aid in the design of effective implants.

A fatigue loading model for investigation of iatrogenic subtrochanteric fractures of the femur

BACKGROUND

Biomechanics of iatrogenic subtrochanteric femur fractures have been examined. Previously-described loading models employed monotonic loading on the femoral head, which is limited in emulating physiological features. We hypothesize that cyclic loading combined with the engagement of abductor forces will reliably cause iatrogenic subtrochanteric fractures.

METHODS

Finite element analysis determined the effects of adding the abductor muscle forces to the hip contact force around holes located in the lateral femoral cortex. Finite element analysis predictions were validated by strain gage measurements using Sawbones™ femurs (Pacific Research Laboratories, Inc., Vashon, Washington, USA) with or without abductor muscle forces. The newly developed physiologically-relevant loading model was tested on cadaveric femurs (N=8) under cyclic loading until failure.

FINDINGS

Finite element analysis showed the addition of the abductor muscle forces increased the maximum surface cortical strain by 107% and the strain energy density by 332% at the lateral femoral cortex. Strain gages detected a 72.9% increase in lateral cortical strain using the combined loading model. The cyclic, combined loading led to subtrochanteric fractures through the drill hole in all cadaveric femurs.

INTERPRETATION

Finite element analysis simulations, strain gage measurements, and cyclic loading of fresh-frozen femurs indicate the inclusion of abductor forces increases the stress and strain at the proximal-lateral femoral cortex. Furthermore, a cyclic loading model that incorporates a hip contact force and abductor muscles force creates the clinically encountered subtrochanteric fractures in vitro. This physiologically-relevant loading model may be used to further study iatrogenic subtrochanteric femur fractures.

Predicting Bone Remodeling in Response to Total Hip Arthroplasty: Computational Study Using Mechanobiochemical Model

Periprosthetic bone loss following total hip arthroplasty (THA) is a serious concern leading to the premature failure of prosthetic implant. Therefore, investigating bone remodeling in response to hip arthroplasty is of paramount for the purpose of designing long lasting prostheses. In this study, a thermodynamic-based theory, which considers the coupling between the mechanical loading and biochemical affinity as stimulus for bone formation and resorption, was used to simulate the femoral density change in response to THA. The results of the numerical simulations using 3D finite element analysis revealed that in Gruen zone 7, after remarkable postoperative bone loss, the bone density started recovering and got stabilized after 9% increase. The most significant periprosthetic bone loss was found in Gruen zone 7 (−17.93%) followed by zone 1 (−13.77%). Conversely, in zone 4, bone densification was observed (+4.63%). The results have also shown that the bone density loss in the posterior region of the proximal metaphysis was greater than that in the anterior side. This study provided a quantitative figure for monitoring the distribution variation of density throughout the femoral bone. The predicted bone density distribution before and after THA agree well with the bone morphology and previous results from the literature.

Biomechanical analysis using infrared thermography of a traditional metal plate versus a carbon fibre/epoxy plate for Vancouver B1 femur fractures

Traditional high-stiffness metal plates for Vancouver B1 femur shaft fractures below the tip of a hip implant can cause stress shielding, bone resorption, and implant loosening. This is the first study to compare the biomechanics of a traditional metal plate versus a low-stiffness carbon fibre/epoxy composite plate for this injury. A total hip replacement was implanted in two previously validated intact artificial femurs. Femurs were fitted with either a metal or composite plate and had a 5 mm fracture gap created to simulate a Vancouver B1 shaft fracture. Femurs were cyclically loaded using 5 Hz at 7° of adduction with an average axial load of 800 N (range = 400–1200 N). Overall mechanical stiffnesses and femur and plate thermographic stresses were obtained. Femur/metal plate stiffness (698 N/mm) was only 12% higher than femur/composite plate stiffness (625 N/mm). The femur with the composite plate had 22.7% higher combined average stress compared to the femur with the metal plate, having specific differences of 29.5% (anterior view), 33.9% (posterior view), 1.0% (medial view), and 26.4% (lateral view). The composite plate itself had an average 21.1% reduction in stress compared to the metal plate. The composite plate reduced stress shielding, yet provided adequate stiffness.

Patient-specific modelling of bone and bone-implant systems: the challenges

In the past three decades, finite element (FE) modelling has provided considerable understanding to the area of musculoskeletal biomechanics. However, most of this understanding has been generated using generic, standardised or idealised models. Patient-specific modelling (PSM) is almost never used for making clinical decisions. Imaging technologies have made it possible to create patient-specific geometries and FE meshes for modelling. While these have brought us closer to PSM, several challenges associated with the definition of material properties, loads, boundary conditions and interaction between components still need to be overcome. This study reviews the current status of PSM with respect to defining material behaviour and prescribing boundary conditions and interactions. With regard to the constitutive modelling of bone, it is seen that imaging is being increasingly used to define elastic properties (isotropic as well as anisotropic). However, the post-elastic and time-dependent behaviour, important for several modelling situations, is mostly obtained from in vitro experiments. Strain-based plasticity, not commonly available in FE codes, appears to have the potential of reducing an element of patient-specificity in modelling the yielding behaviour of bone. PSM of real boundary conditions that include muscles and ligaments continues to remain a challenge; many clinically relevant questions can be, however, answered without their inclusion. Simulation techniques to undertake PSM of interactions between bone and uncemented implants are available. Interference fit employed in both joint replacement fracture treatments induces considerable preload whose inclusion in models is important for the prediction of interface behaviour.

Accounting for patient variability in finite element analysis of the intact and implanted hip and knee: a review

It is becoming increasingly difficult to differentiate the performance of new joint replacement designs using available preclinical test methods. Finite element analysis is commonly used and the majority of published studies are performed on representative anatomy, assuming optimal implant placement, subjected to idealised loading conditions. There are significant differences between patients and accounting for this variability will lead to better assessment of the risk of failure. This review paper provides a comprehensive overview of the techniques available to account for patient variability. There is a brief overview of patient-specific model generation techniques, followed by a review of multisubject patient-specific studies performed on the intact and implanted femur and tibia. In particular, the challenges and limitations of manually generating models for such studies are discussed. To efficiently account for patient variability, the application of statistical shape and intensity models (SSIM) are being developed. Such models have the potential to synthetically generate thousands of representative models generated from a much smaller training set. Combined with the automation of the prosthesis implantation process, SSIM provides a potentially powerful tool for assessing the next generation of implant designs. The potential application of SSIM are discussed along with their limitations.

Influence of femur size and morphology on load transfer in the resurfaced femoral head: A large scale, multi-subject finite element study

Femoral resurfacing has become an increasingly popular procedure, especially for young, active patients. The procedure is known to alter load transfer through the proximal femur and this has been linked with the most commonly observed complication, neck fracture. An intriguing observation noted by registry data and clinical studies is an inverse relationship between implant size and revision rate. While computational analysis has become an established part of biomedical engineering, the majority of work uses a single or small set of bone models, with a single implant size, due to the constraints of time and data availability. Therefore, it has been infeasible to run a study incorporating natural inter-patient variability or the performance of smaller implants could not be meaningfully studied. In previous work a statistical model of the whole femur was used to generate large numbers of unique, realistic, FE-ready femur models describing both geometry and material properties. The current study demonstrates a methodology for virtually implanting and performing stress analysis of cemented femoral resurfacing components, with model specific sizing and orientation. Automated analysis of 400 generated femurs, in both implanted and intact configurations showed the strain changes induced by resurfacing. This produced a statistically meaningful number of results and allowed the examination of outliers. Results showed increased femoral neck strain changes potentially increasing the risk of neck fracture, associated with smaller, less dense femurs and smaller implant sizes; agreeing with clinical observations. The study demonstrates a methodology for more comprehensive analyses, based on populations rather than individuals.

Quantitative CT with finite element analysis: towards a predictive tool for bone remodelling around an uncemented tapered stem

PURPOSE

We used quantitative CT in conjunction with finite element analysis to provide a new tool for assessment of bone quality after total hip arthroplasty in vivo. The hypothesis of this prospective five-year study is that the combination of the two modalities allows 3D patient-specific imaging of cortical and cancellous bone changes and stress shielding.

METHOD

We tested quantitative CT in conjunction with finite elements on a cohort of 29 patients (31 hips) who have been scanned postoperatively and at one year, two years and five years follow-up. The method uses cubic Hermite finite element interpolation for efficient mesh generation directly from qCT datasets. The element Gauss points that are used for the geometric interpolation functions are also used for interpolation of osteodensitometry data.

RESULTS

The study showed changes of bone density suggestive of proximal femur diaphysis load transfer with osteointegration and moderate metaphyseal stress shielding. Our model revealed that cortical bone initially became porous in the greater trochanter, but this phenomenon progressed to the cortex of the lesser trochanter and the posterior aspect of the metaphysis. The diaphyseal area did not experience major change in bone density for either cortical or cancellous bone.

CONCLUSION

The combination of quantitative CT with finite element analysis allows visualization of changes to bone density and architecture. It also provides correlation of bone density/architectural changes with stress patterns enabling the study of the effects of stress shielding on bone remodelling in vivo. This technology can be useful in predicting bone remodeling and the quality of implant fixation using prostheses with different design and/or biomaterials.

Femoral head density on CT scans of patients following hip fracture fixation by expandable proximal peg or dynamic screw

Computed tomography (CT) is currently considered to be an accurate method for evaluating bone density. We evaluated the CT measurements of bone density using the Hounsfield units (HUs) in 23 patients who had been operated in the past for an extra-capsular hip fracture. Twelve patients were treated with a dynamic hip screw and 11 with a proximal femoral expandable hip nail. All the CTs had been performed for non-orthopedic purposes. Bone density with a region of interest (ROI) could be assessed for both hips. We compared the bone density between the operated versus the non-operated sides as well as between the two surgical groups. Bone density was higher in the hip peg (the femoral component of the expandable nail) side 262.5 (range, 169-351) HU, compared to the opposite non-operated side and to the hip screw group 194 (range, 99-283) HU. The hip screw side had decreased bone density compared to the opposite non-operated side. We were able to define a density index and a difference index: both were higher in the hip peg group. These findings persisted over time. It would be interesting to speculate that increased bone density around an expandable peg provides better fracture stabilization and probably faster healing than a dynamic hip screw.

Bone remodeling in a new biomimetic polymer-composite hip stem

Adaptive bone remodeling is an important factor that leads to bone resorption in the surrounding femoral bone and implant loosening. Taking into account this factor in the design of hip implants is of clinical importance, because it allows the prediction of the bone-density redistribution and enables the monitoring of bone adaptation after prosthetic implantation. In this article, adaptive bone remodeling around a new biomimetic polymer-composite-based (CF/PA12) hip prosthesis is investigated to evaluate the amount of stress shielding and bone resorption. The design concept of this new prosthesis is based on a hollow substructure made of hydroxyapatite-coated, continuous carbon fiber (CF)-reinforced polyamide 12 (PA12) composite with an internal soft polymer-based core. Strain energy density theory coupled with 3D Finite Element models is used to predict bone density redistributions in the femoral bone before and after total hip replacement (THR) using both polymer-composite and titanium (Ti) stems. The result of numerical simulations of bone remodeling revealed that the CF/PA12 composite stem generates a better bone density pattern compared with the Ti-based stem, indicating the effectiveness of the composite stem to reduce bone resorption caused by stress-shielding phenomenon. This may result in an extended lifetime of THR.

A fast and robust patient specific Finite Element mesh registration technique: application to 60 clinical cases

Finite Element mesh generation remains an important issue for patient specific biomechanical modeling. While some techniques make automatic mesh generation possible, in most cases, manual mesh generation is preferred for better control over the sub-domain representation, element type, layout and refinement that it provides. Yet, this option is time consuming and not suited for intraoperative situations where model generation and computation time is critical. To overcome this problem we propose a fast and automatic mesh generation technique based on the elastic registration of a generic mesh to the specific target organ in conjunction with element regularity and quality correction. This Mesh-Match-and-Repair (MMRep) approach combines control over the mesh structure along with fast and robust meshing capabilities, even in situations where only partial organ geometry is available. The technique was successfully tested on a database of 5 pre-operatively acquired complete femora CT scans, 5 femoral heads partially digitized at intraoperative stage, and 50 CT volumes of patients' heads. In the latter case, both skin and bone surfaces were taken into account by the mesh registration process in order to model the face muscles and fat layers. The MMRep algorithm succeeded in all 60 cases, yielding for each patient a hex-dominant, Atlas based, Finite Element mesh with submillimetric surface representation accuracy, directly exploitable within a commercial FE software.

Minimal stress shielding with a Mallory-Head titanium femoral stem with proximal porous coating in total hip arthroplasty

BACKGROUND

As longevity of cementless femoral components enters the third decade, concerns arise with long-term effects of fixation mode on femoral bone morphology. We examined the long-term consequences on femoral remodeling following total hip arthroplasty with a porous plasma-sprayed tapered titanium stem.

METHODS

Clinical data and radiographs were reviewed from a single center for 97 randomly selected cases implanted with the Mallory-Head Porous femoral component during primary total hip arthroplasty. Measurements were taken from preoperative and long-term follow-up radiographs averaging 14 years postoperative. Average changes in the proximal, middle and diaphyseal zones were determined.

RESULTS

On anteroposterior radiographs, the proximal cortical thickness was unchanged medially and the lateral zone increased 1.3%. Middle cortical thickness increased 4.3% medially and 1.2% laterally. Distal cortical thickness increased 9.6% medially and 1.9% laterally. Using the anteroposterior radiographs, canal fill at 100 mm did not correlate with bony changes at any level (Spearman's rank correlation coefficient of -0.18, 0.05, and 0.00; p value = 0.09, 0.67, 0.97). On lateral radiographs, the proximal cortical thickness increased 1.5% medially and 0.98% laterally. Middle cortical thickness increased 2.4% medially and 1.3% laterally. Distal cortical thickness increased 3.5% medially and 2.1% laterally. From lateral radiographs, canal fill at 100 mm correlated with bony hypertrophy at the proximal, mid-level, and distal femur (Spearman's rank correlation coefficient of 0.85, 0.33, and 0.28, respectively; p value = 0.001, 0.016, and 0.01, respectively).

CONCLUSION

Stress shielding is minimized with the Mallory-Head titanium tapered femoral stem with circumferential proximal plasma-sprayed coating in well-fixed and well-functioning total hip arthroplasty. Additionally, the majority of femora demonstrated increased cortical thickness in all zones around the stem prosthesis.

LEVEL OF EVIDENCE

Therapeutic Level III.

Computational modeling for bedside application

Advances in computer power, novel diagnostic and therapeutic medical technologies, and an increasing knowledge of pathophysiology from gene to organ systems make it increasingly feasible to apply multiscale patient-specific modeling based on proven disease mechanisms. Such models may guide and predict the response to therapy in many areas of medicine. This is an exciting and relatively new approach, for which efficient methods and computational tools are of the utmost importance. Investigators have designed patient-specific models in almost all areas of human physiology. Not only will these models be useful in clinical settings to predict and optimize the outcome from surgery and non-interventional therapy, but they will also provide pathophysiologic insights from the cellular level to the organ system level. Models, therefore, will provide insight as to why specific interventions succeed or fail.

Development and Validation of Patient-Specific Finite Element Models of the Hemipelvis Generated From a Sparse CT Data Set

To produce a patient-specific finite element (FE) model of a bone such as the pelvis, a complete computer tomographic (CT) or magnetic resonance imaging (MRI) geometric data set is desirable. However, most patient data are limited to a specific region of interest such as the acetabulum. We have overcome this problem by providing a hybrid method that is capable of generating accurate FE models from sparse patient data sets. In this paper, we have validated our technique with mechanical experiments. Three cadaveric embalmed pelves were strain gauged and used in mechanical experiments. FE models were generated from the CT scans of the pelves. Material properties for cancellous bone were obtained from the CT scans and assigned to the FE mesh using a spatially varying field embedded inside the mesh while other materials used in the model were obtained from the literature. Although our FE meshes have large elements, the spatially varying field allowed them to have location dependent inhomogeneous material properties. For each pelvis, five different FE meshes with a varying number of patient CT slices (8–12) were generated to determine how many patient CT slices are needed for good accuracy. All five mesh types showed good agreement between the model and experimental strains. Meshes generated with incomplete data sets showed very similar stress distributions to those obtained from the FE mesh generated with complete data sets. Our modeling approach provides an important step in advancing the application of FE models from the research environment to the clinical setting.

Numerical investigations of stress shielding in total hip prostheses

Aseptic loosening of the prosthesis is still a problem in artificial joint implants. The loosening can be caused by, among other factors, resorption of the bone surrounding the prosthesis owing to stress shielding. In order to find out the influence of the prosthesis type on post-operative stress shielding, a static finite element analysis of a femur provided with the conventional uncemented stem BICONTACT and of one with the femoral neck prosthesis SPIRON was carried out. Strain energy densities and maximal principal strain distributions were calculated and compared with the physiological situation. Here, stress shielding was demonstrated in both periprosthetic femora. To determine the areas of the stress shielding, the bone in each FE model was subdivided into three regions of interest (ROI): proximal, diaphyseal, and distal. The numerical computations show stress shielding in the proximal ROI of both periprosthetic femora. Diaphyseally, the femoral neck prosthesis SPIRON, in contrast to the conventional uncemented long-stem prosthesis BICONTACT, causes no decrease in the strain distribution and thus no stress shielding. Distally, no change in the load distribution of either periprosthetic femur could be found, compared with the physiological situation.

Determination of suitable sample sizes for multi-patient based finite element studies

Finite element analysis is used extensively to assess joint replacements, but the majority of these are single sample studies. Recent investigations have suggested that such studies are unable to account for natural inter-patient variation in bone geometry and material property distribution. Recent developments in computer tomography based analyses make multiple sample studies possible; the question remains how many femurs are required to perform a study which accounts for such variations. This work investigates the factors that should be considered in answering this question. It explores sample sizing techniques when comparing strain distribution in the intact and implanted femur and when comparing two or more implant designs in a group of femurs. An example analysis of the effect of femoral head resurfacing was undertaken. Two sample sizing calculations were utilised, one based on achieving the desired precision in results, the other based on determining if a significant difference exists between two designs. The analysis shows that reasonable statistical precision can be achieved with a group of femurs. The study was also able to determine a suitable sample size for the analysis of a statistically significant difference between two groups of femurs with varying design parameters. The study concluded that while sample sizing is recommended for an accurate analysis, consideration must be made for the practicality of such a task.

Influence of changes in stem positioning on femoral loading after THR using a short-stemmed hip implant

BACKGROUND

Short-stemmed hip implants were introduced to conserve proximal bone mass and may facilitate the use of minimally invasive surgery, in which smaller incisions limit access to the joint. This limited access may increase the risk of surgical mal-positioning of the implant, however the sensitivity of femoral loading to such mal-positioning of a short-stemmed implant has not been studied.

METHODS

Finite element models were developed of a femur and a short-stemmed implant positioned to reproduce the intact hip centre, as well as with the implant placed in increased anteversion or offset. The effect of these surgical variables on femoral loading was examined for walking and stair climbing using loads from a validated musculoskeletal model. Results of the implanted models were compared with an intact model to evaluate stress shielding.

FINDINGS

Implant position had little influence on cortical strains along the length of the diaphysis, although strains decreased by up to 95% at the neck resection level compared to the intact femur. In the proximal Gruen zones I and VII strain energy density among the implanted models varied by up to 0.4 kJ/m(3) (28%) and 0.6 kJ/m(3) (24%) under walking and stair climbing, respectively. All implanted models showed characteristic proximal stress shielding, indicated by a decrease in strain energy density of up to 5.4 kJ/m(3) (69%) compared to the intact femur.

INTERPRETATION

Small changes in stem placement would likely have little influence on the internal loading of the femur after bone ingrowth has been achieved, however a reduction in strain energy density and therefore stress shielding was seen even for a short-stemmed implant, which may have consequences for longer-term bone remodelling.

Physiologically based boundary conditions in finite element modelling

Finite element analysis has been used extensively in the study of bone loading and implant performance, such as in the femur. The boundary conditions applied vary widely, generally producing excessive femoral deformation, and although it has been shown that the muscle forces influence femoral deflections and loading, little consideration has been given to the displacement constraints. It is hypothesised that careful application of physiologically based constraints can produce physiological deformation, and therefore straining, of the femur. Joint contact forces and a complete set of muscle forces were calculated based on the geometry of the Standardised Femur using previously validated musculoskeletal models. Five boundary condition cases were applied to a finite element model of the Standardised Femur: (A) diaphyseally constrained with hip contact and abductor forces; (B) case A plus vasti forces; (C) case A with complete set of muscle forces; (D) distally constrained with all muscle forces; (E) physiological constraints with all muscle forces. It was seen that only the physiological boundary conditions, case E, produced physiological deflections (< 2.0mm) of the femoral head in both the coronal and sagittal planes, which resulted in minimal reaction forces at the constrained nodes. Strains in the mid-diaphysis varied by up to 600 micro-strain under walking loads and 1000 micro-strain under stair climbing loads. The mode of loading, as indicated by the strain profiles on the cortex also varied substantially under these boundary conditions, which has important consequences for studies that examine localised bone loading such as fracture or bone remodelling simulations.