The majority of manual wheelchair users (MWCU) develop shoulder pain or injuries, which is often caused by impingement. Because propulsion mechanics are influenced by the recovery hand pattern used, the pattern may affect shoulder loading and susceptibility to injury. Shoulder muscle weakness is also correlated with shoulder pain, but how shoulder loading changes with specific muscle group weakness is unknown. Musculoskeletal modeling and simulation were used to compare glenohumeral joint contact forces (GJCFs) across hand patterns and determine how GJCFs vary when primary shoulder muscle groups are weakened. Experimental data were analyzed to classify individuals into four hand pattern groups. A representative musculoskeletal model was then developed for each group and simulations generated to portray baseline strength and six muscle weakness conditions. Three-dimensional GJCF peaks and impulses were compared across hand patterns and muscle weakness conditions. The semicircular pattern consistently had lower shear (anterior-posterior and superior-inferior) GJCFs compared to other patterns. The double-loop pattern had the highest superior GJCFs, while the single-loop pattern had the highest anterior and posterior GJCFs. These results suggest that using the semicircular pattern may be less susceptible to shoulder injuries such as subacromial impingement. Weakening the internal rotators and external rotators resulted in the greatest increases in shear GJCFs and decreases in compressive GJCF, likely due to decreased force from rotator cuff muscles. These findings suggest that strengthening specific muscle groups, especially the rotator cuff, is critical for decreasing the risk of shoulder overuse injuries.