Deep learning models (DLM) are efficient replacements for computationally intensive optimization techniques. Musculoskeletal models (MSM) typically involve resource-intensive optimization processes for determining joint and muscle forces. Consequently, DLM could predict MSM results and reduce computational costs. Within the total shoulder arthroplasty (TSA) domain, the glenohumeral joint force represents a critical MSM outcome as it can influence joint function, joint stability, and implant durability. Here, we aimed to employ deep learning techniques to predict both the magnitude and direction of the glenohumeral joint force. To achieve this, 959 virtual subjects were generated using the Markov-Chain Monte-Carlo method, providing patient-specific parameters from an existing clinical registry. A DLM was constructed to predict the glenohumeral joint force components within the scapula coordinate system for the generated subjects with a coefficient of determination of 0.97, 0.98, and 0.98 for the three components of the glenohumeral joint force. The corresponding mean absolute errors were 11.1, 12.2, and 15.0 N, which were about 2% of the maximum glenohumeral joint force. In conclusion, DLM maintains a comparable level of reliability in glenohumeral joint force estimation with MSM, while drastically reducing the computational costs.