
==== Front
J Exp Orthop
J Exp Orthop
10.1002/(ISSN)2197-1153
JEO2
Journal of Experimental Orthopaedics
2197-1153
John Wiley and Sons Inc. Hoboken

10.1002/jeo2.70000
JEO270000
Original Paper
Original Paper
Finite element analysis part 1 of 2: Influence of short stem implant polyethylene configuration on glenohumeral joint biomechanics
Nourissat Geoffroy http://orcid.org/0000-0002-2578-3051
1 gnourissat@wanadoo.fr

Housset Victor 1 2
Daudet Jean‐Marie 3
Fradet Léo 4
Bianco Rohan‐Jean 4
Srikumaran Uma 5
1 Groupe Maussins Clinique Maussins Nollet‐Ramsay Santé Paris France
2 Hôpital Henri Mondor Université Paris‐Est Créteil Créteil France
3 FX Shoulder Solutions Viriat France
4 Philomec Inc. Montréal Québec Canada
5 Department of Orthopedic Surgery The Johns Hopkins University Baltimore Maryland USA
* Correspondence Geoffroy Nourissat, Clinique Maussins Nollet‐Ramsay Santé, 67 Rue de Romainville, Paris 75019, France.
Email: gnourissat@wanadoo.fr

19 9 2024
7 2024
11 3 10.1002/jeo2.v11.3 e7000006 7 2024
09 12 2023
20 7 2024
© 2024 The Author(s). Journal of Experimental Orthopaedics published by John Wiley & Sons Ltd on behalf of European Society of Sports Traumatology, Knee Surgery and Arthroscopy.
https://creativecommons.org/licenses/by/4.0/ This is an open access article under the terms of the http://creativecommons.org/licenses/by/4.0/ License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Abstract

Purpose

Stress shielding in short‐stem arthroplasty can cause critical metaphyseal bone loss. If the size and shape of the humeral shaft are important factors, it is unknown whether the shape of the polyethylene component in reverse shoulder arthroplasty (RSA) affects bone stress around or within the stem. We explored the impact of polyethylene shape on humeral and scapular stress distribution using a finite element model.

Methods

We developed a shoulder‐specific finite element model. A defined set of muscle forces was applied to simulate movements. An intact rotator cuff state and a superior deficient rotator cuff state were modelled. We used the FX V135 short stem in three conditions: total shoulder arthroplasty (TSA), and RSA with symmetrical and asymmetrical polyethylene (145°/135°). We measured biomechanical markers related to bone stress for different implant sizes. Joint kinematics and the mechanical behaviour of the implant were compared.

Results

Rupture of the supraspinatus muscle produced a functionally limited shoulder. The placement of an anatomic TSA with an intact rotator cuff restored function similar to that of a healthy shoulder. RSA in the rotator cuff‐deficient shoulder restored function regardless of stem size and polyethylene shape. While stem size had an impact on the stress distribution in the bone and implant, it did not show significant potential for increasing or decreasing overall stress. For the same stem, stress distribution at the humerus is different between TSA and RSA. Polyethylene shape did not alter the transmission of stress to the bone in RSA. Asymmetric polyethylene produced a greater abduction range of motion.

Conclusions

In terms of bone stress distribution, smaller stems seemed more appropriate for TSA, while larger stems may be more appropriate for RSA. Polyethylene shape resulted in different ranges of motion but did not influence bone stress.

Level of Evidence

Diagnostic Tests or Criteria; Level IV.

anatomic total shoulder arthroplasty
bone stress
finite element model
polyethylene shape
reverse shoulder arthroplasty
source-schema-version-number2.0
cover-dateJuly 2024
details-of-publishers-convertorConverter:WILEY_ML3GV2_TO_JATSPMC version:6.4.8 mode:remove_FC converted:19.09.2024
Nourissat, G. , Housset, V. , Daudet, J.‐M. , Fradet, L. , Bianco, R.‐J. & Srikumaran, U. (2024) Finite element analysis part 1 of 2: Influence of short stem implant polyethylene configuration on glenohumeral joint biomechanics. Journal of Experimental Orthopaedics, 11 , e70000. 10.1002/jeo2.70000
==== Body
pmcAbbreviations

RSA reverse shoulder arthroplasty

TSA total shoulder arthroplasty

BACKGROUND

Humeral stress shielding is a major concern in short stem anatomic total shoulder arthroplasty (TSA) and reverse shoulder arthroplasty (RSA) because it can cause critical metaphyseal bone loss and fractures around the stem [15, 17, 40]. Several elements may contribute to stress shielding, including the size and shape of the humeral shaft [16, 30], but it is unknown whether the shape of the polyethylene component (i.e., humeral insert) influences bone stress around or within the humeral stem. Polyethylene shape changes the position of the humeral head relative to the glenoid. This geometrical modification can have an impact on abduction initiation and maximal range of motion. This needs also to be quantified in parallel with bone and implant stress. Finite element models have been used to study the mechanisms of shoulder mobility and stability in both healthy and rotator cuff‐deficient configurations [2, 18, 45]. Concerning RSA, finite element studies have evaluated risk factors for augmented stresses at the interface of humeral components and the risk of bone resorption and stress shielding. Stress distribution is different between TSA and RSA [2, 18, 45]. No study has evaluated the difference of the stress shielding of the same short stem design in TSA and RSA configuration. Furthermore, to our knowledge, no study has evaluated the role of the shape of the polyethylene humeral liner and short stem size in RSA.

Our objective was to explore the impact of polyethylene component shape and stem size on postoperative shoulder biomechanics in short humeral stems with a focus on bone stress and to evaluate if the stress shielding around the stem is the same between anatomic and reverse configuration. Our hypothesis was that both of the tested parameters have different influences on postoperative glenohumeral biomechanics.

MATERIALS AND METHODS

Finite element model description

A finite element model of a glenohumeral joint was created according to the geometric measurements of a patient with no known disorders or deformities (BodyParts3D, Integrated Database Center for Life Sciences, licensed by CC Attribution‐ShareAlike 2.1). The model was morphed to match the 50th percentile bone geometry of a healthy male patient previously described by other authors [27, 36, 39]. The model included trabecular and cortical bone for the humerus, scapula and clavicle; cartilage and labrum of the glenohumeral joint; glenohumeral ligaments; and 12 muscles involved in scapulohumeral movements, including the rotator cuff and deltoid. The ligament and tendon origins and attachments were determined from the literature [14, 48] and confirmed by two independent senior surgeons (Uma Srikumaran and Geoffroy Nourissat) (Figure 1). The regional thickness of the cortical bone wall and cartilage was represented according to published morphometric measurements [10, 23, 24, 34, 38, 41, 43].

Figure 1 Meshed model of the humerus, scapula and clavicle including the 12 muscle bundle elements (1D active springs) and tendinous insertions (2D shell elements).

Comparisons between the model and literature references were made for scapular thickness [9], glenoid version and inclination [28, 32, 35] and glenoid cavity and humeral head curvatures [49]. These comparisons showed that the current model represented normal human morphology. The model was meshed with tetrahedral elements of 1 mm for the cartilage and labrum and 2–3 mm for the bone structures. A mesh size sensitivity analysis was performed to ensure stability of the model.

Spring 1D elements were used to model the ligaments based on descriptions of their origin and insertion areas. The resting lengths of the ligaments were defined through preliminary range of motion simulations and were found to be similar to measurements achieved by Yang et al. [48]. Tendons were modelled using triangular 2D elements. Muscles were modelled with 1D active spring elements with specific pulling forces depending on each muscle bundle (detailed hereafter). The tendons and ligaments were modelled using linear elastic mechanical properties. The bone structures were modelled using isotropic elastoplastic material laws specific to cortical and trabecular bone based on similar work by Astier et al. [4].

Finite element model function

The simulation approach consisted of imposing muscle forces and observing glenohumeral biomechanics throughout abduction. This approach enabled us to represent the history of bone loading and to identify key stages for all observed mechanical metrics. A linear ramp of increasing force was applied in each muscle element through active stimuli to simulate abduction movement. The force ratio between each individual muscle was based on muscle cross‐sectional areas [29]. A factor was used on the forces applied in the rotator cuff muscles to reach 60 N in the rotator cuff when the middle deltoid force reached 150 N [6]. To simulate abduction, only the forces in the middle deltoid, anterior deltoid, subscapularis, supraspinatus, infraspinatus and teres minor were activated [3, 6, 19, 20, 25, 26, 33, 40, 47]. The muscle activation was increased until a scapulohumeral angle of 60° was reached, or the maximum middle deltoid force of 150 N was reached, based on muscle strength limits identified by Baumgartner et al. [6].

The distal half of the scapula and the clavicle were considered a rigid body, and their movements were fixed in all directions. The distal two‐thirds of the humerus was also considered as a rigid body. A vertical force of 25 N was applied to the distal humerus to account for upper limb weight [13]. Nonpenetration contact interfaces were defined between the humerus and the scapula, as well as between soft tissues and bones, using a point/surface penalty method with a Coulomb‐type friction coefficient of 0.2 and a minimal gap of 0.05 mm [7]. The thickness of each tendon and muscle was taken into account in these contact interfaces to provide an accurate representation of muscle moment arms [8]. The same model geometry was used to simulate two patient configurations: a healthy patient with an intact rotator cuff and a rotator cuff‐deficient patient. Deficiency was represented by deactivating forces in the supraspinatus muscle.

Configurations and methods of measurement

The instrumented analysis was performed with the FX V135 (FX Shoulder Solutions) implants for two stem sizes (T14 and T16) in the following configurations (Figure 2): (1) anatomic TSA in the healthy patient and (2) RSA with symmetrical polyethylene (i.e., cylindrical with parallel faces, associated with a 135° bone cut) and asymmetrical polyethylene (i.e., cylindrical with angled faces, associated with a 145° bone cut) in the rotator cuff‐deficient patient.

Figure 2 Design of experiment and description of implant placement for the total shoulder arthroplasty (TSA) configuration (top row) and reverse shoulder arthroplasty (RSA) configurations (middle and bottom row) and their corresponding glenoid counterpart.

Virtual surgical process included removal of ligaments, cartilage, labrum and bone resection consistent with TSA and RSA surgeries. RSA also included removal of the superior cuff. Positioning of the stem and baseplate was blindly verified by two independent senior surgeons (Uma Srikumaran and Geoffroy Nourissat) to ensure the best positioning as requested by the manufacturer. A nonpenetration contact interface was applied between the humeral and glenoid implant components, as well as between bone and implant. A total of eight simulations were performed using an explicit dynamic Finite Element Model Solver (Radioss, release 2021.1, Altair Engineering, Inc.).

For each configuration, cortical and trabecular bone von Mises stress fields in the humerus and scapula were observed and compared. Quantitative measurements were joint kinematics (scapulohumeral angle), maximal von Mises stress for bone and implant components, contact forces between individual components of the model, and glenohumeral contact force resultant vector. Finally, the volume of humeral and scapular trabecular bone where stress exceeded 0.4 MPa was also measured to provide insight into bone remodelling and stress shielding areas, which is induced by volumetric mechanical loading of trabecular bone or lack thereof, respectively. The value of 0.4 MPa was chosen to facilitate comparison between configurations.

RESULTS

Glenohumeral joint kinematics and kinetics

The abduction kinematic analysis of the healthy and cuff‐deficient models revealed that a rupture of the supraspinatus produces a functionally limited shoulder, as expected, with only 17° of scapulohumeral angle reached at maximum muscle force (150 N in the middle deltoid). The scapulohumeral‐angle‐over‐time curves (Figure 3) of the uninstrumented models show an initial stabilization portion from 0 to 50 ms, as the coaptor muscles and weight of the arm were not given any initial state before 0 ms. Scapulohumeral angles increased with muscle activation. The healthy model reached a scapulohumeral angle of 60° (equivalent to 90° of abduction) with a total muscle force of 374 N.

Figure 3 Kinematic scapulohumeral angle and force activation curve on healthy and pseudo‐paralytic model (black) and instrumented models (red and blue) and evolution of abduction motion for configuration 1—total shoulder arthroplasty (TSA)—T14 (bottom row).

The placement of both TSA and RSA (with a deficient rotator cuff) enabled full restoration of function and range of motion similar to that of a healthy shoulder. RSA in the rotator cuff‐deficient shoulder restored function regardless of stem size and polyethylene shape. After RSA, less muscle activation was required to reach 60° of scapulohumeral angle compared with TSA configurations. The symmetric polyethylene reached a 60° scapulohumeral angle with slightly lower muscle forces compared with the asymmetric polyethylene. TSA produced a glenohumeral kinematic profile similar to that of the healthy model. For each implant configuration, stem size did not affect scapulohumeral angle kinematics (Table 1).

Table 1 Summary of measured mechanical metrics for all simulated configurations obtained at final humeral elevation (60° glenohumeral angle).

Configuration	Head‐Glenoid contact force (N)	Volume of stressed trabecular bone (>0.4 MPa) (cm3)	Von Mises stress (MPa)	
Humerus	Scapula	Humerus	Scapula	Implant HDPE (male)	Implant CrCo (female)	Implant other (TA6V)	
Cortical bone	Trabecular bone	Cortical bone	Trabecular bone	
1—Anat T14	311.17	20.72	3.68	37.32	7.89	18.53	3.17	7.81	163.25	228.08	
2—Anat T16	287.69	11.52	4.68	18.32	7.35	19.48	3.79	7.60	92.37	410.81	
3—Rev135 T14	219.22	1.50	0.42	5.47	1.61	15.75	5.24	3.11	25.69	340.12	
4—Rev135 T16	215.98	1.65	0.50	5.38	2.24	16.16	5.16	2.81	25.63	479.98	
5—Rev145 T14	226.15	2.78	1.28	5.99	2.68	17.03	6.23	2.58	22.77	581.38	
6—Rev145 T16	218.64	3.30	0.49	6.05	3.02	16.91	4.29	2.50	20.64	308.49	
Note: Head‐glenoid contact force relates to the total contact force between humeral and scapular implants. The implant stress data was separated by material, with HDPE corresponding to the polyethylene component, CrCo to the cobalt‐chrome component, and TA6V to all other components which are made of titanium.

John Wiley & Sons, Ltd.

More force was needed to reach a 60° scapulohumeral angle in the TSA model compared with the RSA model, resulting in more contact forces between the stem and the humeral trabecular bone (Figure 4a). The direction and contact point patterns were different for TSA and RSA configurations (Figure 4b–d), leading to different load distributions in the stem and adjacent bone. In the TSA configurations, the glenohumeral contact force resultant vectors were directed towards the centre of the shaft during all abduction motion. At maximum contact forces (60° scapulohumeral angle), the forces were directed towards the longitudinal axis direction of the stem. The contact forces between the trabecular bone and stem were evenly distributed around the stem shaft (Figure 4c). In RSA configurations, the glenohumeral contact force resultant vectors on the polyethylene were directed towards the centre of the stem in the early stages of abduction motion and then moved towards the medial part of the humerus, towards the calcar, as the abduction angle increased, leading to eccentric loading at maximal contact forces (60° scapulohumeral angle). The contact forces between the trabecular bone and stem were mainly located on smaller areas around the stem crown and medial to the distal tip of the stem (Figure 4e).

Figure 4 Head‐glenoid contact forces change between total shoulder arthroplasty (TSA) and reverse shoulder arthroplasty (RSA) but not in different configuration of RSA (a), glenohumeral reaction force directions at 10° and 60° abduction angle for TSA and RSA configurations (b, d), and Von Mises stress in the humeral trabecular bone superior to 0.1 MPa at 60° scapulohumeral angle (c, e).

Humerus stress analysis

In the TSA configurations, the higher stem size led to lower maximal local stress in the cortical and trabecular humeral bone. The stressed trabecular bone volume was also lower because there was less trabecular bone in the canal, and the forces were distributed over a larger surface area. Smaller stems led to higher stress in the proximal epiphysis close to the stem holes and metaphysis (Figure 5a,c,e).

Figure 5 Maximal local bone stress (chart) in cortical (a, b) and trabecular (c, d) bone of scapula and humerus, and volume of trabecular bone with stress superior to 0.4 MPa (e, f).

In the RSA configurations, increasing stem size did not affect cortical bone stress but did increase the trabecular bone maximal local stress and volume of stressed bone (>0.4 MPa) of the humerus, specifically around the medial part of the stem crown. Larger stems led to higher stress in the medial proximal epiphysis. The asymmetric polyethylene led to higher trabecular bone maximal local stress and volume of stressed bone (>0.4 MPa) (Table 1).

The lowest levels of trabecular bone stress were observed around the distal half of the humeral stem. Trabecular bone stress was higher for instrumented configurations compared with preoperative models.

Scapula stress analysis

The stress analysis on the scapula cortical bone showed slightly higher stress in TSA configurations compared with RSA because of higher muscle force activation (Figure 5a–f). There was no clear trend concerning the influence of stem size on cortical bone. In TSA configurations, the trabecular bone was less locally stressed and had a larger stressed bone volume because the polyethylene glenoid had a larger contact surface and was less rigid compared with the baseplate and screw system of the RSA configurations (Figure 5b,d,f). In RSA configurations, the high‐stress areas were located around the superior and posterior screw insertions. Polyethylene shape did not appear to influence scapular trabecular bone stress.

Implant stress

Stress analysis in the implant showed that the polyethylene glenoid (in TSA) was up to three times as stressed as the polyethylene (in RSA) because the muscle activation forces were higher, and the contact area was more localized. The larger stem size (T16) showed slightly lower high local stress compared with the T14 stem. The high‐stressed areas were located at the contact with the opposed sphere/glenosphere (Figure 6a). The maximal stress level was below the yield limit of the material.

Figure 6 Maximal local von Mises stress in implant polyethylene (a), titanium (b) and CoCr (c) components.

The stress analysis in the other metallic parts of the implants (titanium alloy and cobalt‐chromium alloy) showed that the high‐stressed areas were located at the connections (stem and taper or baseplate and screw connections). The maximal stress level was below the yield limit of the material (Figure 6b,c). There was no clear trend on the effect of the stem size and the metallic part maximal local stress.

DISCUSSION

Our study demonstrates that biomechanical stress around the humeral stem differs between RSA and TSA conditions. The first important finding is that the main difference in biomechanics between TSA and RSA is caused by the glenohumeral contact force vector, which we quantified and described for the specific configuration studied here (Figure 4b,d). TSA resulted in more loading around the stem shaft, and RSA resulted in more loading around the calcar. This is an important finding suggesting that clinical survival analysis of humeral short stems must be conducted in anatomical and reverse consideration. It should be not recommended to extrapolate radiological stress shielding occurring in TSA configuration identical to RSA consideration.

The second important finding is that smaller stems seem to be more appropriate for TSA, in which stress distribution around the implant allowed a floating implant, whereas larger stems generated higher local stress concentration in the proximal lateral humeral head. Larger stems seem to be more appropriate for RSA, because most trabecular bone stress was located medially to the implant, and greater canal fill allowed a better distribution of stress throughout the bone.

An asymmetric polyethylene component induced higher stresses in the metal structures of the implant compared with a symmetric polyethylene component. Asymmetrical polyethylene resulted in a greater abduction range of motion than symmetric polyethylene components. Higher polyethylene stress in RSA compared with TSA was expected because of low congruence between the radii of glenoid and humeral components and can be a predictor of earlier material wear.

In this finite element model, a superior rotator cuff deficiency demonstrated that RSA restores function when the supraspinatus is torn and leads to lower glenohumeral reaction forces. Ackland et al. [1] used musculoskeletal models of the glenohumeral joint to evaluate the effect of rotator cuff tear severity on joint and muscle forces after RSA and to quantify the stresses at the glenoid and humeral components. They simulated abduction and flexion after RSA in different cases of intact cuff and different types of tears of the cuff. The largest calculated glenohumeral joint contact forces, muscle forces and implant stresses were higher with increasing rotator cuff deficiency, which is consistent with our results.

In the case of supraspinatus deficiency in TSA, eccentric loading and an increased humeral head translation may be observed. It has been suggested that a downward inclination of the glenoid component may balance the supraspinatus deficiency, which would be interesting to evaluate in a future study.

Although we were able to analyse the effect of stem size on bone stress distribution, the two stem sizes tested here were similar. Further evaluation of a wider range of canal fill could provide additional insight into patient‐specific osseointegration; based on bone density or humeral metaphysis individual anatomy. However, the current results are consistent with previous literature. Langohr et al. [31] performed finite element analysis to quantify the effect of varying the size of the short stem humeral components having two different diametral sizes on the changes in bone stress from the preoperative intact to the reconstructed state for loading states consistent with 45° and 75° of abduction. They found that the smaller short‐stem implant produced humeral trabecular and cortical bone stresses that were closer to the intact state. The effect of humeral implant size and canal fill on bone stresses and interface contact in the proximal humerus was also evaluated by Synnott et al. [44], who compared three generic short stems with different cross‐sectional thickness (thinner, medium and thicker) to determine the effect of stem thickness on bone contact, bone stresses and bone resorption caused by stress shielding. They found that increasing the size had no significant effects on bone‐to‐implant contact during loading, but the thinner implant with the lowest canal fill ratio produced significantly less change in stress from the intact state for both cortical and trabecular bone. The thinner implant also resulted in less stress shielding and bone resorption compared with the two other designs. Although these studies align with our results, our conclusions may not apply to all implant designs.

In a finite element analysis by Barth et al. [5], the authors analysed the influence of proximal humeral stem geometry on stress distribution and torsional stability after TSA. We found that the stem shape of the FX V135 has a large proximal base/crown, and the only shape difference is in the distal part of the stem. Thus, the stability of the stem might not change much between stem sizes.

The credibility of the finite element model used in this study was assessed by comparing several simulation outputs with phenomena and measurements reported in the literature. The glenohumeral joint forces measured were consistent with the literature [11, 37, 46], although the ranges described vary widely depending on the study or the patient. Local kinematics of the healthy preoperative model showed a superoinferior translation of the contact area of the humeral head on the glenoid, which is also consistent with the literature [12, 22, 42]. Finally, while bone stress distribution is specific to each configuration studied, we verified that our results were in the same order of magnitude compared with previous studies [12, 22]. Also, bone and implant stresses were below yield stress, which is expected because the motion is physiologic.

Despite these elements of verification and validation, our study has some limitations. Although our model enabled 3D mobility of the humerus, the configurations we studied were tested only in abduction, and the simulation of other physiological movements should be explored to provide a better understanding of implant osseointegration. Also, bone stress is not a direct indicator of bone grown or osteolysis, and further analysis should be performed to link stress distribution and levels to osseointegration, as proposed by Fernandes et al. [21].

CONCLUSION

Our finite element study enabled quantification of differences in biomechanics between TSA and RSA. Based on the bone stress distribution analysis, we conclude that smaller stems may be more appropriate for TSA, while larger stems may be more appropriate for RSA. The current study demonstrated that different polyethylene shape on the same short stem does not increase humeral stresses but results in different ranges of motion and stresses in the metal structures of the implant.

AUTHOR CONTRIBUTIONS

All authors developed the methodology, wrote the original draft and revised the manuscript. Geoffroy Nourissat and Uma Srikumaran supervised the study, validated the results and revised the manuscript. Léo Fradet, Rohan‐Jean Bianco and Jean‐Marie Daudet reconceptualized the study, developed the methodology, managed the project and provided resources. Léo Fradet and Rohan‐Jean Bianco performed statistical analyses, developed the methodology, programmed the software, supervised the project and wrote the original draft. Geoffroy Nourissat and Uma Srikumaran conceptualized the study, developed the methodology and revised the manuscript.

CONFLICTS OF INTEREST STATEMENT

Geoffroy Nourissat, MD, PhD: Receiving reimbursements, fees, funding, or salary from an organization that may in any way gain or lose financially from the publication of the manuscript, either now or in the future. Victor Housset: Receiving reimbursements, fees, funding, or salary from an organization that may in any way gain or lose financially from the publication of the manuscript, either now or in the future. Jean‐Marie Daudet: Holding stocks or shares in an organization that may in any way gain or lose financially from the publication of the manuscript, either now or in the future. Léo Fradet, PhD, Eng.: Holding stocks or shares in an organization that may in any way gain or lose financially from the publication of the manuscript, either now or in the future. Rohan‐Jean Bianco, PhD.: Holding stocks or shares in an organization that may in any way gain or lose financially from the publication of the manuscript, either now or in the future. Uma Srikumaran, MD: Receiving reimbursements, fees, funding, or salary from an organization that may in any way gain or lose financially from the publication of the manuscript, either now or in the future.

ETHICS STATEMENT

Not applicable.

ACKNOWLEDGEMENTS

For editorial assistance, we thank Rachel Walden, MS, in the Editorial Services group of The Johns Hopkins Department of Orthopaedic Surgery. This study was funded by Fx Shoulder. The current work was founded by Philomec Inc and Fx Shoulder Solutions.

DATA AVAILABILITY STATEMENT

Due to commercial restrictions, supporting data are not available. The data sets used and/or analysed during the current study are available from the corresponding author upon reasonable request.
==== Refs
REFERENCES

1 Ackland, D.C. , Robinson, D.L. , Wilkosz, A. , Wu, W. , Richardson, M. , Lee, P. et al. (2019) The influence of rotator cuff tears on muscle and joint‐contact loading after reverse total shoulder arthroplasty. Journal of Orthopaedic Research, 37 (1 ), 211–219. Available from: 10.1002/jor.24152 30325061
2 Adams, C.R. , Baldwin, M.A. , Laz, P.J. , Rullkoetter, P.J. & Langenderfer, J.E. (2007) Effects of rotator cuff tears on muscle moment arms: a computational study. Journal of Biomechanics, 40 (15 ), 3373–3380. Available from: 10.1016/j.jbiomech.2007.05.017 17597135
3 Antony, N.T. & Keir, P.J. (2010) Effects of posture, movement and hand load on shoulder muscle activity. Journal of Electromyography and Kinesiology, 20 (2 ), 191–198. Available from: 10.1016/j.jelekin.2009.04.010 19473855
4 Astier, V. , Thollon, L. , Arnoux, P.J. , Mouret, F. & Brunet, C. (2008) Development of a finite element model of the shoulder: application during a side impact. International Journal of Crashworthiness, 13 (3 ), 301–312. Available from: 10.1080/13588260801933741
5 Barth, J. , Garret, J. , Geais, L. , Bothorel, H. , Saffarini, M. , Shoulder Friends Institute , Godenèche, A. (2019) Influence of uncemented humeral stem proximal geometry on stress distributions and torsional stability following total shoulder arthroplasty. Journal of Experimental Orthopaedics, 6 (1 ), 8. Available from: 10.1186/s40634-019-0178-4 30805747
6 Baumgartner, D. , Tomas, D. , Gossweiler, L. , Siegl, W. , Osterhoff, G. & Heinlein, B. (2014) Towards the development of a novel experimental shoulder simulator with rotating scapula and individually controlled muscle forces simulating the rotator cuff. Medical & Biological Engineering & Computing, 52 (3 ), 293–299. Available from: 10.1007/s11517-013-1120-z 24170552
7 Bianco, R.‐J. , Arnoux, P.‐J. , Wagnac, E. , Mac‐Thiong, J.‐M. & Aubin, C.‐É. (2017) Minimizing pedicle screw pullout risks: a detailed biomechanical analysis of screw design and placement. Clinical Spine Surgery: A Spine Publication, 30 (3 ), E226–E232. Available from: 10.1097/BSD.0000000000000151
8 Bouaicha, S. , Ernstbrunner, L. , Jud, L. , Meyer, D.C. , Snedeker, J.G. & Bachmann, E. (2018) The lever arm ratio of the rotator cuff to deltoid muscle explains and predicts pseudoparalysis of the shoulder: the Shoulder Abduction Moment index. The Bone & Joint Journal, 100–B (12 ), 1600–1608. Available from: 10.1302/0301-620X.100B12.BJJ-2018-0493.R1
9 Burke, C.S. , Roberts, C.S. , Nyland, J.A. , Radmacher, P.G. , Acland, R.D. & Voor, M.J. (2006) Scapular thickness—implications for fracture fixation. Journal of Shoulder and Elbow Surgery, 15 (5 ), 645–648. Available from: 10.1016/j.jse.2005.10.005 16979064
10 Cambra‐Moo, O. , Nacarino Meneses, C. , Rodríguez Barbero, M.Á. , García Gil, O. , Rascón Pérez, J. , Rello‐Varona, S. et al. (2014) An approach to the histomorphological and histochemical variations of the humerus cortical bone through human ontogeny. Journal of Anatomy, 224 (6 ), 634–646. Available from: 10.1111/joa.12172 24660964
11 Charlton, I.W. & Johnson, G.R. (2006) A model for the prediction of the forces at the glenohumeral joint. Proceedings of the Institution of Mechanical Engineers, Part H: Journal of Engineering in Medicine, 220 (8 ), 801–812. Available from: 10.1243/09544119JEIM147 17236514
12 Comenda, M. , Quental, C. , Folgado, J. , Sarmento, M. & Monteiro, J. (2019) Bone adaptation impact of stemless shoulder implants: a computational analysis. Journal of Shoulder and Elbow Surgery, 28 (10 ), 1886–1896. Available from: 10.1016/j.jse.2019.03.007 31255444
13 Debski, R.E. , McMahon, P.J. , Thompson, W.O. , Woo, S.L.Y. , Warner, J.J.P. & Fu, F.H. (1995) A new dynamic testing apparatus to study glenohumeral joint motion. Journal of Biomechanics, 28 (7 ), 869–874. Available from: 10.1016/0021-9290(95)95276-B 7657685
14 Dekker, T.J. , Aman, Z.S. , Peebles, L.A. , Storaci, H.W. , Chahla, J. , Millett, P.J. et al. (2020) Quantitative and qualitative analyses of the glenohumeral ligaments: an anatomic study. The American Journal of Sports Medicine, 48 (8 ), 1837–1845. Available from: 10.1177/0363546520917665 32441993
15 Denard, P.J. , Haidamous, G. , Gobezie, R. , Romeo, A.A. & Lederman, E. (2020) Short‐term evaluation of humeral stress shielding following reverse shoulder arthroplasty using press‐fit fixation compared with cemented fixation. Journal of Shoulder and Elbow Surgery, 29 (5 ), 906–912. Available from: 10.1016/j.jse.2019.09.042 31911215
16 Denard, P.J. , Noyes, M.P. , Walker, J.B. , Shishani, Y. , Gobezie, R. , Romeo, A.A. et al. (2018) Proximal stress shielding is decreased with a short stem compared with a traditional‐length stem in total shoulder arthroplasty. Journal of Shoulder and Elbow Surgery, 27 (1 ), 53–58. Available from: 10.1016/j.jse.2017.06.042 28865965
17 Denard, P.J. , Raiss, P. , Gobezie, R. , Edwards, T.B. & Lederman, E. (2018) Stress shielding of the humerus in press‐fit anatomic shoulder arthroplasty: review and recommendations for evaluation. Journal of Shoulder and Elbow Surgery, 27 (6 ), 1139–1147. Available from: 10.1016/j.jse.2017.12.020 29422391
18 Elwell, J.A. , Athwal, G.S. & Willing, R. (2018) Development and validation of a muscle wrapping model applied to intact and reverse total shoulder arthroplasty shoulders. Journal of Orthopaedic Research, 36 (12 ), 3308–3317. Available from: 10.1002/jor.24131 30144134
19 Favre, P. , Senteler, M. , Hipp, J. , Scherrer, S. , Gerber, C. & Snedeker, J.G. (2012) An integrated model of active glenohumeral stability. Journal of Biomechanics, 45 (13 ), 2248–2255. Available from: 10.1016/j.jbiomech.2012.06.010 22818663
20 Favre, P. , Sheikh, R. , Fucentese, S.F. & Jacob, H.A.C. (2005) An algorithm for estimation of shoulder muscle forces for clinical use. Clinical Biomechanics, 20 (8 ), 822–833. Available from: 10.1016/j.clinbiomech.2005.04.007 15975696
21 Fernandes, P. , Rodrigues, H. & Jacobs, C. (1999) A model of bone adaptation using a global optimisation criterion based on the trajectorial theory of Wolff. Computer Methods in Biomechanics and Biomedical Engineering, 2 (2 ), 125–138. Available from: 10.1080/10255849908907982 11264822
22 Filardi, V. (2020) Stress distribution in the humerus during elevation of the arm and external abduction. Journal of Orthopaedics, 19 , 218–222. Available from: 10.1016/j.jor.2020.02.003 32071517
23 Fox, J.A. , Cole, B.J. , Romeo, A.A. , Meininger, A.K. , Williams, J.M. , Glenn, R.E. et al. (2008) Articular cartilage thickness of the humeral head: an anatomic study. Orthopedics, 31 (3 ), 216. Available from: 10.3928/01477447-20080301-11 19292258
24 Giannicola, G. , Spinello, P. , Scacchi, M. & Gumina, S. (2017) Cartilage thickness of distal humerus and its relationships with bone dimensions: magnetic resonance imaging bilateral study in healthy elbows. Journal of Shoulder and Elbow Surgery, 26 (5 ), e128–e136. Available from: 10.1016/j.jse.2016.10.012 28131685
25 Halder, A.M. , Itoi, E. & An, K.N. (2000) Anatomy and biomechanics of the shoulder. Orthopedic Clinics of North America, 31 (2 ), 159–176. Available from: 10.1016/S0030-5898(05)70138-3 10736387
26 Hughes, R.E. & An, K.N. (1996) Force analysis of rotator cuff muscles. Clinical Orthopaedics and Related Research, 330 , 75–83. 10.1097/00003086-199609000-00010
27 Jacobson, A. , Gilot, G.J. , Hamilton, M.A. , Greene, A. , Flurin, P.‐H. , Wright, T.W. et al. (2015) Glenohumeral anatomic study. a comparison of male and female shoulders with similar average age and BMI. Bulletin of the Hospital for Joint Diseases, 73 (Suppl 1 ), S68–S78.26631200
28 Kadavkolan, A.S. & Jawhar, A. (2018) Glenohumeral joint morphometry with reference to anatomic shoulder arthroplasty. Current Orthopaedic Practice, 29 (1 ), 71–83. Available from: 10.1097/BCO.0000000000000552
29 Karlsson, D. & Peterson, B. (1992) Towards a model for force predictions in the human shoulder. Journal of Biomechanics, 25 (2 ), 189–199. Available from: 10.1016/0021-9290(92)90275-6 1733994
30 Kim, S.C. , Park, J.H. , Bukhary, H. & Yoo, J.C. (2022) Humeral stem with low filling ratio reduces stress shielding in primary reverse shoulder arthroplasty. International Orthopaedics, 46 (6 ), 1341–1349. Available from: 10.1007/s00264-022-05383-4 35353240
31 Langohr, G.D.G. , Reeves, J. , Roche, C.P. , Faber, K.J. & Johnson, J.A. (2020) The effect of short‐stem humeral component sizing on humeral bone stress. Journal of Shoulder and Elbow Surgery, 29 (4 ), 761–767. Available from: 10.1016/j.jse.2019.08.018 31711829
32 Lewis, G.S. & Armstrong, A.D. (2011) Glenoid spherical orientation and version. Journal of Shoulder and Elbow Surgery, 20 (1 ), 3–11. Available from: 10.1016/j.jse.2010.05.012 20932782
33 Longosz, A.F. , Serwadda, D. , Nalugoda, F. , Kigozi, G. , Franco, V. , Gray, R.H. et al. (2014) Impact of HIV subtype on performance of the limiting antigen‐avidity enzyme immunoassay, the Bio‐Rad avidity assay, and the BED capture immunoassay in Rakai, Uganda. AIDS Research and Human Retroviruses, 30 (4 ), 339–344. Available from: 10.1089/aid.2013.0169 24083837
34 Majed, A. , Thangarajah, T. , Southgate, D. , Reilly, P. , Bull, A. & Emery, R. (2019) Cortical thickness analysis of the proximal humerus. Shoulder & Elbow, 11 (2 ), 87–93. Available from: 10.1177/1758573217736744 30936946
35 Matsumura, N. , Ogawa, K. , Ikegami, H. , Collin, P. , Walch, G. & Toyama, Y. (2014) Computed tomography measurement of glenoid vault version as an alternative measuring method for glenoid version. Journal of Orthopaedic Surgery and Research, 9 (1 ), 17. Available from: 10.1186/1749-799X-9-17 24618285
36 Somesh, M.S. , Prabhu, L.V. , Shilpa, K. , Pai, M.M. , Krishnamurthy, A. & Murlimanju, B.V. (2011) Morphometric study of the humerus segments in Indian population. International Journal of Morphology, 29 , 1174–1180. Available from: 10.4067/S0717-95022011000400017
37 Nikooyan, A.A. , Veeger, H.E.J. , Westerhoff, P. , Bolsterlee, B. , Graichen, F. , Bergmann, G. et al. (2012) An EMG‐driven musculoskeletal model of the shoulder. Human Movement Science, 31 (2 ), 429–447. Available from: 10.1016/j.humov.2011.08.006 22244106
38 Pitocchi, J. , Wirix‐Speetjens, R. , van Lenthe, G.H. & Pérez, M.Á. (2020) Integration of cortical thickness data in a statistical shape model of the scapula. Computer Methods in Biomechanics and Biomedical Engineering, 23 (10 ), 642–648. Available from: 10.1080/10255842.2020.1757082 32364819
39 Polguj, M. , Jędrzejewski, K.S. , Podgórski, M. & Topol, M. (2011) Correlation between morphometry of the suprascapular notch and anthropometric measurements of the scapula. Folia Morphologica, 70 (2 ), 109–115.21630232
40 Quental, C. , Folgado, J. , Ambrósio, J. & Monteiro, J. (2013) Multibody system of the upper limb including a reverse shoulder prosthesis. Journal of Biomechanical Engineering, 135 (11 ), 111005. Available from: 10.1115/1.4025325.24008920
41 Ruckstuhl, H. , Krzycki, J. , Petrou, N. , Vanwanseele, B. & Stussi, E. (2008) A quantitative study of humeral cartilage in individuals with spinal cord injury. Spinal Cord, 46 (2 ), 129–134. Available from: 10.1038/sj.sc.3102084 17563754
42 Soslowsky, L.J. , Flatow, E.L. , Bigliani, L.U. & Mow, V.C. (1992) Articular geometry of the glenohumeral joint. Clinical Orthopaedics and Related Research, (285 ), 181–190.
43 Sprecher, C.M. , Schmidutz, F. , Helfen, T. , Richards, R.G. , Blauth, M. & Milz, S. (2015) Histomorphometric assessment of cancellous and cortical bone material distribution in the proximal humerus of normal and osteoporotic individuals. Medicine, 94 (51 ), e2043. Available from: 10.1097/MD.0000000000002043.26705200
44 Synnott, S. , Langohr, G.D.G. , Reeves, J.M. , Johnson, J.A. & Athwal, G.S. (2021) The effect of humeral implant thickness and canal fill on interface contact and bone stresses in the proximal humerus. JSES International, 5 (5 ), 881–888. Available from: 10.1016/j.jseint.2021.05.006 34505100
45 Webb, J.D. , Blemker, S.S. & Delp, S.L. (2014) 3D finite element models of shoulder muscles for computing lines of actions and moment arms. Computer Methods in Biomechanics and Biomedical Engineering, 17 (8 ), 829–837. Available from: 10.1080/10255842.2012.719605 22994141
46 Westerhoff, P. , Graichen, F. , Bender, A. , Halder, A. , Beier, A. , Rohlmann, A. et al. (2009) In vivo measurement of shoulder joint loads during activities of daily living. Journal of Biomechanics, 42 (12 ), 1840–1849. Available from: 10.1016/j.jbiomech.2009.05.035 19643418
47 Wuelker, N. , Wirth, C.J. , Plitz, W. & Roetman, B. (1995) A dynamic shoulder model: reliability testing and muscle force study. Journal of Biomechanics, 28 (5 ), 489–499. Available from: 10.1016/0021-9290(94)E0006-O 7775486
48 Yang, C. , Goto, A. , Sahara, W. , Yoshikawa, H. & Sugamoto, K. (2010) In vivo three‐dimensional evaluation of the functional length of glenohumeral ligaments. Clinical Biomechanics, 25 (2 ), 137–141. Available from: 10.1016/j.clinbiomech.2009.10.009 19942330
49 Zumstein, V. , Kraljević, M. , Hoechel, S. , Conzen, A. , Nowakowski, A. & Müller‐Gerbl, M. (2014) The glenohumeral joint—a mismatching system? A morphological analysis of the cartilaginous and osseous curvature of the humeral head and the glenoid cavity. Journal of Orthopaedic Surgery and Research, 9 , 34. Available from: 10.1186/1749-799X-9-34 24886613
