The masticatory region is an important human motion system that is essential for basic human tasks like mastication, speech or swallowing. An association between temporomandibular disorders (TMDs) and high temporomandibular joint (TMJ) stress has been suggested, but in vivo joint force measurements are not feasible to directly test this assumption. Consequently, biomechanical computer simulation remains as one of a few means to investigate this complex system. To thoroughly examine orofacial biomechanics, we developed a novel, dynamic computer model of the masticatory system. The model combines a muscle driven rigid body model of the jaw region with a detailed finite element model (FEM) disk and elastic foundation (EF) articular cartilage. The model is validated using high-resolution MRI data for protrusion and opening that were collected from the same volunteer. Joint stresses for a clenching task as well as protrusive and opening movements are computed. Simulations resulted in mandibular positions as well as disk positions and shapes that agree well with the MRI data. The model computes reasonable disk stress patterns for dynamic tasks. Moreover, to the best of our knowledge this model presents the first ever contact model using a combination of EF layers and a FEM body, which results in a clear decrease in computation time. In conclusion, the presented model is a valuable tool for the investigation of the human TMJ and can potentially help in the future to increase the understanding of the masticatory system and the relationship between TMD and joint stress and to highlight potential therapeutic approaches for the restoration of orofacial function.

No abstract

Graphical abstract

Increased mechanical loading of the temporomandibular joint (TMJ) is often connected with the onset and progression of temporomandibular joint disorders (TMD). The potential role of occlusal factors and sleep bruxism in the onset of TMD are a highly debated topic in literature, but ethical considerations limit in vivo examinations of this problem. The study aims to use an innovative in silico modeling approach to thoroughly investigate the connection between morphological parameters, bruxing direction and TMJ stress. A forward-dynamics tracking approach was used to simulate laterotrusive and mediotrusive tooth grinding for 3 tooth positions, 5 lateral inclination angles, 5 sagittal tilt angles and 3 force levels, giving a total of 450 simulations. Muscle activation patterns, TMJ disc von Mises stress as well as correlations between mean muscle activations and TMJ disc stress are reported. Computed muscle activation patterns agree well with previous literature. The results suggest that tooth inclination and grinding position, to a smaller degree, have an effect on TMJ loading. Mediotrusive bruxing computed higher loads compared to laterotrusive simulations. The strongest correlation was found for TMJ stress and mean activation of the superficial masseter. Overall, our results provide in silico evidence that TMJ disc stress is related to tooth morphology.

Our formulation is not only valuable for shoulder simulations, but could be used in various clinical situations (e.g. for different joints and rehabilitation therapy tasks) where the direction and/or magnitude of reaction forces are of interest.