Predicting the subject-specific primary stability of cementless implants during pre-operative planning: Preliminary validation of subject-specific finite-element models

Does preclinical analysis based on static loading underestimate post-surgery stem micromotion in THA as opposed to dynamic gait loading?

The success of cementless hip stems depends on the primary stability of the implant quantified by the amount of micromotion at the bone-stem interface. Most finite element (FE)-based preclinical studies on post-surgery stem stability rely on static analysis. Hence, the effect of dynamic gait loading on bone-stem relative micromotion remains virtually unexplored. Furthermore, there is a paucity of research on the primary stability of grooved stems as opposed to plain stem design. The primary aim of this FE study was to understand whether transient dynamic gait had any incremental effect on the net micromotion results and to further draw insights into the effects of grooved texture vis-à-vis a plain model on micromotion and proximal load transfer in host bone. Two musculoskeletal loading regimes corresponding to normal walking (NW) and stair climbing (SC) were considered. Although marginally improved load transfer was predicted proximally for the grooved construct under static loading, the micromotion values (max: NW ~ 7 μm; SC ~ 10 μm) were found to be considerably less in comparison to plain stem (max: NW ~ 50 μm; SC ~ 20 μm). For both physiological load cases, a significant surge in micromotion values was predicted in dynamic analyses as opposed to static analyses for the grooved stem (~ 390% greater). For the plain model, the increase in these values from static to dynamic loading is relatively moderate yet clinically significant (~ 230% greater). This suggests that the qualitative similarities notwithstanding, there were significant dissimilarities in the quantitative trends of micromotion for different cases under both analyses.

The influence of loading configurations on numerical evaluation of failure mechanisms in an uncemented femoral prosthesis

The clinical relevance of numerical predictions of failure mechanisms in femoral prosthesis could be impaired due to simplification of musculoskeletal loading. This study investigated the extent to which loading configurations affect the preclinical analysis of an uncemented femoral implant. Patient-specific, CT-scan based FE models of intact and implanted femurs were developed and analysed using three loading configurations, which comprised of load cases representing daily activities. First loading configuration consisted of two load cases, each of walking and stair climbing. The second consisted of more number of load cases for each of these activities. The third included load cases of additional activities of standing up and sitting down. Failure criteria included maximum principal strains, interface debonding, implant-bone relative displacement and adaptive bone remodelling. Simplified loading configurations led to a reduction (100-1500 με) around cortical principal strains. The area prone to interface debonding were observed in the proximo-medial part of implant and was maximum when all activities were considered. This area was reduced by 35%, when simplified loading configurations were chosen. Interfacial area of 88%-96% experienced implant-bone relative displacements below 40 μm; however maximum of 110 μm was observed at the calcar region. Lack of consideration of variety of activities overestimated (30%-50%) bone resorption around the lateral part of the implant; hence, these bone remodelling results were less clinically relevant. Considering a variety daily activities along with an adequate number of load cases for each activity seemed necessary for pre-clinical evaluations of reconstructed femur.

COMPARATIVE ANALYSIS FOR THREE FIXTURES OF PAUWELLS-II BY THE BIOMECHANICAL FINITE ELEMENT METHOD

Little is known about why and how biomechanics govern the hypothesis that three-Lag-Screw (3LS) fixation is a preferred therapeutic technique. A series models of surgical internal-fixation for femoral neck fractures of Pauwells-II will be constructed by an innovative approach of finite element so as to determine the most stable fixation by comparison of their biomechanical performance. Seventeen sets of CT scanned femora were imported onto Mimics extracting 3D models; these specimens were transferred to Geomagic Studio for a simulative osteotomy and kyrtograph; then, they underwent UG to fit simulative solid models; three sorts of internal fixators were expressed virtually by Pro-Engineer. Processed by Hypermesh, all compartments were assembled onto three systems actually as “Dynamic hip screw (DHS), 3LS and DHS+LS”. Eventually, numerical models of Finite Elemental Analysis (FEA) were exported to AnSys for solution. Three models for fixtures of Pauwells-II were established, validated and analyzed with the following findings: Femoral-shaft stress for [Formula: see text](3LS) is the least; Internal-fixator stress (MPa) for [Formula: see text]; Integral stress (MPa) for [Formula: see text]; displacement of femoral head (mm) for a[Formula: see text](DHS+LS) = 0.735; displacement of femoral shaft (mm) for [Formula: see text]; and displacement of fixators for [Formula: see text]. Mechanical comparisons for other femoral parks are insignificantly different, and these data can be abstracted as follows: the stress of 3LS-system was checked to be the least, and an interfragmentary displacement of DHS+LS assemblages was assessed to be the least”. A 3LS-system should be recommended to clinically optimize a Pauwells-II facture; if treated by this therapeutic fixation, breakage of fixators or secondary fracture is supposed to occur rarely. The strength of this study is that it was performed by a computer-aided simulation, allowing for design of a preoperative strategy that could provide acute correction and decrease procedure time, without harming to humans or animals.

A comparative assessment of two designs of hip stem using rule-based simulation of combined osseointegration and remodelling

The stem–bone interface of cementless total hip arthroplasty undergoes an adaptive process of bone ingrowth until the two parts become osseointegrated. Another important phenomenon associated with aseptic loosening of hip stem is stress-shielding induced adverse bone remodelling. The objective of this study was to preclinically assess the relative performances of two distinct designs of hip stems by addressing the combined effect of bone remodelling and osseointegration, based on certain rule-based criteria obtained from the literature. Premised upon non-linear finite element analyses of patient-specific implanted femur models, the study attempts to ascertain in silico outcome of the hip stem designs based on an evolutionary interfacial condition, and to further comment on the efficacy of the rule-based technique on the prediction of peri-prosthetic osseointegration. One of the two hip stem models was a trade-off design obtained from an earlier shape optimization study, and the other was based on TriLock stem (DePuy). Both designs predicted similar long-term osseointegration (∼89% surface), although trade-off stem predicted higher post-operative osseointegration. Proximal bone resorption was found higher for TriLock (by ∼110%) as compared to trade-off model. The rule-based technique predicted clinically coherent osseointegration around both stems and appears to be an alternative to expensive mechanobiology-based schemes.

Comparative analysis for five fixations of Pauwels-I by the biomechanical finite-element method

Background: Little is known about how biomechanics govern the five fixtures such as DHS, MLS, DHS + LS, LP, and HA are accepted as common therapeutic techniques. Aims and objectives: A series of numerical models for a femoral neck fracture of Pauwels-I will be constructed by innovative approach of finite element in order to determine the most optimized option in comparison with biomechanical performance. Method: Twenty sets of computer tomography scanned femora were imported onto Mimics to extract 3 D models; these specimens were transferred to Geomagic-Studio for a simulative osteotomy and kyrtograph; then, they underwent UG to fit simulative solid models; 5 sorts of fixture were then expressed by Pro-Engineer virtually. After processing with HyperMesh, all compartments (fracture model + internal implant) were assembled onto 5 systems: "Dynamic Hip Screw (DHS), Multiple Lag screw (MLS), DHS + LS, femoral Locking Plate (LP) and HemiArthroplasty (HA)." Eventually, numerical models of the finite-elemental analysis were exported to AnSys to determine the solution. Result: Four models of fixation and a simulation of HA for Pauwels-I were established, validated, and analyzed with the following findings: In term of displacement, these 5 fixtures ranged between 0.3801 and 0.7536 mm have no significant difference; in term of stress, the averages of peaks for integral assemblage are b(MLS) = 43.5766 ≈< d(LP) = 43.6657 ≈< e(Ha) = 43.6657 < c(DHS + LS) = 66.5494 < a(DHS) = 105.617 in MPa indicate that MLS, LP and HA are not significantly different, but less than DHS + LS or DHS in each. Conclusion: A fixture of MLS or LP with optional HA should be recommended to clinically optimize a Pauwels-I facture.

Modal analysis for the assessment of cementless hip stem primary stability in preoperative THA planning

This numerical vibration finite element (FE) study introduces resonance three-dimensional planning (RP3D) to assess preoperatively the primary stability of a cementless stem for total hip arthroplasty. Based on a patient's CT-scan and a numerical model of a stem, RP3D aims at providing mechanical criteria indicative of the achievable primary stability. We investigate variations of the modal response of the stem to changes of area and apparent stiffness of the bone-implant interface. The model is computationally cheap as it does not include a mesh of the bone. The apparent stiffness of the bone is modeled by springs attached to the nodes of the stem's mesh. We investigate an extended range of stiffness values while, in future works, patient's specific Hounsfield values could be used to define stiffness. We report modal frequencies, shapes, and a ratio of elastic potential energies (rEPE) that quantifies the proximal motion that should be minimum for a stable stem. The modal response exhibits a clear transition between loose and tight contact as area and stiffness of the interface increase. rEPE thresholds that could potentially discriminate preoperatively between stable and unstable stems are given for a Symbios SPS^® size C stem.

Finite Element Analysis of porously punched prosthetic short stem virtually designed for simulative uncemented Hip Arthroplasty

Background

There is no universal hip implant suitably fills all femoral types, whether prostheses of porous short-stem suitable for Hip Arthroplasty is to be measured scientifically.

Methods

Ten specimens of femurs scanned by CT were input onto Mimics to rebuild 3D models; their *stl format dataset were imported into Geomagic-Studio for simulative osteotomy; the generated *.igs dataset were interacted by UG to fit solid models; the prosthesis were obtained by the same way from patients, and bored by punching bears designed by Pro-E virtually; cements between femora and prosthesis were extracted by deleting prosthesis; in HyperMesh, all compartments were assembled onto four artificial joint style as: (a) cemented long-stem prosthesis; (b) porous long-stem prosthesis; (c) cemented short-stem prosthesis; (d) porous short-stem prosthesis. Then, these numerical models of Finite Element Analysis were exported to AnSys for numerical solution.

Results

Observed whatever from femur or prosthesis or combinational femora-prostheses, “Kruskal-Wallis” value p > 0.05 demonstrates that displacement of (d) ≈ (a) ≈ (b) ≈ (c) shows nothing different significantly by comparison with 600 N load. If stresses are tested upon prosthesis, (d) ≈ (a) ≈ (b) ≈ (c) is also displayed; if upon femora, (d) ≈ (a) ≈ (b) < (c) is suggested; if upon integral joint, (d) ≈ (a) < (b) < (c) is presented.

Conclusions

Mechanically, these four sorts of artificial joint replacement are stabilized in quantity. Cemented short-stem prostheses present the biggest stress, while porous short-stem & cemented long-stem designs are equivalently better than porous long-stem prostheses and alternatives for femoral-head replacement. The preferred design of those two depends on clinical conditions. The cemented long-stem is favorable for inactive elders with osteoporosis, and porously punched cementless short-stem design is suitable for patients with osteoporosis, while the porously punched cementless short-stem is favorable for those with a cement allergy. Clinically, the strength of this study is to enable preoperative strategy to provide acute correction and decrease procedure time.

Biomechanical advantages of robot-assisted pedicle screw fixation in posterior lumbar interbody fusion compared with freehand technique in a prospective randomized controlled trial-perspective for patient-specific finite element analysis

BACKGROUND CONTEXT

There have been conflicting results on the surgical outcome of lumbar fusion surgery using two different techniques: robot-assisted pedicle screw fixation and conventional freehand technique. In addition, there have been no studies about the biomechanical issues between both techniques.

PURPOSE

This study aimed to investigate the biomechanical properties in terms of stress at adjacent segments using robot-assisted pedicle screw insertion technique (robot-assisted, minimally invasive posterior lumbar interbody fusion, Rom-PLIF) and freehand technique (conventional, freehand, open approach, posterior lumbar interbody fusion, Cop-PLIF) for instrumented lumbar fusion surgery.

STUDY DESIGN

This is an additional post-hoc analysis for patient-specific finite element (FE) model.

PATIENT SAMPLE

The sample is composed of patients with degenerative lumbar disease.

OUTCOME MEASURES

Intradiscal pressure and facet contact force are the outcome measures.

METHODS

Patients were randomly assigned to undergo an instrumented PLIF procedure using a Rom-PLIF (37 patients) or a Cop-PLIF (41), respectively. Five patients in each group were selected using a simple random sampling method after operation, and 10 preoperative and postoperative lumbar spines were modeled from preoperative high-resolution computed tomography of 10 patients using the same method for a validated lumbar spine model. Under four pure moments of 7.5 Nm, the changes in intradiscal pressure and facet joint contact force at the proximal adjacent segment following fusion surgery were analyzed and compared with preoperative states.

RESULTS

The representativeness of random samples was verified. Both groups showed significant increases in postoperative intradiscal pressure at the proximal adjacent segment under four moments, compared with the preoperative state. The Cop-PLIF models demonstrated significantly higher percent increments of intradiscal pressure at proximal adjacent segments under extension, lateral bending, and torsion moments than the Rom-PLIF models (p=.032, p=.008, and p=.016, respectively). Furthermore, the percent increment of facet contact force was significantly higher in the Cop-PLIF models under extension and torsion moments than in the Rom-PLIF models (p=.016 under both extension and torsion moments).

CONCLUSIONS

The present study showed the clinical application of subject-specific FE analysis in the spine. Even though there was biomechanical superiority of the robot-assisted insertions in terms of alleviation of stress increments at adjacent segments after fusion, cautious interpretation is needed because of the small sample size.

Full-field measurement of micromotion around a cementless femoral stem using micro-CT imaging and radiopaque markers

A good primary stability of cementless femoral stems is essential for the long-term success of total hip arthroplasty. Experimental measurement of implant micromotion with linear variable differential transformers is commonly used to assess implant primary stability in pre-clinical testing. But these measurements are often limited to a few distinct points at the interface. New techniques based on micro-computed tomography (micro-CT) have recently been introduced, such as Digital Volume Correlation (DVC) or markers-based approaches. DVC is however limited to measurement around non-metallic implants due to metal-induced imaging artifacts, and markers-based techniques are confined to a small portion of the implant. In this paper, we present a technique based on micro-CT imaging and radiopaque markers to provide the first full-field micromotion measurement at the entire bone-implant interface of a cementless femoral stem implanted in a cadaveric femur. Micromotion was measured during compression and torsion. Over 300 simultaneous measurement points were obtained. Micromotion amplitude ranged from 0 to 24µm in compression and from 0 to 49µm in torsion. Peak micromotion was distal in compression and proximal in torsion. The technique bias was 5.1µm and its repeatability standard deviation was 4µm. The method was thus highly reliable and compared well with results obtained with linear variable differential transformers (LVDTs) reported in the literature. These results indicate that this micro-CT based technique is perfectly relevant to observe local variations in primary stability around metallic implants. Possible applications include pre-clinical testing of implants and validation of patient-specific models for pre-operative planning.

Integrative Performance Analysis of a Novel Bone Level Tapered Implant

Primary mechanical stability, as measured by maximum insertion torque and resonance frequency analysis, is generally considered to be positively associated with successful secondary stability and implant success. Primary implant stability can be affected by several factors, including the quality and quantity of available bone, the implant design, and the surgical procedure. The use of a tapered implant design, for instance, has been shown to result in good primary stability even in clinical scenarios where primary stability is otherwise difficult to achieve with traditional cylindrical implants—for example, in soft bone and for immediate placement in extraction sockets. In this study, bone-type specific drill procedures are presented for a novel Straumann bone level tapered implant that ensure maximum insertion torque values are kept within the range of 15 to 80 Ncm. The drill procedures are tested in vitro using polyurethane foam blocks of variable density, ex vivo on explanted porcine ribs (bone type 3), and finally in vivo on porcine mandibles (bone type 1). In each test site, adapted drill procedures are found to achieve a good primary stability. These results are further translated into a finite element analysis model capable of predicting primary stability of tapered implants. In conclusion, we have assessed the biomechanical behavior of a novel taper-walled implant in combination with a bone-type specific drill procedure in both synthetic and natural bone of various types, and we have developed an in silico model for predicting primary stability upon implantation.

Assessment of cartilage contact pressure and loading in the hip joint during split posture

PURPOSE

Given the crucial role of the mechanical behavior in the degenerative process of the hip joint, analyzing the contact mechanics in the articular layers during physical activities could contribute to understanding the pathology. Indeed, the development process of hip osteoarthritis generally evolves over a long time period, and therefore analyzing the mechanical behavior of the hip joint during extreme repetitive movements will be helpful to analyze degeneration causes. The aim of the study was to investigate the link between the excessive movements and the development of hip osteoarthritis.

METHODS

To individualize the analysis, we used a subject-specific and noninvasive approach based on finite element analysis and magnetic resonance imaging (MRI) data. The contact pressure distribution and loading conditions on the acetabular cartilage were assessed on eleven professional dancer subjects performing a split movement. This movement is frequently practiced (repetitive) by dancers during their daily exercises. Moreover, split postures are mostly characterized by high anatomical angles with subluxation (excessive). To ensure the motion accuracy, MRI data of the subjects were acquired in neutral and split positions performed inside the MRI scanner. Based on the reconstructed bone models from the MRI data, a motion tracking approach was used to compute the transformation between the two poses. To evaluate the contact during the split movement and to quantify the role of the labrum in the hip joint mechanics, additional simulations of two daily activities (walking and stand-up) were performed. Finally, a clinical study based on morphological and radiological analysis of the subjects was performed and validated by orthopedic surgeons and radiological experts to evaluate the proposed approach.

RESULTS

The reconstructed split movement was characterized by high anatomical angles with a subluxation on the left hip. Consequently, strong deformations and pressures were observed during the simulation. The comparison of the simulation results of split posture and daily activities showed higher pressure and lower contact area during extreme movements. Moreover, the presence of labrum absorbed part of load and consequently decreased the predicted contact pressure and contact area on the acetabular cartilage.

CONCLUSION

The comparison of the simulation results of the split posture and daily activities, as well as the correlation between the results of the analysis on extreme movement results and the clinical analysis performed by medical experts, strongly suggests that repetitive extreme movement could lead to early hip osteoarthritis.

Patient-specific bone modeling and analysis: the role of integration and automation in clinical adoption

Patient-specific analysis of bones is considered an important tool for diagnosis and treatment of skeletal diseases and for clinical research aimed at understanding the etiology of skeletal diseases and the effects of different types of treatment on their progress. In this article, we discuss how integration of several important components enables accurate and cost-effective patient-specific bone analysis, focusing primarily on patient-specific finite element (FE) modeling of bones. First, the different components are briefly reviewed. Then, two important aspects of patient-specific FE modeling, namely integration of modeling components and automation of modeling approaches, are discussed. We conclude with a section on validation of patient-specific modeling results, possible applications of patient-specific modeling procedures, current limitations of the modeling approaches, and possible areas for future research.

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.

Comparison of various functionally graded femoral prostheses by finite element analysis

This study is focused on finite element analysis of a model comprising femur into which a femoral component of a total hip replacement was implanted. The considered prosthesis is fabricated from a functionally graded material (FGM) comprising a layer of a titanium alloy bonded to a layer of hydroxyapatite. The elastic modulus of the FGM was adjusted in the radial, longitudinal, and longitudinal-radial directions by altering the volume fraction gradient exponent. Four cases were studied, involving two different methods of anchoring the prosthesis to the spongy bone and two cases of applied loading. The results revealed that the FG prostheses provoked more SED to the bone. The FG prostheses carried less stress, while more stress was induced to the bone and cement. Meanwhile, less shear interface stress was stimulated to the prosthesis-bone interface in the noncemented FG prostheses. The cement-bone interface carried more stress compared to the prosthesis-cement interface. Stair climbing induced more harmful effects to the implanted femur components compared to the normal walking by causing more stress. Therefore, stress shielding, developed stresses, and interface stresses in the THR components could be adjusted through the controlling stiffness of the FG prosthesis by managing volume fraction gradient exponent.

Inter-individual variability of bone density and morphology distribution in the proximal femur and T12 vertebra

Bone geometry, density and microstructure can vary widely between subjects. Knowledge about this variation in a population is of interest in particular for the design of orthopedic implants and interventions. The goal of this study is to investigate the local variability of bone density and microstructural parameters between subjects using a novel inter-subject image registration approach. Human proximal femora of 29 and T12 vertebrae of 20 individuals were scanned using a HR-pQCT and a micro-CT system, respectively. A pre-defined iso-anatomic mesh template was morphed to each micro-CT scan. For each element bone volume fraction and other morphological parameters (Tb.Th, Tb.N, Tb.Sp, SMI, DA) were determined and assigned to the element. A coefficient of variation (CV) was calculated for each parameter at each element location of the 29 femora and 20 T12 vertebrae. Contour plots of the CV distribution revealed very detailed information about the inter-individual variation in bone density and morphology. It is also shown that analyzing large sub-volumes, as commonly done in previous studies, would miss much of this variation. Detailed quantitative information of bone morphological parameters for each sample in the femur and the T12 database and their inter-individual variability are available from the mesh templates as supplementary data (http://w3.bmt.tue.nl/nl/fe_database/). We expect that these results can help to optimize implants and orthopedic procedures by taking local bone morphological parameter variations into account.

Estimation of change of bone structures after total hip replacement using bone remodeling simulation

BACKGROUND

The principal cause of femoral stem failure is the loosening of the total hip replacement due to bone resorption in the vicinity of the prosthesis (stress shielding). Bone rebuilds its structure continuously according to the daily mechanical stimuli. Therefore, surgical intervention alters the mechanical condition of bone severely. In this study, we propose a method to predict the change of bone structure after total hip replacement using bone remodeling simulation.

METHOD

The bone-stem complex structure model after total hip replacement was reconstructed based on CT-images used for preoperative planning by orthopedic experts. The bone remodeling simulation was conducted under the daily loading condition using our previous remodeling model, and the average equivalent stresses in the Gruen zone were evaluated.

FINDINGS

The predicted bone loss relevant to stress shielding was consistent to follow-up clinical data. Moreover, the remodeling simulation when using the stems of different size for the same patient could detect the size-dependent change of stress in the Gruen zone. In particular, the zone under the neck of the stem showed significant changes of stress and large bone loss, accompanying the risk of loosening or fracture.

INTERPRETATIONS

Prediction of bone structure changes after total hip replacement gives us significant information for longevity of prosthesis. Simulation results showed that the present computational framework could be considered to have potential in preoperative planning of total hip replacement.

Total hip arthroplasty by using a cementless ultrashort stem: A subject-specific finite element analysis for a young patient clinical case

In this article, a subject-specific finite element analysis has been developed to study a clinical case of a surgically misaligned hip prosthesis with an ultrashort stem. It was set out to study the strain energy density pattern, comparing the results obtained with computed tomography images. The authors developed two other numerical models: the first one analyzes the stress and strain distributions in the healthy femur (without prosthesis) and the second one analyzes the same boneimplant biomechanical system of the clinical case but assuming the prosthesis in the proper position. The misaligned prosthesis produced an overload at the proximal posterior plane of the femur, as confirmed by computed tomography images, which detect the formation of new bone. The numerical model of the correctly positioned prosthesis demonstrated that the bone is not overloaded and that the position of neutral axis does not significantly shift from the physiological condition.

SOA-based digital library services and composition in biomedical applications

Carefully collected, high-quality data are crucial in biomedical visualization, and it is important that the user community has ready access to both this data and the high-performance computing resources needed by the complex, computational algorithms that will process it. Biological researchers generally require data, tools and algorithms from multiple providers to achieve their goals. This paper illustrates our response to the problems that result from this. The Living Human Digital Library (LHDL) project presented in this paper has taken advantage of Web Services to build a biomedical digital library infrastructure that allows clinicians and researchers not only to preserve, trace and share data resources, but also to collaborate at the data-processing level.

Validation of FE micromotions and strains around a press-fit cup: introducing a new micromotion measuring technique

Finite element (FE) analysis provides an useful tool with which to analyze the potential performance of implantations in a variety of surgical, patient and design scenarios. To enable the use of FE analysis in the investigation of such implants, models must be experimentally validated. Validation of a pelvic model with an implanted press-fit cup in terms of micromotion and strain is presented here. A new method of micromotion has been introduced to better describe the overall movement of the cup within the pelvis. The method uses a digitizing arm to monitor the relative movement between markers on the cup and the surrounding acetabulum. FE analysis was used to replicate an experimental set up using a synthetic hemi-pelvis with a press-fitted all-metal cup, subject to the maximum loading observed during normal walking. The work presented here has confirmed the ability of FE models to accurately describe the mechanical performance of the press-fitted acetabulum and surrounding bone under typical loading conditions in terms of micromotion and strain distribution, but has demonstrated limitations in its ability to predict numerical micromotion values. A promising digitizing technique for measuring acetabular micromotions has also been introduced.

Application of finite element analysis in pre-operative planning for deformity correction of abnormal hip joints--a case series

In experimental and clinical research, it is difficult to directly measure responses in the human body, such as contact pressure and stress in a joint, but finite element analysis (FEA) enables the examination of in vivo responses by contact analysis. Hence, FEA is useful for pre-operative planning prior to orthopaedic surgeries, in order to gain insight into which surgical options will result in the best outcome. The present study develops a numerical simulation technique based on FEA to predict the surgical outcomes of osteotomy methods for the treatment of slipped capital femoral epiphyses. The correlation of biomechanical parameters including contact pressure and stress, for moderate and severe cases, is investigated. For severe slips, a base-of-neck osteotomy is thought to be the most reliable and effective surgical treatment, while any osteotomy may produce dramatic improvement for moderate slips. This technology of pre-operative planning using FEA can provide information regarding biomechanical parameters that might facilitate the selection of optimal osteotomy methods and corresponding surgical options.

Automated preoperative planning of femoral stem in total hip arthroplasty from 3D CT data: atlas-based approach and comparative study

Atlas-based methods for automated preoperative planning of the femoral stem implant in total hip arthroplasty are described. Statistical atlases are constructed from a number of past preoperative plans prepared by experienced surgeons in order to represent the surgeon's expertise of the planning. Two types of atlases are considered. One is a statistical distance map atlas, which represents surgeon's preference of the contact pattern between the femoral canal (host bone) and stem (implant) surfaces. The other is an optimal reference plan, which is selected as the best representative plan expected to minimize the deviation from the surgeon's preferred contact pattern. These atlases are fitted to the patient data to automatically generate the preoperative plan of the femoral stem. In this paper, we formulate a general framework of atlas-based implant planning, and then describe the methods for construction and utilization of the two proposed atlases. In the experiments, we used 40 cases to evaluate the proposed methods and compare them with previous methods by defining the errors as differences between automated and surgeon's plans. By using the proposed methods, the positional and orientation errors were significantly reduced compared with the previous methods and the size error was superior to inter-surgeon variability in size selection using 2D templates on an X-ray image reported in previous work.

Torsional stability of uncemented femoral stems in oncologic reconstructions

Modular oncology implants using uncemented fixation represent a popular reconstruction technique for limb salvage patients. Initial stability is critical to facilitate bony ingrowth of host bone into the stem of a press-fit oncologic modular rotating-hinge total knee arthroplasty (TKA). The impact of stem design on initial stability has not been defined. The goal of this study was to evaluate the initial stability of 3 different stem designs as defined by torsional load to failure. An analysis of imaging was also performed. The pilot study consisted of 5 femora in each of 3 groups based on stem design. The specimen was mounted on a multi-axis biomechanical test frame equipped with a Vicon 3D motion analysis 4-camera system (Vicon Motion Systems, Lake Forest, California) to track the relative motion between the implant and the femur. Torsional force was applied until failure. The straight-fluted stem design had the highest average torsional stiffness (18.3±8.2 Nm/deg) and average torque at 150 μm of implant micromotion (23.2±10.6 Nm) of the 3 stem types tested.The results of this study will help to guide surgical decision making in limb salvage cases. Further investigation is warranted.

Experimental versus Computational Analysis of Micromotions at the Implant—Bone Interface

In total hip arthroplasty, micromotions at the implant—bone interface influence the long-term survival of the prosthesis. These micromotions are often measured using sensors that are fixed to the implant and bone at points that are remote from the interface. Given that the implant—bone system is not rigid, errors may be introduced. It is not possible to assess the magnitude of these errors with the currently available experimental methods. However, this problem can be investigated using the finite element method (FEM).

The hypothesis that the actual interface micromotions differ from those measured in the experimental manner was tested using a case-specific FE model, validated against deflection experiments. The FE model was used to simulate an ‘experimental’ method to measure micromotions. This ‘experimental’ method was performed by mimicking the distance between the measurement points; the implant point was selected at the interface while the bony point was at the outer surface of bone.

No correlation was found between the micromotions computed at the interface and when using remote reference points. Moreover, the magnitudes of micromotions computed with the latter method were considerably greater. By reducing the distance between the reference points the error decreased, but the correlation stayed unchanged. Care needs to be taken when interpreting the results of micromotion measurement systems that use bony reference points at a distance from the actual interface.

Finite Element Modeling of Resurfacing Hip Prosthesis: Estimation of Accuracy Through Experimental Validation

Metal-on-metal hip resurfacing is becoming increasingly popular, and a number of new devices have been recently introduced that, in the short term, appear to have satisfactory outcome but many questions are still open on the biomechanics of the resurfaced femur. This could be investigated by means of finite element analysis, but, in order to be effective in discerning potential critical conditions, the accuracy of the models’ predictions should be assessed. The major goal of this study was to validate, through a combined experimental-numerical study, a finite element modeling procedure for the simulation of resurfaced femurs. In addition, a preliminary biomechanical analysis of the changes induced in the femoral neck biomechanics by the presence of the device was performed, under a physiologic range of hip joint reaction directions. For this purpose, in vitro tests and a finite element model based on the same specimen were developed using a cadaver femur. The study focused on the Conserve Plus, one of the most common contemporary resurfacing designs. Five loading configurations were identified to correspond to the extremes of physiological directions for the hip joint. The agreement between experimental measurements and numerical predictions was good both in the prediction of the femoral strains (R2>0.9), and in the prosthesis micromotions (error<20 μm), giving confidence in the model predictions. The preliminary biomechanical analysis indicated that the strains in the femoral neck are moderately affected by the presence of the prosthesis, apart from localized strain increments that can be considerable, always predicted near the stem. Low micromotions and contact pressure were predicted, suggesting a good stability of the prosthesis. The model accuracy was good in the prediction of the femoral strains and moderately good in the prediction of the bone-prosthesis micromovements. Although the investigated loading conditions were not completely physiological, the preliminary biomechanical analysis showed relatively small changes for the proximal femur after implantation. This validated model can support realistic simulations to examine physiological load configurations and the effects of variations in prosthesis design and implantation technique.

Finite element simulation of cement-bone interface micromechanics: a comparison to experimental results

Recently, experiments were performed to determine the micromechanical behavior of the cement-bone interface under tension-compression loading conditions. These experiments were simulated using finite element analysis (FEA) to test whether the micromechanical response of the interface could be captured in micromodels. Models were created of experimental specimens based upon microcomputed tomography data, including the complex interdigitated bone-cement morphology and simulated frictional contact at the interface. The models were subjected to a fully reversed tension-compression load, mimicking the experimental protocol. Similar to what was found experimentally, the simulated interface was stiffer in compression than in tension, and the majority of displacement was localized to the cement-bone interface. A weak correlation was found between the FEA-predicted stiffness and the stiffness found experimentally, with average errors of 8 and 30% in tension and compression, respectively. The hysteresis behavior found experimentally was partially reproduced in the simulation by including friction at the cement-bone interface. Furthermore, stress analysis suggested that cement was more at risk of fatigue failure than bone, concurring with the experimental observation that more cracks were formed in the cement than in the bone. The current study provides information that may help explain the load transfer mechanisms taking place at the cement-bone interface.

Subject specific finite element analysis of implant stability for a cementless femoral stem

BACKGROUND

The primary stability of a cementless implant is crucial to ensure long term stability through osseointegration. In the present study we have examined how subject specific finite element models can be used to evaluate the stability of a cementless femoral stem.

METHODS

Micromotion on the bone-implant interface of a cementless stem was measured experimentally in six human cadaver femurs. Subject specific finite element models were built from computed tomography of the same femurs, and used to simulate the same load scenario used experimentally.

FINDINGS

Both experimental measurements and numerical analyses showed a tendency of increased rotational stability for bigger implants. Good correlation was found between measurements and calculated values of axial rotation (R(2)=0.74, P<0.001). The finite element models produced interface micromotion of the same magnitude as measured experimentally, with micromotion generally below 40 microm. Bigger femoral stems were found to decrease the micromotion in the experimental measurements. This tendency could not be recognised in the interface micromotion from the finite element models.

INTERPRETATION

The finite element models showed limited success in predicting interfacial micromotion, but reproduced a similar pattern of rotational stability for the implants as seen experimentally. Since rotation in retroversion is often the main concern when studying implant stability, subject specific finite element models could be employed for pre-clinical evaluation of implants.

Multiscale modelling of the skeleton for the prediction of the risk of fracture

BACKGROUND

The development of a multiscale model of the human musculoskeletal system able to accurately predict the risk of bone fracture is still a grand challenge. The aim of this paper is to present the Living Human Project, to describe the final system and to review the achievements obtained so far. The Living Human musculoskeletal supermodel is conceived as the interconnection of five interdependent sub-models: the continuum, the boundary condition, the constitutive equation, the remodelling history and the failure criterion sub-models.

METHODS

Methods are available to develop accurate subject-specific finite element models of bones that can incorporate the subject's tissue-density distribution and empirically derived constitutive laws. Anatomo-functional musculoskeletal models can be registered with gait analysis data to predict muscle and joint forces acting on the patient's skeleton during gait. These are the boundary conditions for the continuum models that showed an average error of 12% in the prediction of the failure load. Still, the entire supermodel is defined as a collection of procedural macros to predict the risk of fracture and should be improved.

FINDINGS

Even with these limitations, the organ-level model already found some clinically relevant applications, especially in the analysis of joint prostheses. Also, the body-organ level multiscale model finds some clinical applications in paediatric skeletal oncology. The tissue- and the cell-level models are not yet fully validated. Thus, they cannot be safely used in clinical applications.

INTERPRETATION

The continuum sub-model is the most mature model available. More powerful methods are needed for the generation of anatomo-functional musculoskeletal models. Muscle force prediction should be improved, investigating new probabilistic approaches to identify the neuro-motor strategy. The changes of the tissue properties in the various regions of the skeleton and predictive remodelling models should be included. An adequate information technology infrastructure should be developed to support collaborative work and integration of different sub-models.