The finite element method described in this study provides an easy method to simulate the kinetics of multibody mechanisms. It is used in order to develop a musculoskeletal model of the shoulder mechanism. Each relevant morphological structure has been represented by an appropriate element. For the shoulder mechanism two special-purpose elements have been developed: a SURFACE element representing the scapulothoracic gliding plane and a CURVED-TRUSS element to represent muscles which are wrapped around bony contours. The model contains four bones, three joints, three extracapsular ligaments, the scapulothoracic gliding plane and 20 muscles and muscle parts. In the model, input variables are the positions of the shoulder girdle and humerus and the external load on the humerus. Output variables are muscle forces subject to an optimization procedure in which the mechanical stability of the glenohumeral joint is one of the constraints. Four different optimization criteria are compared. For 12 muscles, surface EMG is used to verify the model. Since the optimum muscle length and force-length relationship are unknown, and since maximal EMG amplitude is length dependent, verification is only possible in a qualitative sense. Nevertheless, it is concluded that a detailed model of the shoulder mechanism has been developed which provides good insight into the function of morphological structures.