Methods

Computed tomography (CT) scans of 10 humeri were processed in Mimics 17 to create three-dimensional (3D) cortical and cancellous solid bone models. Endoprostheses of different lengths manufactured by Stryker were modelled using Solidworks 2020. The FE models were divided into four groups namely group A consisting of the intact humerus and groups B, C and D composed of humerus-prosthesis assemblies with a body length of 40, 100 and 120 mm respectively and were meshed using linear 4-noded tetrahedral elements in 3matic 13. The models were then imported into Abaqus CAE 6.14. Isotropic linear elastic behaviour with an elastic modulus of 13400, 2000 and 208 000 MPa were assigned to the cortical bone, cancellous bone and prosthesis respectively and a Poisson's ratio of 0.3 was assumed for each material. To represent the lifting of heavy objects and twisting motion, a tensile load of 200 N for axial loading and a 5 Nm torsional load for torsional loading was applied separately to the elbow joint surface with the glenohumeral joint fixed and with all contact interfaces defined as fully bonded. A comparative analysis against literature was performed to validate the intact model. Statistical analysis of the peak von Mises stress values collected from predicted stress contour plots was performed using a one-way repeated measure of analysis of variance (with a Bonferroni post hoc test) using SPSS Statistics 26. The average change in stress of the resected models from the intact state were then determined.