Most glenoid implants rely on centrally located large fixation features to avoid perforation of the glenoid vault in its peripheral regions [1]. Upon revision of such components there may not be enough bone left for the reinsertion of an anatomical prosthesis, resulting in a large cavity that resembles a sink hole. Multiple press-fit small pegs would allow for less bone resection and strong anchoring in the stiffer and denser peripheral subchondral bone [2], whilst producing a more uniform stress distribution and increased shear resistance per unit volume [3] and avoiding the complications from the use of bone cement. This study assessed the best combination of anchoring strength, assessed as the ratio between push in and pull out forces (Pin/Pout), and spring-back, measured as the elastic displacement immediately after insertion, for five different small press-fitted peg configurations (Figure 1, left) manufactured out of UHMWPE cylinders (5 mm diameter and length).

16 specimens for each configuration were tested in two types of Sawbones solid bone substitute: hard (40 PCF, 0.64 g/cm3, worst-case scenario of Pin) and soft (15 PCF, 0.24 g/cm3, worst-case scenario of spring-back and Pout). Two different interference-fits, Ø, were studied by drilling holes with 4.7 mm and 4.5 mm diameter (Ø 0.3 and Ø 0.5, respectively). A maximum Pin per peg of 50 N was defined, in order to avoid fracture of the glenoid bone during insertion of multiple pegs. The peg specimens were mounted into the single-axis screw-driven Instron through a threaded fixture. A schematic of the experimental set up is made available (Figure 1, centre). The peg was pushed in vertically for a maximum of 5 mm at a 1 mm/s rate, under displacement control, recording Pin. The spring-back effect was assessed by switching to load control and reducing the load to zero. The peg was then pulled out at a rate of 1 mm/s, recording Pout. The test profile is depicted in Figure 1 (right).

Average Pout/Pin, spring back (in mm) and force-displacement curves for all 80 specimens tested are shown in Figure 2. These were split into groups according to the type of bone substitute and interference-fit, with the right column showing the average values for the Pin. High repeatability among samples of the same configuration tested is noted. Configurations #1, #3 and #5 all exceed the maximum Pin per peg for at least one type of bone. Configuration #2 has the lowest Pin of all (best thread aspect ratio), followed by configuration #4 (thinner threads). The peg configurations #4 and #2 had the highest Pin/Pout. The peg configurations with lowest spring-back after insertion were configuration #2 and #4. Interference fit of Ø 0.3 mm was shown to reduce Pin below maximum limit of 50 N without great influence in spring-back.

Common post-operative problems in shoulder arthroplasty such as glenoid loosening and joint instability can be reduced by improvements in glenoid design shape, material choice and fixation method [1]. Innovation in shoulder replacement is usually carried out by introducing incremental changes to functioning implants [2], possibly overlooking other successful design combinations.

We propose an automated framework for parametric analysis of implant design in order to efficiently assess different possible glenoid configurations. Parametric variations of reference geometries of a glenoid implant were automatically generated in SolidWorks. The different implants were aligned and implanted with repeatability using Rhino. The glenoid-bone models were meshed in Abaqus, and boundary conditions and loading applied via a custom-made Python script. Finally, another MATLAB script integrated and automated the different steps, extracted and analysed the results.

This study compared the influence of reference shape (keel vs. 2-pegged) and material on the von Mises stresses and tensile and compressive strains of glenoid components with bearing surface thickness and fixation feature width of 3, 4, 5 or 6 mm. A total of 96 different glenoid geometries were implanted into a bone cube (E = 300 MPa, ν = 0.3). Fixed boundary conditions were applied at the distal surface of the cube and a contact force of 1000 N was distributed between the central nodes on the bearing surface. The implants were assigned UHMWPE (E = 1 GPa, ν = 0.46), Vitamin E PE (E = 800 MPa, ν = 0.46), CFR-PEEK (E = 18 GPa, ν = 0.41) or PCU (E = 2 GPa, ν = 0.38) material properties and the bone-implant surface was tied (Figure 1). The von Mises stresses, compressive and tensile strains for the different models were extracted. The influence of design parameters in the mechanical environment of the implant could be assessed. In this particular example, the 95^th percentile values of the tensile and compressive strains induced by modifications in reference shape could be evaluated for all the different geometries simultaneously in form of radar plots. 2-pegged geometries (green) consistently produced lower tensile and compressive strains than the keeled (blue) configurations (Figure 2). Vitamin E PE and PCU glenoids also produced lower maximum von Mises stresses values than CFR-PEEK and UHMWPE designs (Figure 3).

The developed method allows for simple, direct, rapid and repeatable comparison of different design features, material choices or fixation methods by analysing how they influence the mechanical environment of the bone surrounding the implant. Such tool can provide invaluable insight in implant design optimisation by screening through multiple potential design modifications at an early design evaluation stage and highlighting the best performing combinations. Future work will introduce physiological bone geometries and loading, a wider variety of reference geometries and fixation features, and look at bone/interface strength and osteointegration predictions.

Recently finite element studies have incorporated bone remodelling algorithms in an attempt to simulate bone's mechano-adaptation to loading conditions. In order to simplify these analyses, bone is usually considered to be isotropic, which does not explain the directionality of its internal structures; neither the orthotropic properties measured at the continuum level. Furthermore, simplified loading is usually applied to the bone models, which can result in an unrealistic remodelling stimulus. However, free boundary condition modelling of the femoral and pelvic constructs has been shown to produce more physiological stress and strain distributions.

This paper describes the application of a 3D remodelling algorithm (with bone modelled as a strain-adaptive continuum with local orthotropic material properties) to a free boundary model of the femoral construct, where the hip and knee joints, as well as muscles and ligaments crossing the joints were included explicitly. Two load cases were analysed: single leg stance and standing up.

Material properties and directionality distributions were produced for the whole femur, showing good agreement with observed structures from clinical studies. This indicates that the loading conditions modelled correspond to those experienced in vivo. In addition, the impact of the different load cases in bone structure modelling could be compared. Observations of the material properties distribution and orientation for standing up indicate that it promotes changes in bone stiffness in the anterior regions of the femoral neck and cortical shaft and the posterior side of the condyles.

Development of this approach to modelling and bone structure prediction can lead to a better understanding of bone's mechanical behaviour and to the development and public release of orthotropic heterogeneous models for different constructs. These can be applied in many areas of interest in orthopaedic biomechanics, such as the study of bone-implant interfaces, improvement of the currently used surgical tools and techniques and the influence of certain activities in affecting local bone strength and mineralisation.