Introduction. Density-modulus relationships are often used to map the mechanical properties of bone based on CT- intensity in finite element models (FEMs). Although these relationships are thought to be site-specific, relationships developed for alternative anatomic locations are often used regardless of bone being modeled. Six relationships are commonly used in finite element studies of the shoulder; however, the accuracy of these relationships have yet to be compared. This study compares each of these six relationships ability to predict apparent strain energy density (SED. app. ) in trabecular bone cores from the glenoid. Methods. Quantitative-CT (QCT) (0.625 mm isotropic voxels), and µ-CT scans (0.032 mm isotropic voxels) were obtained for fourteen cadaveric scapulae (7 male, 7 female). Micro finite element models (µ-FEMs) were created from 98 virtual ‘cores’ using direct conversion to hexahedral elements. Two µ-FEM cases were considered: homogeneous tissue modulus of 20 GPa, and heterogeneous tissue modulus scaled by CT intensity of the µ-CT images (196 models). Each µ- FEM model was compressively loaded to 0.5% apparent strain and apparent strain energy density (SED. app. ) was calculated. Additionally, each of the six density-modulus relationships were used to map heterogeneous material properties to co- registered QCT-derived models (588 models in total). The loading and boundary conditions were replicated in the QCT-FEMs and the SED. app. was calculated and compared to the µ-FEM SED. app. To account for more samples than donors, restricted maximum likelihood estimation (REML) linear regression compared µ-FEM SED. app. and QCT-FEM SED. app. for each relationship. Results. When considering comparisons between QCT-FEMs and µ-FEMs with a homogeneous tissue modulus, near absolute statistical agreement (Y=X) was observed between the µ-FEMs and the QCT-FEMs using the Morgan et al. (2003) pooled relationship. Not surprisingly, due to the similarity between the two relationships, the Gupta & Dan (2004) and Carter and Hayes (1977) models showed near identical REML linear regression fit parameters. All relationships other than the Morgan et al. (2003) pooled relationship, greatly underestimated the µ-FEM apparent strain energy density (SED. app. ) when considering a homogeneous tissue modulus in the µ-FEMs. The same result with the pooled relationship did not hold true when heterogeneous tissue modulus was considered in the µ-FEMs. The Büchler et al., (2002) relationship most accurately predicted the SED. app. for this comparison. Interestingly, the Gupta & Dan (2004) and Carter and Hayes (1977) relationships again showed near identical REML linear regression fit parameters. DISCUSSION. This study compared six common density-modulus relationships used to map mechanical properties of bone in shoulder FE studies. It was found that when considering a homogeneous tissue modulus for µ-FEMs, relationships pooled from alternative anatomic locations may accurately predict the mechanical properties of glenoid trabecular bone. However, when considering a heterogeneous tissue modulus, this did not hold true. Further studies to determine if these relationships can be translated to whole bones may provide insight into the predictive capabilities of using pooled density-modulus equations in the mapping of mechanical properties in future FEMs of the shoulder

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

Introduction. Non-union is debilitating, costly and affects 2–8% of intramedullary fixed fractures. Clinical data suggest that percutaneous interfragmentary screws offer a less invasive alternative to exchange nailing. This study aimed to assess their efficiency with biomechanical analyses. Materials and Methods. A tibia was prepared for finite element analysis by creating a fracture of AO classification 42A2b, prior to reaming and insertion of an intramedullary nail. A callus was modelled as granulation tissue and gait loads were applied. The model was validated against published data and with sensitivity studies. The effects of weightbearing, fracture gap and angle, percutaneous screws and exchange nailing were compared through quantification of interfragmentary motion and strain, with the latter used to gauge healing performance via mechano-regulation theory. Results. Axial interfragmentary motion increased with increasing weightbearing, however, shear decreased at 25–50% weightbearing, leading to superior healing performance. Fracture gap had minimal effect on axial motion, but larger gaps gave greater shear, compromising healing. Elevated fracture obliquity culminated in more shear and inferior healing. Exchange nailing reduced axial motion by ∼30%, but had little effect on shear. Conversely, percutaneous screws had negligible effect on axial motion, but reduced shear by ∼15%, with three screws having a similar net effect on healing as exchange nailing from 10 to 11mm. Conclusions. This study provides new insight into fracture healing biomechanics and discovered that partial weightbearing, less oblique fractures and percutaneous screws all reduce shear, enhancing healing

Introduction. On the basis of a proposal by Noble, the marrow cavity form can be classified into three categories: stovepipe, normal, and champagne-fluted. In the present study, three typical finite element femoral models were created using CT data based on Noble's three categories. The purpose was to identify the relationship of stress distribution of the surrounding areas between femoral bone marrow cavity form and hip stem. The results shed light on whether the distribution of the high-stress area reflects the stem design concept. In order to improve the results of THA, researchers need to consider the instability of a stem design based on the pressure zone and give feedback on future stem selection. Methods. To develop finite element models, two parts (cortical bone and stem) were constructed using four-node tetrahedral elements. The model consisted of about 40,000 elements. The material characteristics were defined by the combination of mass density, elastic coefficient, and Poisson's ratio. Concerning the analysis system, HP Z800 Workstation(HP, Japan) was used as hardware and LS-DYNA Ver. 971 (Livermore Software Technology Corporation, USA) as software. The distal end of the femur was constrained in all directions. On the basis of ISO 7206 Part 4,8 that specifies a method of endurance testing for joint prostheses, the stem was tilted 10°, and a 500 N resultant force in the area around the hip joint was applied to the head at an angle of 25° with the long axis. Automatic contact with a consideration of slip was used. Von Mises stress during a 1.0 s period after loading was analyzed, and stress distribution in the stem and its maximum value were calculated. Result. The maximum stress at marrow cavity form of normal was shown to be 72 MPa. The stress of champagne-fluted was evenly distributed from proximal to distal, and the maximum stress was 67 MPa. For stovepipe, the maximum proximal stress was shown to be 120 MPa; moreover, stress concentration was observed. Discussion. The design concept for a Zweymüller-type stem can distribute load across a wide range of cortical bone from the middle position to the distal femur. It is determined using this concept that a wide range of stress was absorbed at the middle position and distal femur in the champagne-fluted and normal cases. On the other hand, the contact pressure zone of stovepipe could not meet the expected level at the distal femur. The method of this research involves controlling the stress conditions within the stem design. At this point, it is considered possible for the stability of various stem designs to be predicted and the stability to be assessed positively. On the basis of Noble's categories, three types of finite element model were made, and stress distribution measurement and finite element analyses were performed. The results indicate that Zweymüller stem has clinical validity for securing force in the champagne-fluted and stovepipe types from the stress distribution

INTRODUCTION. In theory, Finite Element Analysis (FEA) is an attractive method for elucidating the mechanics of modular implant junctions, including variations in materials, designs, and modes of loading. However, the credence of any computational model can only be established through validation using experimental data. In this study we examine the validity of such a simulation validated by comparing values of interface motion predicted using FEA with values measured during experimental simulation of stair-climbing. MATERIALS and METHODS. Two finite element models (FEM) of a modular implant assembly were created for use in this study, consisting of a 36mm CoCr femoral head attached to a TiAlV rod with a 14/12 trunnion. Two head materials were modelled: CoCr alloy (118,706 10-noded tetrahedral elements), and alumina ceramic (124,710 10-noded tetrahedral elements). The quasi-static coefficients of friction (µ. s. ) of the CoCr-TiAlV and Ceramic-TiAlV interfaces were calculated from uniaxial assembly (2000N) and dis-assembly experiments performed in a mechanical testing machine (Bionix, MTS). Interface displacements during taper assembly and disassembly were measured using digital image correlation (DIC; Dantec Dynamics). The assembly process was also simulated using the computational model with the friction coefficient set to µ. s. and solved using the Siemens Nastran NX 11.0 Solver. The frictional conditions were then varied iteratively to find the value of µ providing the closest estimate to the experimental value of head displacement during assembly. To validate the FEA model, the relative motion between the head and the trunnion was measured during dynamic loading simulating stair-climbing. Each modular junction was assembled in a drop tower apparatus and then cyclically loaded from 230–4300N at 1 Hz for a total of 2,000 cycles. The applied load was oriented at 25° to the trunnion axis in the frontal plane and 10° in the sagittal plane. The displacement of the head relative to the trunnion during cyclic loading was measured by a three-camera digital image correlation (DIC) system. The same loading conditions were simulated using the FEA model using the optimal value of µ derived from the initial head assembly trials. RESULTS. For both head materials, the predicted values of axial displacement of the head on the trunnion closely approximated the measured values derived from DIC measurements, with differences of −0.17% to +6.5%, respectively. Larger differences were calculated for individual components of motion for the stair climbing activity. However, the predicted magnitude of interface motion was still within 10% of the observed values, ranging from −7% to −5%. CONCLUSIONS. Our simulations closely approximated physical testing using complex loading, coming within 7% of the target values. By generating a validated computational model of a modular junctions with varying head materials, we will be able to simulate additional activities of daily living to determine micromotion and areas of peak pressure and contact stresses generated. For any figures or tables, please contact authors directly

Introduction. On the basis of a proposal by Noble, the marrow cavity form can be classified into three categories: normal, champagne-fluted and stovepipe. In the present study, three typical finite element femoral models were created using CT data based on Noble's three categories. The purpose was to identify the relationship of stress distribution of the surrounding areas between femoral bone marrow cavity form and hip stems. The results shed light on whether the distribution of the high-stress area reflects the stem design concept. In order to improve the results of THA, researchers need to consider the instability of a stem design based on the stress distributioin and give feedback on future stem selection. Methods. As analyzing object, we selected SL-PLUS and BiCONTACT stems. To develop finite element models, two parts (cortical bone and stem) were constructed using four-node tetrahedral elements. The model consisted of about 60,000 elements. The material characteristics were defined by the combination of mass density, elastic coefficient, and Poisson's ratio. Concerning the analysis system, HP Z800 Workstation was used as hardware and LS-DYNA Ver. 971 as software. The distal end of the femur was constrained in all directions. On the basis of ISO 7206 Part 4,8 that specifies a method of endurance testing for joint prostheses, the stem was tilted 10°, and a 1500 N resultant force in the area around the hip joint was applied to the head at an angle of 25° with the long axis. Automatic contact with a consideration of slip was used. Result. The maximum stress on femur implanted a SL-PLUS with marrow cavity form of normal, champagne-fluted and stovepipe were shown to be 90MPa, 90MPa and 45MPa. The maximum stress on a BiCONTACT with marrow cavity form of normal, champagne-fluted and stovepipe were shown to be 45MPa, 90MPa and 15MPa. Discussion. The design concept for aZweymüller-type stem can distribute load across a wide range of cortical bone from the middle position to the distal femur. It is determined using this concept that a wide range of stress was absorbed at the middle position and distal femur in the champagne-fluted and normal cases. On the other hand, the contact pressure zone of stovepipe could not meet the expected level at the distal femur. The method of this research involves controlling the stress conditions within the stem design. At this point, it is considered possible for the stability of various stem designs to be predicted and the stability to be assessed positively. On the basis of Noble's categories, three types of finite element model were made, and stress distribution measurement and finite element analyses were performed. The results indicate that Zweymüller stem has clinical validity for securing force in the champagne-fluted and stovepipe types from the stress distribution

Introduction. Relative motion at the modular head-neck junction of hip prostheses can lead to severe surface damage through mechanically-assisted corrosion. One factor affecting the mechanical performance of modular junctions is the frictional resistance of the mating surfaces to relative motion. Low friction increasing forces normal to the head-neck interface, leading to a lower threshold for slipping during weight-bearing. Conversely, a high friction coefficient is expected to limit interface stresses but may also allow uncoupling of the interface in service. This study was performed to examine this trade-off using finite element models of the modular head-neck junction. Methods. A finite element model (FEM) of the trunnion/ head assembly of a total hip prosthesis was initially created and experimentally validated. CAD models of a stem trunnion (taper size: 12/14mm) and a prosthetic femoral head (diameter: 28mm) were discretized into elements for finite element analysis (FEA). The trunnion (Ti6Al4V) was modelled with a hexahedral mesh (33,648 elements) and the femoral head (CoCrMo) with a tetrahedral mesh (51,182 elements). A friction-based sliding contact interface was defined between the mating surfaces. The model was loaded in 2 stages: (i) an assembly load of 4000N applied along the trunnion axis, and (ii) 500N applied along the trunnion axis in combination with a torque of 10Nm. A linear static solution was set up using Siemens NX-Nastran solver. Multiple simulations were executed by modulating the frictional coefficient at the taper-bore interface from 0.05 to 0.15 in increments of 0.01, the coefficient of 0.1 serving as the control case (Swaminathan and Gilbert, 2012). Results. The vertical and tangential displacements of the nodes on the taper of the trunnion relative to the femoral head demonstrated a strong inverse dependence upon the coefficient of friction at the interface (Fig. 1). A similar trend was observed with respect to the peak interface pressure (Fig. 2). The peak von Mises stress, however, increases with increasing coefficient of friction (Fig. 2). A Fisher's R to Z correlation test was performed on each output variable to determine its correlation with coefficient of friction. The coefficient of friction correlated significantly (p<0.0001) with both tangential displacement (r = −0.990) and vertical displacement (r = −0.974). Peak von Mises stress (r = 0.995) and peak contact pressure (r = −0.984) were also found to be significantly (p<0.0001) correlated to the coefficient of friction. Discussion. A higher coefficient of friction at the taper-bore interface led to lower contact pressure and sliding at the modular junction. However, higher coefficients of friction also led to increased von Mises stresses within the bore and the trunnion increasing the risk of yielding and fatigue failure. The current results strongly indicate that factors affecting the frictional coefficient at the interface likely influence the occurrence of and severity of mechanically-assisted corrosion in THA. Significance. The results from this study will help us set tolerances for the interlocking mechanism, identifying the minimum frictional coefficient required to obtain stable implant mechanics

Introduction. Geometric variations of the hip joint can give rise to repeated abnormal contact between the femur and acetabular rim, resulting in cartilage and labrum damage. Population-based geometric parameterisation can facilitate the flexible and automated in silico generation of a range of clinically relevant hip geometries, allowing the position and size of cams to be defined precisely in three dimensions. This is advantageous compared to alpha angles, which are unreliable for stratifying populations by cam type. Alpha angles provide an indication of cam size in a single two-dimensional view, and high alpha angles have been observed in asymptomatic individuals. Parametric geometries can be developed into finite element models to assess the potential effects of morphological variations in bone on soft tissue strains. The aim of this study was to demonstrate the capabilities of our parameterisation research tool by assessing impingement severity resulting from a range of parametrically varied femoral and acetabular geometries. Methods. Custom made MATLAB (MathWorks) and Python codes. [1]. were used to generate bone surfaces, which were developed into finite element models in Abaqus (SIMULIA). Parametric femoral surfaces were defined by a spherical proximal head and ellipse sections through the neck/cam region. This method produced surfaces that were well fitted to bone geometry segmented from CT scans of cam patients and capable of producing trends in results similar to those found using segmented models. A simplified spherical geometry, including the labrum and acetabular cartilage, represented the acetabulum. Femoral parameters were adjusted to define relevant variations in cam size and position. Two radii (small and large cams) and two positions (anterior and superior cams) were defined resulting in four models. Alpha angles of these parametric femurs were measured in an anterior-posterior view and a cross-table lateral view using ImageJ (NIH). A further model was developed using a femur with a medium cam size and position, and the level of acetabular coverage and labrum length were varied. Bones were modelled as rigid bodies and soft tissues were modelled as transversely isotropic linearly elastic materials. With the acetabulum fully constrained in all cases, the femurs were constrained in translation and rotated to simulate flexion followed by internal rotation to cause impingement against the labrum. Results and Discussion. Models generated using the parametric approach showed that potential for tissue damage, indicated through local strain, was not predicted by measured alpha angle, but resulted from cam extent and position as defined by the ellipses. When variations were made to the acetabular rim, an increase in bone coverage had the greatest effect on impingement severity, indicated by strain in the cartilage labral-junction. An increase in labral length increased labral displacement, but had less effect on cartilage-labral strain. Patient specific models currently require full image segmentation, but there is potential to further develop these parametric methods to assess likely impingement severity based on a series of measures of the neck and acetabulum when three-dimensional imaging of patients is available

Mechanical wear and corrosion lead to the release of metal particulate debris and subsequent release of metal ions at the trunnion-taper surface. In order to quantify the amount of volume loss to ultimate locations in the surrounding joint space, finite element analysis of the modular head-stem junction is being carried out. The key purpose being to determine a set of optimum design changes that offer the least material loss at the taper-trunnion junction using optimization algorithms such as the gradient based local search (Sequential Quadratic Programming–SQP) and global search (Non-Dominated Sorting Genetic Algorithm-II–NSGA-II). In a broader sense, the principal goal is to work toward the minimization of wear debris produced in the hip joint, thereby resulting in a longer prosthetic lifetime. A numerical approach that simulates wear in modular hip prostheses with due consideration to the taper-trunnion junction on metal-on-metal contacts is proposed. A quasi-static analysis is performed considering realistic loading stages in the gait cycle, and nonlinear contact analysis is to be employed. The technique incorporates a measured wear rate as an input to the finite element model. The simulation of wear is performed by progressively changing nodal coordinates to simulate the wear loss that occurs during surface interaction. The geometry of the worn surface is updated under gait loading. With a given geometry and gait loading, the linear and volumetric wear increases with the number of gait cycles. The continuous wear propagation is discretized and an approximation scheme known as surrogates is to be developed using Artificial Neural Networks (ANN) to reduce the expensive computational simulations during optimization. The model is employed in the optimization schemes coded in MATLAB and linked to the finite element model developed in ANSYS batch mode. The objective function of the optimization problem is to minimize the volumetric wear at taper-trunnion interface under some constraints. By minimizing the volumetric wear, the chance of failure of modular hip implants is also minimized. The FE model developed to reproduce fretting wear is validated through in vitro wear simulations. The important taper design variables considered to have impact on the fretting corrosion performance include; medial-lateral offset, neck length, taper head diameter, trunnion length and diameter, included angle for the head/neck tapers, angle of mismatch or variation in taper trunnion angle, etc. It is expected from clinical outcomes that increased offset and large taper diameter has serious implications in the fretting corrosion behavior primarily because these variables control the bending stresses and strains along the length of the taper. During cyclic loading of the taper, the higher the strain range, the higher will be the relative micromotion at the point of engagement between the stem and head tapers. This research is carried out with the objective to optimize the effects of these geometrical factors at the mating taper interfaces. The developed models have great potentiality for accurate assessment of wear in a range of metal-on-metal (MoM) hip prostheses at the femoral head taper-trunnion junction while substantially reducing the wear and failure rate of prostheses

Introduction. Circular frame fixation has become a cornerstone of non-union and deformity management since its inception in the 1950s. As a consequence of modularity and heterogenous patient and injury factors, the prediction of the mechanobiological environment within a defect is subject to wide variations in practice. Given these wide range of confounding variables, clinical and cadaveric experimentation is close to impossible and frame constructs are based upon clinician experience. The Finite Element Analysis (FEA) method provides a powerful tool to numerically analyse mechanics. This work aims to develop an FEA model of a tibial defect and predict the mechanical response within the construct. Materials and Methods. The geometry of a tibia was acquired via CT and a series of bone defects were digitally created in the tibial diaphysis. A 4-ring, 10-wire Ilizarov fixator was constructed using 180mm stainless steel rings and 1.8mm stainless steel wires tensioned to 1200N. An axial load (800N) was applied to simulate single leg stance of an 80kg patient. The magnitude of displacement was measured for defects with varying sizes (5–40mm). A numerical analysis was performed in large-strain regime using open-source FEA library (MoFEM). Results. Defect size did not effect displacement, but significantly influenced strain. Measured displacements were 5.72–5.78mm, however strain ranged from 14.5–100%. Moreover, it was found that bone material properties also have no significant impact on the results. Conclusions. Accounting for FEA assumptions, this model predicted a strain environment which was above expected favourable range for bone healing. The addition of graft within the environment is likely to change the mechanobiological environment which warrants further investigation. We plan to develop this model to answer further research questions in the limb reconstruction discipline and validate its accuracy with mechanical data. We believe the presented approach can be a useful tool for investigating the performance circular frames

The Birmingham Hip Mid Head Resection (BMHR) was designed to accommodate patients with lower quality bone in the proximal half of the femoral head. It is a relatively new conservative hip implant with promising early results. Finite element modelling may provide an insight into mid-term results. A cadaveric femur was CT scanned and 3D geometry of the intact femur constructed. The correctly sized BMHR implants (with and without visual stop) were positioned and these verified by a surgeon; hence constructing the post-operative models. Walking loads were applied and contact surfaces defined. Stress analyses were performed using the finite element method and contact examined. Also, a strain-adaptive bone remodelling analysis was run using 45% gait hip loading data. Virtual DEXA images were computed and were analysed in seven regions of the bone surrounding the implants. The BMHR was found to be mechanically stable with all surfaces indicating micromotion less than the critical 150 microns. Stress distribution was similar to the intact femur, with the exception of the head-neck region where some stress/strain shielding occurs. This is mirrored in the bone remodelling results, which show some bone resorption in this region. The visual stop, which is designed to ensure that the stem is not overdriven during implantation, did not affect the stress/strain results; only on a very local scale. There is minimal data available in the literature regarding conservative hip implants and no data regarding the BMHR. This study is the first to look at the mechanical response of the bone to this implant

INTRODUCTION. Mechanical overloading of the knee can occur during activities of daily living such as stair climbing, jogging, etc. In this finite element study we aim to investigate which parameters could detrimentally influence peri-implant bone in the tibial reconstructed knee. Bone quality and patient variables are potential factors influencing knee overloading (Zimmerman 2016). METHODS. Finite element (FE) models of post-mortem retrieved tibial specimens (n=7) from a previous study (Zimmerman 2016) were created using image segmentation (Mimics Materialise v14) of CT scan data (0.6 mm voxel resolution). Tibial tray and polyethylene inserts were recreated from CT data and measurements of the specimens (Solidworks 2015). Specimens with varying implant geometry (keel/pegged) were chosen for this study. A cohesive layer between bone and cement was included to simulate the behavior of the bone–cement interface using experimentally obtained values. The FE models predict plasticity of bone according to Keyak (2005). Models were loaded to 10 body weight (BW) and then reduced to 1 BW to mimic experimental measurements. Axial FE bone strains at 1 BW were compared with experimental Digital Image Correlation (DIC) bone strains on cut sections of the specimens. After validation of the FE models using strain data, models were rotated and translated to the coordinate system defined in Bergmann (2014). Four loading cases were chosen – walking, descending stairs, sitting down and jogging. Element strains were written to file for post-processing. The bone in all FE models was divided into regions of equal thickness (10 mm) for comparison of strains. RESULTS. Results are shown for two specimens at present. Strain-maps of the specimen cut section compare reasonably well with FE cutting-plane strains. The FE models however show some regions of high strain in certain locations which do not correspond with the experimental results (Figure 1). Plasticity predicted by the models at 10 BW is shown in Figure 2. Median bone strains for two loading cases are shown as a function of distance below the tibial tray in Figure 3. This figure shows that specimen 1 is less likely to be overloaded during jogging when compared with specimen 2. Both specimens remain below the 7300 με threshold for compressive yield. DISCUSSION. Using functioning knee replacement tibial specimens, we study which factors influence bone overloading. Validation using DIC strain measurements is challenging due to the large plasticity regions predicted by the material model used here. The present results were obtained using plasticity relationships from Keyak (2005) for the proximal femur. To further improve on our results, plasticity-bone density relationships for the proximal tibia (Keyak 1996) will be included. Proximal tibial bone has been shown to be stiffer than femoral bone (Morgan 2003). Despite these limitations, FE models provide valuable information on the risk of overloading during daily living activities. The study will be expanded to include an analysis of implant geometry, bone quality and other loading cases. For any figures or tables, please contact the authors directly

Distal radius fractures are the most common osteoporotic fractures among women. The treatment of these fractures has been shifting from a traditional non-operative approach to surgery, using volar locking plate (VLP) technology. Surgery, however, is not without risk, complications including failure to restore an anatomic reduction, fracture re-displacement, and tendon rupture. The VLP implant is also marked by bone loss due to stress-shielding related to its high stiffness relative to adjacent bone. Recently, a novel internal, composite-based implant, with a stiffness less than the VLP, was designed to eradicate the shortcomings associated with the VLP implant. It is unclear, however, what effect this less-stiff implant will have upon adjacent bone density distributions long-term. The objective of this study was to evaluate the long-term effects of the two implants (the novel surgical implant and the gold-standard VLP) by using subject-specific finite element (FE) models integrated with an adaptive bone formation/resorption algorithm. Specimen: One fresh-frozen human forearm specimen (female, age = 84 years old) was imaged using CT and was used to create a subject-specific FE model of the radius. Finite element modeling: In order to simulate a clinically relevant (unstable) fracture of the distal radius, a wedge of bone was removed from the model, which was approximately 10 mm wide and centered 20 mm proximal to the tip of the radial styloid. Bone remodeling algorithm: A strain-energy density (SED) based bone remodeling theory was used to account for bone remodeling. With this approach, bone density decreased linearly when SED per bone density was less than 67.5 µJ/g and increased when it was more than 232.5 µJ/g. When it was in the lazy zone (67.5 to 232.5 µJ/g), no changes in density occurred. Boundary conditions: A 180 N quasi-static force representing the scaphoid, and a 120 N quasi-static force representing the lunate was applied to the radius. The midshaft of the radius was constrained. FE outcomes: To examine the effects of stress shielding associated with each implant, the long-term changes of bone density within proximal transverse cross-sections of radius were inspected. The regional density analysis focused on three transverse cross-sections. The transverse cross-sections were positioned proximal to the subchondral plate, and were distanced 50 (cross-section A), 57 (cross-section B), and 64 mm (cross-section C) from the subchondral endplate. For both implants in all three cross-sections, cortical bone was reserved completely at the volar side. On the dorsal side, the cortical bone was completely resorbed in the VLP model. In all cross-sections, the averaged resultant density was higher for the “novel implant”. The difference ranged from 33% (cross-section A) to 36% (cross-section C) in favor of the “novel implant”. On average, the density values of the novel implant were 34% higher in transverse cross-sections (A, B, and C). This study showed that the novel implant offered higher density distributions compared to the VLP, which suggests that the novel implant may be superior to the VLP in terms of avoiding stress shielding

Background. Use of a baseplate with a smaller diameter in reverse shoulder arthroplasty has been recommended, especially in patients with a small glenoid or insufficient bony stock due to severe glenoid wear. However, effect of a smaller baseplate on stability of the glenoid component has not been evaluated. The purpose of this study was to determine whether a smaller baseplate (25 mm) is beneficial to the initial primary stability of the glenoid component compared to that with a baseplate of a commonly used size (29 mm) by finite element analysis. Methods. Computed tomography (CT) scans of fourteen scapulae were acquired from cadavers with no apparent deformity or degenerative change. Glenoid diameter corresponding to the diameter of the inferior circle of glenoid was measured using a caliper and classified into the small and large glenoid groups based on 25mm diameter. CT slices were used to construct 3-dimensional models with Mimics (Materialise, Leuven, Belgium). A corresponding 3D Tornier Aequalis® Reversed Shoulder prosthesis model was generated by laser scanning (Rexcan 3D Laser Scanner, Solutionix, Seoul, Korea). Glenoid components with 25mm and 28mm diameter of the baseplate were implanted into the scapular of small and large glenoid group, respectively. Finite element models were constructed using Hypermesh 11.0 (Altair Engineering, Troy, MI, USA) and a reverse engineering program (Rapidform 3D Systems, Inc., Rock Hill, SC, USA). Abaqus 6.10 (Dassault Systemes, Waltham, MA) was used to simulate 30. o. , 60. o. , and 90. o. glenohumeral abduction in the scapular plane. Single axial loads of 686N (1 BW) at angles of 30. o. , 60. o. , and 90. o. abduction were applied to the center of the glenosphere parallel to the long axis of the humeral stem. Relative micromotion at the middle and inferior thirds bone–glenoid component interface, and distribution of bone stress under the glenoid component and around the screws were analyzed. Wilcoxon's rank-sum test was used for statistical comparison and p < 0.05 was considered as a minimum level of statistical significance. Results. In small glenoid group, micromotion at the middle and inferior thirds of the glenoid-glenosphere interface at angles of 30. o. and 60. o. abduction were significantly greater in the 29mm baseplate than in the 25mm baseplate. There was no significant difference in micromotion at angle of 90. o. abduction between 25mm and 29mm baseplate. In large glenoid group, there was no statistically significant difference in micromotion between 25mm and 29mm baseplate at all angles of abduction. In small glenoid group, maximum bone stress was measured at the point of cortical engagement of the inferior screw and was statistically greater in the 29mm baseplate than in the 25mm baseplate. In large glenoid group, there was no statistically significant difference of maximum bone stress around the inferior screw between 25mm and 29mm baseplates. Conclusions. Use of a baseplate with a smaller diameter (25 mm) in reverse shoulder arthroplasty is suitable for improving the primary stability of the glenoid component, especially in small glenoid

Introduction. Total hip replacement is an established surgical procedure done to alleviate hip pain due to joint diseases. However, this procedure is avoided in yonger patients with higher functional demands due to the potential for early failure. An ideal prosthesis will have have a high endurance against impact loading, with minimal micromotion at the bone cement interface, and a reduced risk of fatigue failure, with a favourable stress distribution pattern in the femur. We study the effect of varying the material properties and design element in a standard cemented total hip using Finite Element Analysis. Methods. A patient-specific 3D model of femur will be constructed from CT scan data, while a Summit® Cemented Hip System (DePuy Orthopedic) will be used to as a control for comparative evaluation. We vary the material stiffness of different parts of the prosthesis(see Fig.1) to formulate a design concept for a new total hip prosthesis design; and use Finite Element Method to predict the micromotion of the hip prosthesis at the bone cement interface, as well as the stress distribution in the the femur. Result. Validation of computational protocol was being done by comparing the principal maximum strain of the femoral cortex along the diaphysis, and the amount of deflection, with published literature, similarly, contact modelling validation was also done. Model 1–4 induced lower peak Von Mises stress in the cement, which takes a much lower value than any of the cement mechanical limits postulated. Therefore, the risk of cement failure is greatly reduced in Model 1–4. However, the effect of varying stiffness in different regions is not significant in terms of load transmission to the cement. Micromotion at the bone-cement interface was studied via two approaches: Peak micromotion at the bone cement interface; and the micromotion data at 12 Regions of Interest (ROI)s. Both results showed that model 2 and 3 are capable of reducing micromotion at bone-cement interface, in comparison with the Summit® Cemented Hip System. By comparing the Von Mises Stress distribution in the proximal femur; model 1 is found to result in a significantly reduced stress shielding effect, while model 2–4 are also favourable in comparison to the standard Summit® prosthesis in terms of stress distribution in the femur. Figure 2 shows the effects of the performance of model 1–4, presented as percentage difference from the Summit® prosthesis. Model 1 is unfavourable, despite its favourable stress distribution, because its peak and overall micromotion at the bone-cement interface is greatly increased. Conclusion. Model 2 and 3 have favourable design elements. They both have reduced micromotion at the bone-cement interface; and a favourable stress distribution in the femur. Further refining and testing of model 2 and 3 should done, as these models may provide information which may be useful in improving the performance of the current range of total hip replacement prostheses

Squeaking ceramics bearing surfaces have been recently recognised as a problem in total hip arthroplasty. The position of the acetabular cup has been alluded to as a potential cause of the squeaking, along with particular combinations of primary stems and acetabular cups. This study has used the finite element method to investigate the propensity of a new large diameter preassembled ceramic acetabular cup to squeaking due to malpositioning. A verified three-dimensional FE model of a cadaveric human pelvis was developed which had been CT scanned, and the geometry reconstructed; this was to be used to determine the behaviour of large diameter acetabular cup system with a thin delta ceramic liner in the acetabulum. The model was generated using ABAQUS CAE pre-processing software. The bone model incorporated both the geometry and the materials properties of the bone throughout based on the CT scan. Finite element analysis and bone material assignment was performed using ABAQUS software and a FORTRAN user subroutine. The loading applied simulated edge loading for rising from a chair, heel-strike, toe off and stumbling. All results of the analysis were used to determine if the liner separated from the shell and if the liner was toggling out of the shell. The results were also examined to see if there was a propensity for the liner to demobilise and vibrate causing a squeaking sound under the prescribed loading regime. This study indicates that there is a reduction in contact area between the ceramic liner and titanium shell if a patient happens to trip or stumble. However, since the contact between the liner and the shell is not completely lost the propensity for it to squeak is highly unlikely

Introduction. In total hip arthroplasty, press-fit anchorage is one of the most common fixation methods for acetabular cups and mostly ensures sufficient primary stability. Nevertheless, implants may fail due to aseptic loosening over time, especially when the surrounding bone is affected by stress-shielding. The use of acetabular cups made of isoelastic materials might help to avoid stress-shielding and osteolysis. The aim of the present numerical study was to determine whether a modular acetabular cup with a shell made of polyetheretherketone (PEEK) may be an alternative to conventional titanium shells (Ti6Al4V). For this purpose, a 3D finite element analysis was performed, in which the implantation of modular acetabular cups into an artificial bone stock using shells made of either PEEK or Ti6Al4V, was simulated with respect to stresses and deformations within the implants. Methods. The implantation of a modular cup, consisting of a shell made of PEEK or Ti6Al4V and an insert made of either ceramic or polyethylene (PE), into a bone cavity made of polyurethane foam (20 pcf), was analysed by 3D finite element simulation. A two-point clamping cavity was chosen to represent a worst-case situation in terms of shell deformation. Five materials were considered; with Ti6Al4V and ceramic being defined as linear elastic and PE and PEEK as plastic materials. The artificial bone stock was simulated as a crushable foam. Contacts were generated between the cavity and shell (μ = 0.5) and between the shell and insert (μ = 0.16). In total, the FE models consisted of 45,282 linear hexahedron elements and the implantation process was simulated in four steps: 1. Displacement driven insertion of the cup; 2. Relief of the cup; 3. Displacement driven placement of the insert; 4. Load driven insertion of the insert (maximum push-in force of 500 N). The FE model was evaluated with respect to the radial deformations of the shell and insert as well as the principal stresses in case of the ceramic inserts. The model was experimentally validated via comparison of nominal strains of the titanium shells. Results. The maximum radial deformation of the shell made of PEEK was 581 μm (insertion) and 470 μm (relief) and therefore multiple times higher compared to the Ti6Al4V shell (42 μm and 21 μm). As a result, larger deformations occurred at the PE and ceramic inserts in combination with the PEEK shell. Partially, the deformations were above an usual clearance of 100 μm. When the ceramic insert was combined with the shell made of PEEK, maximum principal stresses in the ceramic insert amounted to 30 MPa and were clearly lower than approved bending strength of the ceramic material (948 MPa). Conclusion. The examined acetabular shell made of PEEK was intensively deformed during insertion compared to the geometrically identical Ti6Al4V shell and is therefore not suitable for modular acetabular cups. In future studies it should be clarified to what extent acetabular cups with shells made of carbon fiber reinforced PEEK materials with higher stiffness lead to reduced deformations during the insertion procedure

A common location for radius fracture is the proximal radial head. With the arm in neutral position, the fracture usually happens in the anterolateral quadrant (Lacheta et al., 2019). If traditional surgeries are not enough to induce bone stabilization and vascularization, or the fracture can be defined grade III or grade IV (Mason classification), a radial head prosthesis can be the optimal compromise between bone saving and recovering the “terrible triad”. A commercially available design of radial head prosthesis such as Antea (Adler Ortho, Milan, Italy) is characterized by flexibility in selecting the best matching size for patients and induced osteointegration thanks to the Ti-Por. ®. radial stem realized by 3D printing with laser technique (Figure 1). As demonstrated, Ti-Por. ®. push-out resistance increased 45% between 8 −12 weeks after implantation, hence confirming the ideal bone-osteointegration. Additional features of Antea are: bipolarity, modularity, TiN coating, radiolucency, hypoallergenic, 10° self-aligning. The osteointegration is of paramount importance for radius, in fact the literature is unfortunately reporting several clinical cases for which the fracture of the prosthesis happened after bone-resorption. Even if related to an uncommon activity, the combination of mechanical resistance provided by the prosthesis and the stabilization due to the osteointegration should cover also accidental movements. Based upon Lacheta et al. (2019), after axial compression-load until radii failure, all native specimens survived a compression-load of 500N, while the failure happened for a mean compression force of 2560N. The aim of this research study was to test the mechanical resistance of a radial head prosthesis obtained by 3D printing. In detail, a finite element analysis (FEA) was used to understand the mechanical resistance of the core of the prosthesis and the potential bone fracture induced in the radius with simulated bone- resorption (Figure 2a). The critical level was estimated at the height for which the thickness of the core is the minimum (Figure 2b). Considered boundary conditions:. - Full-length prosthesis plus radius out of the cement block equal to 60mm (Figure 2a);. - Bone inside the cement equal to 60mm (Figure 2b);. - Load inclined 10° epiphysiary component (Figure 2c);. - Radius with physiological or osteoporotic bone conditions;. - Load (concentrated in the sphere simulating full transmission from the articulation) of 500N or 1300N or 2560N. Figure 3 shows the results in terms of maximum stress on the core of the prosthesis and the risk of fracture (Schileo et al., 2008). According to the obtained results, the radial head prosthesis shows promising mechanical resistance despite of the simulated bone-resorption for all applied loads except for 2560N. The estimated mechanical limit for the material in use is 200MPa. The risk of fracture is in agreement with the experimental findings (Lacheta et al. (2019)), in fact bone starts to fail for the minimum reported failure load, but only for osteoporotic conditions. The presented FEA aimed at investigating the behavior of a femoral head prostheses made by 3D printing with simulated bone-resorption. The prosthesis shows to be a skilled solution even during accidental loads. For any figures or tables, please contact the authors directly

Introduction. In complex primary and revision total knee replacement (TKR) the operating surgeon may encounter proximal tibial bone defects. The correct management of such defects is fundamental to both the initial stability and long-term survival of the prosthesis. Cement or metal augments have been used to address some such type II unconstrained defects [1]. Aim. The aim of this finite element (FE) study was to analyse the comparative behaviour of cement and metal based augments and quantify the stresses within these different augments and underlying cancellous bone. Materials and methods. A three-dimensional FE model was constructed from a computer tomography (CT) scan of the proximal tibia using SIMPLEWARE v3.2 image processing software. The tibial component of a TKR was implanted with either a block or wedge-shaped augment made of either metal or cement. The model was axially loaded with a force of 3600N and testing was conducted with both evenly and eccentrically distributed loads. Results. Upon loading the FE model, the von-Mises stresses in the cancellous bone underneath the augments were found to be higher with cement based augments in comparison their metal based counterparts. This was evident with both block and wedge-shaped augments. The FE model demonstrated that compressive stresses within the metal based augments were greater than those within the cement based augments. This was evident with both block and wedge designs. Upon even loading the maximum recorded compressive stresses within the metal augments were 5 times less than the endurance limit of the material [3]. However, the maximum recorded compressive stresses within cement augments were only half the endurance limit of the material [4] and upon eccentric loading compressive stresses in excess of the endurance limit were recorded. Discussion. The FE model has demonstrated that cement based augments undergo a greater deformation when loaded and therefore transfer greater loads to the underlying cancellous bone. This is a result of the inherent flexibility of the cement based augment in comparison to the stiffer metal counterparts. The greater transference of load to cancellous bone with cement based augments may reduce the possibility of stress shielding. However, the compressive stresses within cement based augments are too close to the endurance limit of the material and with uneven loading even exceed it. This would imply that cement based augments are more prone to fatigue failure than their metal counterparts. Conclusion. This FE study supports the use of metal based augments over cement based augments in augmented and revision TKR surgery

Introduction. Varying degrees of posterior glenoid bone loss occurs in patients with end stage osteoarthritis and can result in increased glenoid retroversion. Ultimately, the goal is to correct retroversion to restore normal biomechanics of the glenohumeral joint. The goal of this study was to identify the optimal augmented glenoid design based on finite element model analysis which will provide key insights into implant loosening mechanisms and stability. Materials and Methods. Two different augmented glenoid designs, posterior wedge and posterior step- were created as a computer model by a computer aided design software (CAD). These implants were virtually implanted to correct 20° glenoid retroversion and the different mechanical parameters were calculated including: the glenohumeral contact pressure, the cement stress, the shear stress, and relative micromotions at the bone cement interface. Results. During abduction, high strain was concentrated around the peg and posterior glenoid bone. Strain was noticeably higher in stepped design (1–2%) than the wedged design (0.4–1.2%). Stepped glenoid models sustained 30% and 70% higher stresses than those experienced by the wedged glenoid implant models at two different corrections. Distractions predicted by the stepped designs were found to be at least twice as much as those by the wedged designs. Similarly, in compression values were 1.5–8 magnitudes higher in stepped designs than those of wedged designs. The wedged design, the amount of micromotion was not affected by the size of the augment (8° and 16°). Discussion. Our study showed that the wedged design experienced less stress compared to stepped design with abduction loading. Notably, the wedged design experienced less stress as the size of the wedge increased to correct a more retroverted arthritic glenoid. The step design also had the highest amount of micromotion which ultimately points to increased failures rate and decreased performace