In vivo bone strain and finite element modeling of a rhesus macaque mandible during mastication

• Published in: Zoology (Jena) Vol. 124; pp. 13 - 29
• Main Authors: Panagiotopoulou, Olga, Iriarte-Diaz, José, Wilshin, Simon, Dechow, Paul C., Taylor, Andrea B., Mehari Abraha, Hyab, Aljunid, Sharifah F., Ross, Callum F.
• Format: Journal Article
• Language: English
• Published: Germany Elsevier GmbH 01.10.2017
• Subjects: Animals; Biomechanical Phenomena; Bone and Bones - physiology; Bone material properties; Bone strain gauges; Chewing; Electromyography; Female; Finite element analysis; Macaca mulatta - physiology; Mandible - physiology; Mastication - physiology; Models, Biological; Musculoskeletal modeling; Bone material properties; Musculoskeletal modeling; Finite element analysis; Chewing; Bone strain gauges
• DOI: 10.1016/j.zool.2017.08.010
• ISSN: 0944-2006; 1873-2720

•First study to combine in vivo and ex vivo data to create and validate a finite element model simulation of animal feeding.•The similarities between the strain regimes recorded in vivo and those generated by the model are marked.•Variation in the bite point is associated with significant variation in strain regimes in the mandible.•Relative timing and magnitude of muscle forces are important determinants of this variation in strain regimes. Finite element analysis (FEA) is a commonly used tool in musculoskeletal biomechanics and vertebrate paleontology. The accuracy and precision of finite element models (FEMs) are reliant on accurate data on bone geometry, muscle forces, boundary conditions and tissue material properties. Simplified modeling assumptions, due to lack of in vivo experimental data on material properties and muscle activation patterns, may introduce analytical errors in analyses where quantitative accuracy is critical for obtaining rigorous results. A subject-specific FEM of a rhesus macaque mandible was constructed, loaded and validated using in vivo data from the same animal. In developing the model, we assessed the impact on model behavior of variation in (i) material properties of the mandibular trabecular bone tissue and teeth; (ii) constraints at the temporomandibular joint and bite point; and (iii) the timing of the muscle activity used to estimate the external forces acting on the model. The best match between the FEA simulation and the in vivo experimental data resulted from modeling the trabecular tissue with an isotropic and homogeneous Young’s modulus and Poisson’s value of 10GPa and 0.3, respectively; constraining translations along X,Y, Z axes in the chewing (left) side temporomandibular joint, the premolars and the m1; constraining the balancing (right) side temporomandibular joint in the anterior-posterior and superior-inferior axes, and using the muscle force estimated at time of maximum strain magnitude in the lower lateral gauge. The relative strain magnitudes in this model were similar to those recorded in vivo for all strain locations. More detailed analyses of mandibular strain patterns during the power stroke at different times in the chewing cycle are needed.