Effect of glenoid prosthesis design on glenoid bone remodeling: Adaptive finite element based simulation

Glenoid prosthesis loosening is the most common cause for revision total shoulder arthroplasty. Improved glenoid prosthesis design requires looking beyond initial post-implantation static stress analyses. Adaptive bone remodeling simulations based on Wolff’s law are needed for predicting long-term glenoid prosthesis results. This study demonstrates the capability of predicting glenoid bone remodeling produced by changing prosthesis design features. Twelve glenoid prostheses were designed to fit each of six donor human glenoids, using combinations of three peg types and four backing-peg material combinations (polyethylene and or metal). The twelve FE prosthesis models were individually combined, simulating surgical implantation, with the glenoid models. Remodeling simulations, using a validated adaptive bone remodeling simulation, commenced with homogeneous glenoid bone density. To produce bone remodeling, center, posterior-offset, and anterior-offset physiologic loads were consecutively applied to the bone–prosthesis FE models for 300 iterations. Upon completion, region-specific mean predicted bone apparent densities were compared between bone–prosthesis and intact glenoid FE models. Metal fixations significantly increased proximal-center bone density. Polyethylene fixations resulted in bone density approximately equal to intact. Two angled polyethylene peg designs with longer-anterior and shorter-posterior pegs, reflecting natural glenoid shape, best maintained mid and distal glenoid bone density. While these initial results were not validated, they demonstrate the capability and potential of adaptive glenoid bone remodeling simulation. We expect FE glenoid bone remodeling simulations to become powerful and robust tools in the design and evaluation of glenoid prostheses.     

Topological optimization in hip prosthesis design

With particular interest on total hip arthroplasty (THA), optimization of orthopedic prostheses is employed in this work to minimize the probability of implant failure or maximize prosthesis reliability. This goal is often identified with the reduction of stress concentrations at the interface between bone and these devices. However, aseptic loosening of the implant is mainly influenced by bone resorption phenomena revealed in some regions of the femur when a prosthesis is introduced. As a consequence, bone resorption appears due to stress shielding, that is to say the decrease of the stress level in the implanted femur caused by the significant load carrying of the prosthesis due to its higher stiffness. A maximum stiffness topological optimization-based (TO) strategy is utilized for non-linear static finite element (FE) analyses of the femur–implant assembly, with the goal of reducing stress shielding in the femur and to furnish guidelines for re-designing hip prostheses. This is accomplished by employing an extreme accuracy for both the three-dimensional reconstruction of the femur geometry and the material properties maps assigned as explicit functions of the local densities.     

Three-dimensional finite element analysis of glenoid replacement prostheses: a comparison of keeled and pegged anchorage systems

Glenoid component loosening is the dominant cause of failure in total shoulder arthroplasty. It is presumed that loosening in the glenoid is caused by high stresses in the cement layer. Several anchorage systems have been designed with the aim of reducing the loosening rate, the two major categories being "keeled" fixation and "pegged" fixation. However, no three-dimensional finite element analysis has been performed to quantify the stresses in the cement or to compare the different glenoid prosthesis anchorage systems. The objective of this study was to determine the stresses in the cement layer and surrounding bone for glenoid replacement components. A three-dimensional model of the scapula was generated using CT data for geometry and material property definition. Keeled and pegged designs were inserted into the glenoid, surrounded by a 1-mm layer of bone cement. A 90 deg arm abduction load with a full muscle and joint load was applied, following van der Helm (1994). Deformations of the prosthesis, stresses in the cement, and stresses in the bone were calculated. Stresses were also calculated for a simulated case of rheumatoid arthritis (RA) in which bone properties were modified to reflect that condition. A maximum principal stress-based failure model was used to predict what quantity of the cement is at risk of failure at the levels of stress computed. The prediction is that 94 percent (pegged prosthesis) and 68 percent (keeled prosthesis) of the cement has a greater than 95 percent probability of survival in normal bone. In RA bone, however, the situation is reversed where 86 percent (pegged prosthesis) and 99 percent (keeled prosthesis) of the cement has a greater than 95 percent probability of survival. Bone stresses are shown to be not much affected by the prosthesis design, except at the tip of the central peg or keel. It is concluded that a "pegged" anchorage system is superior for normal bone, whereas a "keeled" anchorage system is superior for RA bone.     

Bone ingrowth simulation for a concept glenoid component design

Glenoid component loosening is the major problem of total shoulder arthroplasty. It is possible that uncemented component may be able to achieve superior fixation relative to cemented component. One option for uncemented glenoid is to use porous tantalum backing. Bone ingrowth into the porous backing requires a degree of stability to be achieved directly post-operatively. This paper investigates the feasibility of bone ingrowth with respect to the influence of primary fixation, elastic properties of the backing and friction at the bone prosthesis interface. Finite element models of three glenoid components with different primary fixation configurations are created. Bone ingrowth into the porous backing is modelled based on the magnitude of the relative interface micromotions and mechanoregulation of the mesenchymal stem cells that migrated via the bonded part of the interface. Primary fixation had the most influence on bone ingrowth. The simulation showed that its major role was not to firmly interlock the prosthesis, but rather provide such a distribution of load, that would result in reduction of the peak interface micromotions. Should primary fixation be provided, friction has a secondary importance with respect to bone ingrowth while the influence of stiffness was counter intuitive: a less stiff backing material inhibits bone ingrowth by higher interface micromotions and stimulation of fibrous tissue formation within the backing.     

Finite element modelling of glenohumeral kinematics following total shoulder arthroplasty

Due to the shallowness of the glenohumeral joint, a challenging but essential requirement of a glenohumeral prosthesis is the prevention of joint dislocation. Weak glenoid bone stock and frequent dysfunction of the rotator cuff, both of which are common with rheumatoid arthritis, make it particularly difficult to achieve this design goal. Although a variety of prosthetic designs are commercially available only a few experimental studies have investigated the kinematics and dislocation characteristics of design variations. Analytical or numerical methods, which are predictive and more cost-effective, are, apart from simple rigid-body analyses, non-existent. The current investigation presents the results of a finite element analysis of the kinematics of a total shoulder joint validated using recently published experimental data for the same prostheses. The finite element model determined the loading required to dislocate the humeral head, and the corresponding translations, to within 4% of the experimental data. The finite element method compared dramatically better to the experimental data (mean difference=2.9%) than did rigid-body predictions (mean difference=37%). The goal of this study was to develop an accurate method that in future studies can be used for further investigations of the effect of design parameters on dislocation, particularly in the case of a dysfunctional rotator cuff. Inherently, the method also evaluates the glenoid fixation stresses in the relatively weak glenoid bone stock. Hence, design characteristics can be simultaneously optimised against dislocation as well as glenoid loosening.     

Quantifying the competing relationship between adduction range of motion and baseplate micromotion with lateralization of reverse total shoulder arthroplasty

Lateralizing the center of rotation (COR) of reverse total shoulder arthroplasty (rTSA) could improve functional outcomes and mitigate scapular notching, a commonly occurring complication of the procedure. However, resulting increases in torque at the bone-implant interface may negatively affect initial fixation of the glenoid-side component, especially if only two fixation screws can be placed. Shoulder-specific finite element (FE) models of four fresh-frozen cadaveric shoulders were constructed. Scapular geometry and material property distributions were derived from CT data. Generic baseplates with two and four fixation screws were virtually implanted, after which superiorly-oriented shear loads, accompanied by a compressive load, were applied incrementally further from the glenoid surface to simulate lateralization of the COR. Relationships between lateralization, adduction range of motion (ROM), the number of fixation screws and micromotion of the baseplate (initial implant fixation) were characterized. Lateralization significantly increases micromotion (p=0.015) and adduction ROM (p=0.001). Using two, versus four, baseplate fixation screws significantly increases micromotion (p=0.008). The effect of lateralization and the number of screws on adduction ROM and baseplate fixation is variable on a shoulder-specific basis. Trade-offs exist between functional outcomes, namely adduction ROM, and initial implant fixation and the negative effect of lateralization on implant fixation is amplified when only two fixation screws are used. The possibility of lateralizing the COR in order to improve functional outcomes of the procedure should be considered on a patient-specific basis accounting for factors such as availability and quality of bone stock.     

Analysis of bone–prosthesis interface micromotion for cementless tibial prosthesis fixation and the influence of loading conditions

A lack of initial stability of the fixation is associated with aseptic loosening of the tibial components of cementless knee prostheses. With sufficient stability after surgery, minimal relative motion between the prosthesis and bone interfaces allows osseointegation to occur thereby providing a strong prosthesis-to-bone biological attachment. Finite element modelling was used to investigate the bone–prosthesis interface micromotion and the relative risk of aseptic loosening. It was anticipated that by prescribing different joint loads representing gait and other activities, and the consideration of varying tibial–femoral contact points during knee flexion, it would influence the computational prediction of the interface micromotion. In this study, three-dimensional finite element models were set up with applied loads representing walking and stair climbing, and the relative micromotions were predicted. These results were correlated to in-vitro measurements and to the results of prior retrieval studies. Two load conditions, (i) a generic vertical joint load of 3×body weight with 70%/30% M/L load share and antero-posterior/medial-lateral shear forces, acted at the centres of the medial and lateral compartments of the tibial tray, and (ii) a peak vertical joint load at 25% of the stair climbing cycle with corresponding antero-posterior shear force applied at the tibial–femoral contact points of the specific knee flexion angle, were found to generate interface micromotion responses which corresponded to in-vivo observations. The study also found that different loads altered the interface micromotion predicted, so caution is needed when comparing the fixation performance of various reported cementless tibial prosthetic designs if each design was evaluated with a different loading condition.     

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

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

Multi-factorial analysis of variables influencing the bone loss of an implant placed in the maxilla: Prediction using FEA and SED bone remodeling algorithm

The aim of this study was to investigate the interactions of implant position, implant–abutment connection and loading condition influencing bone loss of an implant placed in the maxilla using finite element (FE) analysis and mathematical bone remodeling theory. The maxilla section contours were acquired using CT images to construct FE models containing RS (internal retaining-screw) and the TIS (taper integrated screwed-in) implants placed in SC (along the axis of occlusal force) and RA (along the axis of residual ridge) positions. The adaptive strain energy density (SED) algorithm was combined with FE approach to study the preliminary bone remodeling around implant systems under different load conditions. The simulated results showed that the implant position obviously influenced the bone loss. An implant placed in the RA position resulted in substantially increased bone loss. Implant receiving a lateral load slightly increased bone loss compared with an axial load. The implant type did not significantly influence bone loss. It was found that buccal site suffered the most bone loss around the implant, followed by distal, lingual and mesial sites. The implant position primarily influenced bone loss and it was found most obviously at the buccal site. Implant placed along the axial load direction of a proposed prosthesis could obtain less bone loss around the implant. Attaining proper occlusal adjustments to reduce the lateral occlusal force is recommended in implant–bone–prosthesis system. Abutments of internal engagement with or without taper-fit did not affect the bone loss in the surrounding bone.     

Bone remodelling of the scapula after a total shoulder arthroplasty

According to Wolff’s law, the changes in stress after a prosthesis implantation may modify the shape and internal structure of bone, thus compromising the long-term prosthesis fixation and, consequently, be a significant factor for glenoid loosening. The aim of the present study is to evaluate the changes in the bone adaptation process of the scapula after an anatomical and reverse total shoulder arthroplasty. Five finite element models of the implanted scapula are developed considering the implantation of three anatomical, cemented, all-polyethylene components; an anatomical, cementless, metal-backed component; and a reverse, all-metal component. The methodology followed to simulate the bone adaptation of the scapula was previously validated for the intact model, prior to the prosthesis implantation. Additionally, the influence of the bone quality on the adaptation process is also investigated by considering an osteoporotic condition. The results show that the stress shielding phenomenon is more concerning in cementless, metal-based components than in cemented, all-polyethylene components, regardless of the bone quality. Consequently, as far as the bone adaptation process of the bone is concerned, cemented, all-polyethylene components are better suited for the treatment of the shoulder joint.     

Acromion-fixation of glenoid components in total shoulder arthroplasty

Successful design of components for total shoulder arthroplasty has proven to be challenging. This is because of the difficulties in maintaining fixation of the component that inserts into the scapula; i.e., the glenoid component. Glenoid components that are fixated to both the glenoid and acromion (a long process extending medially on the dorsal aspect of the scapula) have the possible advantage of greater stability over those that are fixated to the glenoid alone. In this study, a finite element analysis is used to investigate whether or not acromion fixation is advantageous for glenoid components. Full muscle loading and joint reaction forces are included in the finite element model. Reflective photoelasticity of five scapulae is used to obtain experimental data to compare with results from the finite element analysis, and it confirms the structural behaviour of the finite element model. When implanted with an acromion-fixated prosthesis, it is found that high unphysiological stresses occur in the scapula bone, and that stresses in the fixation are not reduced. Very high stresses are predicted in that part of the prosthesis which connects the acromion to the glenoid. It is found that the very high stresses are partly in response to the muscle and joint reaction forces acting at the acromion. It is concluded that, because of the relatively high forces acting at the acromion, fixation to it may not be the way forward in glenoid component design.     

Subject-specific finite element model of the pelvis: development, validation and sensitivity studies

A better understanding of the three-dimensional mechanics of the pelvis, at the patient-specific level, may lead to improved treatment modalities. Although finite element (FE) models of the pelvis have been developed, validation by direct comparison with subject-specific strains has not been performed, and previous models used simplifying assumptions regarding geometry and material properties. The objectives of this study were to develop and validate a realistic FE model of the pelvis using subject-specific estimates of bone geometry, location-dependent cortical thickness and trabecular bone elastic modulus, and to assess the sensitivity of FE strain predictions to assumptions regarding cortical bone thickness as well as bone and cartilage material properties. A FE model of a cadaveric pelvis was created using subject-specific computed tomography image data. Acetabular loading was applied to the same pelvis using a prosthetic femoral stem in a fashion that could be easily duplicated in the computational model. Cortical bone strains were monitored with rosette strain gauges in ten locations on the left hemipelvis. FE strain predictions were compared directly with experimental results for validation. Overall, baseline FE predictions were strongly correlated with experimental results (r2=0.824), with a best-fit line that was not statistically different than the line y=x (experimental strains = FE predicted strains). Changes to cortical bone thickness and elastic modulus had the largest effect on cortical bone strains. The FE model was less sensitive to changes in all other parameters. The methods developed and validated in this study will be useful for creating and analyzing patient-specific FE models to better understand the biomechanics of the pelvis.     

Predicting surface strains at the human distal radius during an in vivo loading task — Finite element model validation and application

Bone strains resulting from physical activity are thought to be a primary driver of bone adaptation, but cannot be directly noninvasively measured. Because bone adapts nonuniformly, physical activity may make an important independent structural contribution to bone strength that is independent of bone mass and density. Our objective was to create and validate methods for subject-specific finite element (FE) model generation that would accurately predict the surface strains experienced by the distal radius during an in vivo loading task, and to apply these methods to a group of 23 women aged 23–35 to examine variations in strain, bone mass and density, and physical activity. Four cadaveric specimens were experimentally tested and specimen-specific FE models were developed to accurately predict periosteal surface strains (root mean square error=16.3%). In the living subjects, when 300 N load was simulated, mean strains were significantly inversely correlated with BMC (r=−0.893), BMD (r=−0.892) and physical activity level (r=−0.470). Although the group of subjects was relatively homogenous, BMD varied by two-fold (range: 0.19–0.40 g/cm³) and mean energy-equivalent strain varied by almost six-fold (range: 226.79–1328.41 με) with a simulated 300 N load. In summary, we have validated methods for estimating surface strains in the distal radius that occur while leaning onto the palm of the hand. In our subjects, strain varied widely across individuals, and was inversely related to bone parameters that can be measured using clinical CT, and inversely related to physical activity history.     

The contribution of the glenoid labrum to glenohumeral stability under physiological joint loading using finite element analysis

The labrum contributes to passive glenohumeral joint stability. Cadaveric studies have demonstrated that this has position and load dependency, which has not been quantified under physiological loads. This study aims to validate subject-specific finite element (FE) models against in vitro measurements of joint stability and to utilise the FE models to predict joint stability under physiological loads. The predicted stability values were within ± one standard deviation of experimental data and the FE models showed a reduction in stability of 10–15% with high, physiological, loads. The developed regression equations provide the first representation of passive glenohumeral stability and will aid surgical decision-making.     

Glenoid cancellous bone strength and modulus.

The objectives of this study were to determine the strength and modulus of glenoid cancellous bone, including regional variations. The motivations were: to select a suitable bone substitute for standardized testing of glenoid prosthesis loosening, to assist in shoulder prosthesis design and to provide input data for finite element analyses. Ten glenoids from eight cadavers (mean age, 81) were tested by in situ indentation. Mean strength ranged from 6.7 to 17 MPa for the ten glenoids, the overall mean being 10.3 MPa. Mean E moduli ranged from 67 to 171 MPa for the individual glenoids, the overall mean being 99 MPa. These values are likely at the lower end of what would be expected for normal bone since strength and modulus decrease with age and the available specimens were older. These values may be appropriate for prosthesis design, however, since mechanical properties are reduced in rheumatoid arthritic bone. Regional trends were very similar for modulus and strength. The strongest region was postero-superior. The central column, correlating with the keel position in many glenoid components, was weaker than both the anterior and posterior regions but deeper. A large drop in strength and modulus below the subchondral layer emphasizes the importance of maintaining this layer during prosthetic replacement.     

Experimental method for the in vitro testing of the initial stability of cementless hip prostheses

Micromotions at the interface between bone and prosthesis are believed to induce bone resorption and ultimately lead to loosening of the implant. Thus the initial stability achieved by a hip prosthesis is an important factor for the long-term function of the implant. Knowing the biological consequences of the mechanical conditions, it appears to be mandatory to measure the extent of these three-dimensional movements. An in vitro dynamic method for measurement of the micromotion of the femoral component of hip prostheses has been developed. Tests in cemented prostheses have confirmed that the use of cement reduces sinkage and rotation manyfold and have yielded reference values for stability. Comparison with two types of cementless prostheses has shown that certain cementless implants may achieve stability comparable to cemented ones in some load directions.     

Analysis of glenoid fixation for a reversed anatomy fixed-fulcrum shoulder replacement

The restoration of pain-free stable function in gleno-humeral arthritic cases in various situations such as rotator cuff deficiency, old trauma and failed total shoulder arthroplasty is a challenging clinical dilemma. The Bayley-Walker shoulder has been designed specifically for very difficult cases where surface replacement devices do not provide sufficient stability. This device is a fixed-fulcrum reverse anatomy implant where the centre of rotation is placed medially and distally with respect to the normal shoulder, to increase the lever arm of the abductor muscles. An important problem in devices of this type is obtaining secure and long-lasting fixation of the glenoid component. In this design, fixation is achieved using a tapered screw for engagement with cortical bone and HA coating for subsequent osseointegration. This study presents the results from a three-dimensional finite element analysis conducted on this component for two load cases at 60 degrees and 90 degrees abduction. The results showed that most of the forces were transmitted from the component to the cortical bone of the scapula, the remaining load being transmitted through cancellous bone. Histology from a retrieved case shows evidence of bone remodelling. The retrieval case obtained some time after implantation showed new bone formation had occurred around the threads onto the HA coating. Fixing the component in this way at multiple locations in cortical bone may overcome the problems of glenoid loosening historically associated with cemented constrained devices.     

Influence of transmucosal height in abutments of single and multiple implant-supported prostheses: a non-linear three-dimensional finite element analysis

The aim of this study was to analyze the influence of three different transmucosal heights of the abutments in single and multiple implant-supported prostheses through the finite element method. External hexagon implants, MicroUnit, and EsthetiCone abutments were scanned and placed in an edentulous maxillary model obtained from a tomography database. The simulations were divided into two groups: (1) one implant with 3.75 × 10 mm placed in the upper central incisor, simulating a single implant-supported fixed prosthesis with an EsthetiCone abutment; and (2) two implants with 3.75 × 10 mm placed in the upper lateral incisors with MicroUnit abutments, simulating a multiple implant-supported prosthesis. Subsequently, each group was subdivided into three models according to the transmucosal height (1, 2, and 3 mm). A static oblique load at an angle of 45 degrees to the long axis of the implant in palatal-buccal direction of 150 and 75 N was applied for multiple and single implant-supported prosthesis, respectively. The implants and abutments were assessed according to the equivalent Von Mises stress analyses while the bone and ceramics were analyzed through maximum and minimum principal stresses. The total deformation values increased in all models, while the transmucosal height was augmented. The transmucosal height of the abutments influences the stress values at the bone, ceramics, implants, and abutments of both the single and multiple implant-supported prostheses, with the transmucosal height of 1 mm showing the lowest stress values.     

Design of partial pelvic replacement with calculation of stresses on the pelvic bone using the finite element method

Although tumors of the pelvic region are relatively rare, with regard to the provision of an individual prosthesis, they make great demands both on the engineer and on the surgeon. In the case of partial pelvic replacement, the main problem is that of fixation. The direction of the preload by the screw and counternut in the anchorage should make possible the introduction of forces into the bone that closely mimic the natural pre-operative situation, in order to keep bone remodeling and stressing to a minimum. A three-dimensional finite element (FE) analysis of the stresses at work in the pelvic bone was performed using an FE model of the pelvis constructed on the basis of computed tomographic data, since with the FE method it is possible to obtain information about deformations, internal stresses and local forces acting on the bone. The results show that the main stressing of the bony pelvis occurs in the region of the acetabulum and the iliosacral joints, and that further points of loading are located in the lower region of the ischium-and this while standing on one leg, on the side of the supporting leg.     

Finite element assessment of block-augmented total knee arthroplasty

Loosening and migration of tibial prostheses have been identified as causes of early total knee replacement (TKR) failure. The problem is made more complex when defects occur in the proximal tibia compromising fixation and alignment. Clinical studies using metal augments have shown these to be an alternative to other means of defect treatment. Finite element (FE) analysis can be used to identify regions that may be prone to loosening and migration. In the current work, 3D FE models of TKR uncontained type-2 defects treated with block augments have been constructed and analysed. It has been shown that a metal augment is the most suitable. The use of bone cement (PMMA) to fill proximal defects is not considered suitable as stresses carried by the cement block exceed those of the fatigue limit of bone cement. It has been shown that the stresses in the proximal cancellous bone of block-augmented models are significantly below levels likely to cause damage due to overloading. Furthermore, the use of stem extensions has been shown to reduce the cancellous bone stresses in the proximal region thus increasing the likelihood of bone resorption. Given this, it is recommended that stem extensions are not required unless necessary to mitigate some other problem.