Revision surgeries for orthopaedic infections are done in two stages – one surgery to implant an antibiotic spacer to clear the infection and another to install a permanent implant. A permanent porous implant, that can be loaded with antibiotics and allow for single-stage revision surgery, will benefit patients and save healthcare resources. Gyroid structures can be constructed with high porosity, without stress concentrations that can develop in other period porous structures [1] [2]. The purpose of this research is to compare the resulting bone and prosthesis stress distributions when porous versus solid stems are implanted into three proximal humeri with varying bone densities, using finite element models (FEM).

Porous humeral stems were constructed in a gyroid structure at porosities of 60%, 70%, and 80% using computer-aided design (CAD) software. These CAD models were analyzed using FEM (Abaqus) to look at the stress distributions within the proximal humerus and the stem components with loads and boundary conditions representing the arm actively maintained at 120˚ of flexion. The stem was assumed to be made of titanium (Ti6Al4V). Three different bone densities were investigated, representing a healthy, an osteopenic, and an osteoporotic humerus, with an average bone shape created using a statistical shape and density model (SSDM) based on 75 cadaveric shoulders (57 males and 18 females, 73 12 years) [3]. The Young's moduli (E) of the cortical and trabecular bones were defined on an element-by-element basis, with a minimum allowable E of 15 MPa. The Von Mises stress distributions in the bone and the stems were compared between different stem scenarios for each bone density model.

A preliminary analysis shows an increase in stress values at the proximal-lateral region of the humerus when using the porous stems compared to the solid stem, which becomes more prominent as bone density decreases. With the exception of a few mesh dependent singularities, all three porous stems show stress distributions below the fatigue strength of Ti-6Al-4V (410 MPa) for this loading scenario when employed in the osteopenic and osteoporotic humeri [4]. The 80% porosity stem had a single strut exceeding the fatigue strength when employed in the healthy bone.

The results of this study indicate that the more compliant nature of the porous stem geometries may allow for better load transmission through the proximal humeral bone, better matching the stress distributions of the intact bone and possibly mitigating stress-shielding effects. Importantly, this study also indicates that these porous stems have adequate strength for long-term use, as none were predicted to have catastrophic failure under the physiologically-relevant loads. Although these results are limited to a single boney geometry, it is based on the average shape of 75 shoulders and different bone densities are considered. Future work could leverage the shape model for probabilistic models that could explore the effect of stem porosity across a broader population. The development of these models are instrumental in determining if these structures are a viable solution to combatting orthopaedic implant infections.

Statistical shape modeling (SSM) and statistical density modeling (SDM) are tools capable of describing the main modes of deviation in the shape and density distribution of the shoulder using a set of uncorrelated variables called principal components (PCs). We hypothesize that the first PC of the SDM, which scales overall density up/down, will be inversely correlated with age and will, on average, be greater for males than females. We also hypothesize that there is a correlation between some PCs of shape and density.

SSM and SDM were developed for scapulae and humeri by segmenting surface meshes from computed tomographic images of 75 cadaveric shoulders. Bones were co-registered and defined by the same surface mesh. Volumetric tetrahedral meshes were defined for one of the specimens serving as base meshes for SDM. Base meshes were morphed to each individual bone's surface and superimposed upon the corresponding CT data to determine image intensity in Hounsfield units at each node. Principal component analysis was performed on the exterior shape and internal density distribution of bones. T-tests were performed to find any differences in PC scores between males and females, and Pearson correlation coefficients were calculated for age and PC scores. Finally, correlation coefficients between each of the PCs of the shape and density models were calculated.

For the humerus, the first three PCs of the SDM were significantly correlated with age (ρ = 0.40, −0.46, and 0.36, all p ≤ 0.007). For the scapula, the first and ninth PCs showed such correlation (ρ = −0.31, and −0.32, all p ≤ 0.02). Statistically significant differences due to sex were found for the second to sixth SDM PCs of the humerus, with differences in average PC scores of 1, 1, −0.7, −0.8, and −0.6 standard deviations, respectively, for males relative to females. For the scapula, the second, fifth and seventh SDM PCs were significantly different between males and females, with average PC scores differing by 1.1, 0.7, and −0.6 standard deviations. Finally, for both bones, the first PC of SSM showed a weak but significant correlation with the second PC of the SDM (ρ = 0.47, p < 0.001 for the humerus, and ρ = 0.39, p < 0.001 for the scapula).

The results of this study suggest that age has a significant influence on the first PC of the SDM, associated with scaling the density in the cortical boundary. Moreover, the negative correlation of age with the second PC of the humerus in SDM which mostly influences the thickness of the cortical boundary implies cortical thinning with age. The second PC of both bones differed significantly between males and females, implying that cortical thickness differs between sexes. Also, there was a significant correlation between the size of the bones and the thickness of the cortical boundary. These findings can help guide the designs of population-based prosthesis components.

During revision total knee arthroplasty (rTKA), proximal tibial bone loss is frequently encountered and can result in a less-stable bone-implant fixation. A 3D printed titanium alloy (Ti6Al4V) revision augment that conforms to the irregular shape of the proximal tibia was recently developed. The purpose of this study was to evaluate the fixation stability of rTKA with this augment in comparison to conventional cemented rTKA.

Eleven pairs of thawed fresh-frozen cadaveric tibias (22 tibias) were potted in custom fixtures. Primary total knee arthroplasty (pTKA) surgery was performed on all tibias. Fixation stability testing was conducted using a three-stage eccentric loading protocol. Static eccentric (70% medial/ 30% lateral) loading of 2100 N was applied to the implants before and after subjecting them to 5×103 loading cycles of 700 N at 2 Hz using a joint motion simulator. Bone-implant micromotion was measured using a high-resolution optical system. The pTKA were removed. The proximal tibial bone defect was measured. One tibia from each pair was randomly allocated to the experimental group, and rTKA was performed with a titanium augment printed using selective laser melting. The contralateral side was assigned to the control group (revision with fully cemented stems). The three-stage eccentric loading protocol was used to test the revision TKAs. Independent t-tests were used to compare the micromotion between the two groups.

After revision TKA, the mean micromotion was 23.1μm ± 26.2μm in the control group and 12.9μm ± 22.2μm in the experimental group. There was significantly less micromotion in the experimental group (p= 0.04). Prior to revision surgery, the control and experimental group had no significant difference in primary TKA micromotion (p= 0.19) and tibial bone loss (p= 0.37).

This study suggests that early fixation stability of revision TKA with the novel 3D printed titanium augment is significantly better then the conventional fully cemented rTKA. The early press-fit fixation of the augment is likely sufficient for promoting bony ingrowth of the augment in vivo. Further studies are needed to investigate the long-term in-vivo fixation of the novel 3D printed augment.

Introduction

Stress shielding of bone around the stem components of total shoulder replacement (TSR) implants can result in bone resorption, leading to loosening and failure. Titanium is an ideal biomaterial for implant stems; however, it is much stiffer than bone. Recent advances in additive manufacturing (AM) have enabled the production of parts with complex geometries from titanium alloys, such as hollow or porous stems. The objective of this computational study is to determine if hollow titanium stems can reduce stress shielding at the proximal humerus. We hypothesize that hollow TSR implant stems will reduce stress shielding in comparison with solid stems and the inner wall thickness of the hollow stem will be a design parameter with a direct effect on bone stresses.

Introduction

Ligament reconstruction following knee soft tissue injuries, such as posterior cruciate ligament (PCL) tears, aim to restore normal joint function and motion; however, persistant pathomechanical joint behavior indicates that there is room for improvement in current reconstruction techniques. Increased attention is being directed towards the roles of secondary knee stabilizers, in an attempt to better understand their contributions to kinematics of knees. The objective of this study is to characterize the relative biomechanical contributions of the posterior oblique ligament (POL) and the deep medial collateral ligament (dMCL) in PCL-deficient knees. We hypothesized that, compared with the POL, the dMCL would have a more substantial role in stabilizing the medial side of the knee, especially in flexion (slack POL).

As a treatment for end-stage elbow joint arthritis, total elbow replacement (TER) results in joint motions similar to the intact joint; however, bearing wear, excessive deformations and/or early fracture may necessitate early revision of failed implant components. Compared to hips, knees and shoulders, very little research has been focused on the evaluation of the outcomes of TER, possible failure mechanisms and the development of optimal designs. The current study aims to develop computational models of TER implants in order to analyze implant behaviour; considering contact stresses, plastic deformations and damage progression.

A geometrical model of a TER assembly was developed based on measurements from a Coonrad-Morrey TER implant (Zimmer, Inc., Warsaw, IN). Ultra high molecular weight polyethylene (UHMWPE) nonlinear elasto-plastic material properties were assigned to the humeral and ulnar bushings. A frictional penalty contact formulation with a coefficient of friction of 0.04 was defined between all of the surfaces of the model to take into account every possible interaction between different implant components in vivo. The loading scenario applied to the model includes a flexion-extension motion, a joint force reaction with variable magnitude and direction and a time varying varus-valgus (VV) moment with a maximum magnitude of 13 N.m, simulating a chair-rise scenario as an extreme loading condition. An explicit dynamic finite element solver was used (ABAQUS Explicit, Dassault Systèmes, Vélizy-Villacoublay, France), due to improved capabilities when performing large deformation analyses. Model results were compared directly with corresponding experimental data. Experimental wear tests were performed on the abovementioned implants using a VIVO (AMTI, Watertown, MA) six degree-of-freedom (6-DOF) joint motion simulator apparatus. The worn TER bushings were scanned after the test using micro computed tomography (µCT) imaging techniques, and reconstructed as 3D models. Comparisons were made based on the sites of damage and deformed geometries between the numerical results and experimental test data. In addition to that, parametric geometrical models were developed using worn geometry of the retrievals in order to account for primary wear and deformations while simulating long-term contact stress and secondary damage progression on the bushings (Fig. 1).

Contact pressure distributions on the humeral and ulnar bushings correlate with the sites of damage as represented by the µCT data and gross observation of clinical retrievals. Furthermore, deformation patterns and kinematics of the components are in good agreement with the experimental results (Fig.2). Excessive plastic deformations are evident in both the numerical and the experimental results close to the regions with high contact pressures. Simulating parametric initially-worn geometries results in the formation of secondary damage zones, as well as redistribution of contact stresses and contact locations (Fig. 3).

The results demonstrate UHMWPE bushing damage due to different loading protocols. Numerical results demonstrate strong agreement with experimental data based on the location of deformation and creep on bushings and exhibit promising capabilities for predicting the damage and failure mechanisms of TER implants.

For figures/tables, please contact authors directly.

Introduction

Hemiarthroplasty is a treatment option for comminuted fractures and non-unions of the distal humerus. Unfortunately, the poor anatomical fit of off-the-shelf distal humeral hemiarthroplasty (DHH) implants can cause altered cartilage contact mechanics. The result is reduced contact area and higher cartilage stresses, thus subsequent cartilage erosion a concern. Previous studies have investigated reverse-engineered DHH implants which reproduce the shape of the distal humerus bone or cartilage at the articulation, but still failed to match native contact mechanics. In this study, design optimization was used to determine the optimal DHH implant shape. We hypothesized that patient-specific optimal implants will outperform population-optimized designs, and both will optimize simple reverse-engineered designs.

Introduction

The intrinsic constraint of a total knee replacement (TKR) implant system is considered an important characteristic which plays a large role in determining stability following surgery. Established techniques for evaluating the constraint of TKR implants, as described in ASTM F 1223-14, do not necessarily map directly to physiologically relevant loading scenarios where instability can occur, and thus give an incomplete picture of the constraint characteristics of a candidate implant design. Sophisticated joint motion simulators now allow for more physiologically representative joint loading (eg. gait), including the contributions of virtual soft tissues. In this study, we employ a function-based constraint measurement technique for evaluating the kinematics of two TKR designs during gait. Furthermore, we employ simulated soft tissues in order to create three “virtual” knees on which the TKR are tested.

Transforaminal lumbar interbody fusion (TLIF) using an implanted cage is the gold standard surgical treatment for disc diseases such as disc collapse and spinal cord compression, when more conservative medical therapy fails. Titanium (Ti) alloys are widely used implant materials due to their superior biocompatibility and corrosion resistance. A new Ti-6Al-4V TLIF cage concept featuring an I-beam cross-section was recently proposed, with the intent to allow bone graft to be introduced secondary to cage implantation. In designing this cage, we desire a clear pathway for bone graft to be injected into the implant, and perfused into the surrounding intervertebral space as much as possible. Therefore, we have employed shape optimization to maximize this pathway, subject to maintaining stresses below the thresholds for fatigue or yielding.

The TLIF I-beam cage (Fig. 1(a)) with an irregular shape was parametrically designed considering a lumbar lordotic angle of 10°, and an insertion angle of 45° through the left or right Kambin's triangles with respect to the sagittal plane. The overall cage dimensions of 30 mm in length, 11 mm in width and 13 mm in height were chosen based on the dimensions of other commercially available cages. The lengths (l_a, l_p) and widths (w_a, w_p) of the anterior and posterior beams determine the sizes of the cage's middle and posterior windows for bone graft injection and perfusion, so they were considered as the design variables for shape optimization. Five dynamic tests (extension/flexion bending, lateral bending, torsion, compression and shear compression, as shown in Fig. 2(b)) for assessing long term cage durability (10⁷ cycles), as described in ASTM F2077, were simulated in ANSYS 15.0. The multiaxial stress state in the cage was converted to an equivalent uniaxial stress state using the Manson-Mcknight approach, in order to test the cage based on uniaxial fatigue testing data of Ti-6Al-4V. A fatigue factor (K) and a critical stress (σ_cr) was introduced by slightly modifying Goodman's equation and von Mises yield criterion, such that a cage design within the safety design region on a Haigh diagram (Fig. 2) must satisfy K ≤ 1 and σ_cr ≤ S_Y = 875 MPa (Ti-6Al-4V yield strength) simultaneously.

After shape optimization, a final design with l_a = 2.30 mm, l_p = 4.33 mm, w_a = 1.20 mm, w_p = 2.50 mm, was converged upon, which maximized the sizes of the cage's windows, as well as satisfying the fatigue and yield strength requirements. In terms of the strength of the optimal cage design, the fatigue factor (K) under dynamic torsion approaches 1 and the critical stress (σ_cr) under dynamic lateral bending approaches the yield strength (S_Y = 875 MPa), indicating that these two loading scenarios are the most dangerous (Table 1). Future work should further validate whether or not the resulting cage design has reached the true global optimum in the feasible design space. Experimental validation of the candidate TLIF I-beam cage design will be a future focus.

For any figures or tables, please contact authors directly (see Info & Metrics tab above).

Introduction

Lateralizing the center of rotation (COR) of reverse total shoulder arthroplasty (rTSA) has the potential to increase functional outcomes of the procedure, namely adduction range of motion (ROM). However, increased torque at the bone-implant interface as a result of lateralization may provoke early implant loosening, especially in situations where two, rather than four, fixation screws are used. The aim of this study was to utilize finite element (FE) models to investigate the effects of lateralization and the number of fixation screws on micromotion and adduction ROM.

Introduction

Reverse total shoulder arthroplasty (RTSA) can partially restore lost range of motion (ROM). Active motion restoration is largely a function of RTSA joint constraint, limiting impingement, and muscle recruitment; however, it may also be a function of implant design. The aim of this computational study was to examine the effects of implant design parameters, such as neck-shaft (N-S) angle and glenoid lateralization, on impingement-free global circumduction range of motion (GC-ROM). GC-ROM summarizes the characteristically complex, wide-ranging envelope of glenohumeral motion into a single quantity for ease of comparison.

Altered distal radioulnar joint contact (DRUJ) mechanics are thought to cause degenerative changes in the joint following injury. Much of the current research examining DRUJ arthrokinematics focuses on the effect of joint malalignment and resultant degenerative changes. Little is known regarding native cartilage contact mechanics in the distal radioulnar joint. Moreover, current techniques used to measure joint contact rely on invasive procedures and are limited to statically loaded positions. The purpose of this study was to examine native distal radioulnar joint contact mechanics during simulated active and passive forearm rotation using a non-invasive imaging approach.

Testing was performed using 8 fresh frozen cadaveric specimens (6 men: 2 women, mean age 62 years) with no CT evidence of osteoarthritis. The specimens were thawed and surgically prepared for biomechanical testing by isolating the tendons of relevant muscles involved in forearm rotation. The humerus was then rigidly secured to a wrist simulator allowing for simulated active and passive forearm rotation. Three-dimensional (3D) cartilage surface reconstructions of the distal radius and ulna were created using volumetric data acquired from computed tomography after joint disarticulation. Using optically tracked motion data and 3D surface reconstructions, the relative position of the cartilage models was rendered and used to measure DRUJ cartilage contact mechanics.

The results of this study indicate that contact area was maximal in the DRUJ at 10 degrees of supination (p=0.002). There was more contact area in supination than pronation for both active (p=0.005) and passive (p=0.027) forearm rotation. There was no statistically significant difference in the size of the DRUJ contact patch when comparing analogous rotation angles for simulated active and passive forearm motion (p=0.55). The contact centroid moved 10.5±2.6 mm volar along the volar-dorsal axis during simulated active supination. Along the proximal-distal axis, the contact centroid moved 5.7±2.4 mm proximal during simulated active supination.

Using the technique employed in this study, it was possible to non-invasively examine joint cartilage contact mechanics of the distal radioulnar joint while undergoing simulated, continuous active and passive forearm rotation. Overall, there were higher contact area values in supination compared with pronation, with a peak at 10 degrees of supination. The contact centroid moved volarly and proximally with supination. There was no difference in the measured cartilage contact area when comparing active and passive forearm rotation. This study gives new insight into the changes in contact patterns at the native distal radioulnar joint during simulated forearm rotation, and has implications for increasing our understanding of altered joint contact mechanics in the setting of deformity.

Failure of reverse total shoulder arthroplasty (rTSA) due to loosening of the metaglene remains a concern. The metaglene is typically affixed to the glenoid via four peripheral bone screws, and the orientations of these screws can affect the stability of the metaglene. The purpose of this finite element analysis (FEA) study was to investigate whether screw orientations should be considered on a patient-specific basis to maximise early fixation.

Three-dimensional geometries of four scapula specimens were obtained by segmenting from CT data in 3D Slicer. A metaglene and four rigidly attached 4.5 mm diameter, 18 mm long cylinders representing screws, were placed on each reamed glenoid. Each screw was placed at one of four orientations, 15° or 7.5° toward or away from the central axis of the metaglene face, while all others were held in the baseline (BL) configuration, where all screws were perpendicular to the metaglene face. Finite element models were created by meshing with linear tetrahedral elements. Material properties of titanium (E=113.8 GPa, v=0.34) were applied to the metaglene and screws. Cortical bone material properties were considered uniform (E=17.5 GPa, ν=0.3) while cancellous bone material properties were non-uniform and mapped on an element-by-element basis using CT attenuation data. The scapula was fully constrained, and a 252 N superiorly oriented shear force was applied to the inferior portion of the metaglene. Contact was modelled at bone-implant and bone-screw interfaces. Displacements of the metaglene with respect to the glenoid were measured. The orientations of each screw that minimised in-plane displacement were used for specimen-specific (SS) configurations. A global (GL) configuration was also defined based on the averages of SS orientations. FE model-predicted metaglene displacements of the SS, GL, and BL screw configurations were compared using paired t-tests.

The average in-plane metaglene displacements for the SS, GL, and BL configurations were 4.8 ± 1.2, 6.5 ± 3.7, and 5.3 ± 1.5 um, respectively. SS configurations significantly decreased displacements by −0.4 ± 0.3 um (−8.5%, p = 0.024) when compared to BL, but the difference of −1.6 ± 3.1 um (25.3%, p = 0.187) was not significant when compared to the GL configuration.

In general, the SS configurations resulted in smaller metaglene displacements than the GL configurations, however the difference was not statistically significant. In one specimen, the GL configuration resulted in abnormally large displacements. These results indicate that, while on average, patient-specific orientations won't yield significantly greater fixation than global configurations; non-patient-specific configurations can, in some cases, yield poor results. Therefore, to ensure optimal fixation for all patients, screw orientations should be considered on a patient-specific basis.