Musculoskeletal models are intended to be used to assist in prevention and treatments of musculoskeletal disorders. To capture important aspects of shoulder dysfunction, realistic simulation of clavicular and scapular movements is crucial. The range of motion of these bones is dependent on thoracic, clavicular and scapular anatomy and therefore different for each individual. Typically, patient or subject measurements will therefore not fit on a model that uses a cadaveric morphology. Up till now, this problem was solved by adjusting measured bone rotations such that they fit on the model, but this leads to adjustments of on average 3.98° and, in some cases, even more than 8°. Two novel methods are presented that decrease this discrepancy between experimental data and simulations. For one method, the model is scaled to fit the subject, leading to a 34 % better fit compared to the existing method. In the other method, the set of possible joint rotations is increased by allowing some variation on motion constraints, resulting in a 42 % better fit. This change in kinematics also affected the kinetics: muscle forces of some important scapular stabilizing muscles, as predicted by the Delft Shoulder and Elbow Model, were altered by maximally 17 %. The effect on the glenohumeral joint contact force was however marginal (1.3 %). The methods presented in this paper might lead to more realistic shoulder simulations and can therefore be considered a step towards (clinical) application, especially for applications that involve scapular imbalance. © 2013 International Federation for Medical and Biological Engineering.