Dogs are an interesting object of investigation because of the wide range of body size, body mass, and physique. In the last several years, the number of clinical and biomechanical studies on dog locomotion has increased. However, the relationship between body structure and joint load during locomotion, as well as between joint load and degenerative diseases of the locomotor system (e.g. dysplasia), are not sufficiently understood. In vivo measurements/records of joint forces and loads or deep/small muscles are complex, invasive, and sometimes ethically questionable. The use of detailed musculoskeletal models may help in filling that knowledge gap. We describe here the methods we used to create a detailed musculoskeletal model with 84 degrees of freedom and 134 muscles. Our model has three key-features: Three-dimensionality, scalability, and modularity. We tested the validity of the model by identifying forelimb muscle synergies of a beagle at walk. We used inverse dynamics and static optimization to estimate muscle activations based on experimental data. We identified three muscle synergy groups by using hierarchical clustering. Predicted activation patterns exhibited good agreement with experimental data for most of the forelimb muscles. We expect that our model will speed up the analysis of how body size, physique, agility, and disease influence joint neuronal control and loading in dog locomotion.