The Mechanical Effect of the Periodontal Ligament on Bone Strain Regimes in a Validated Finite Element Model of a Macaque Mandible

A new biomechanical model of the mammal jaw based on load path analysis

The primary function of the tetrapod jaw is to transmit jaw muscle forces to bite points. The routes of force transfer in the jaw have never been studied but can be quantified using load paths - the shortest, stiffest routes from regions of force application to support constraints. Here, we use load path analysis to map force transfer from muscle attachments to bite point and jaw joint, and to evaluate how different configurations of trabecular and cortical bone affect load paths. We created three models of the mandible of the Virginia opossum, Didelphis virginiana, each with a cortical bone shell, but with different material properties for the internal spaces: (1) a cortical-trabecular model, in which the interior space is modeled with bulk properties of trabecular bone; (2) a cortical-hollow model, in which trabeculae and mandibular canal are modeled as hollow; and (3) a solid-cortical model, in which the interior is modeled as cortical bone. The models were compared with published in vivo bite force and bone strain data, and the load paths calculated for each model. The trabecular model, which is preferred because it most closely approximates the actual morphology, was best validated by in vivo data. In all three models, the load path was confined to cortical bone, although its route within the cortex varied depending on the material properties of the inner model. Our analysis shows that most of the force is transferred through the cortical, rather than trabecular bone, and highlights the potential of load path analysis for understanding form-function relationships in the skeleton.

Comparative assessment of orthodontic clear aligner versus fixed appliance for anterior retraction: a finite element study

BACKGROUND

The aim of this study is to conduct a comparative evaluation of different designs of clear aligners and examine the disparities between clear aligners and fixed appliances.

METHODS

3D digital models were created, consisting of a maxillary dentition without first premolars, maxilla, periodontal ligaments, attachments, micro-implant, 3D printed lingual retractor, brackets, archwire and clear aligner. The study involved the creation of five design models for clear aligner maxillary anterior internal retraction and one design model for fixed appliance maxillary anterior internal retraction, which were subsequently subjected to finite element analysis. These design models included: (1) Model C0 Control, (2) Model C1 Posterior Micro-implant, (3) Model C2 Anterior Micro-implant, (4) Model C3 Palatal Plate, (5) Model C4 Lingual Retractor, and (6) Model F0 Fixed Appliance.

RESULTS

In the clear aligner models, a consistent pattern of tooth movement was observed. Notably, among all tested models, the modified clear aligner Model C3 exhibited the smallest differences in sagittal displacement of the crown-root of the central incisor, vertical displacement of the central incisor, sagittal displacement of the second premolar and second molar, as well as vertical displacement of posterior teeth. However, distinct variations in tooth movement trends were observed between the clear aligner models and the fixed appliance model. Furthermore, compared to the fixed appliance model, significant increases in tooth displacement were achieved with the use of clear aligner models.

CONCLUSIONS

In the clear aligner models, the movement trend of the teeth remained consistent, but there were variations in the amount of tooth displacement. Overall, the Model C3 exhibited better torque control and provided greater protection for posterior anchorage teeth compared to the other four clear aligner models. On the other hand, the fixed appliance model provides superior anterior torque control and better protection of the posterior anchorage teeth compared to clear aligner models. The clear aligner approach and the fixed appliance approach still exhibit a disparity; nevertheless, this study offers a developmental direction and establishes a theoretical foundation for future non-invasive, aesthetically pleasing, comfortable, and efficient modalities of clear aligner treatment.

[Compressive stress in three types of finishing lines with lithium disilicate crowns in permanent teeth: finite element analysis]

Introduction

In oral rehabilitation, the use of ceramic restorations is widely accepted due to its aesthetic capacity to mimic the naturalness of the dental tissue, provide longevity of the material, and present a greater marginal fit compared to crowns with a metal structure. Termination lines are biological preparations whose function is to minimize the cervical opening of the marginal seal. Consequently, analyzing the behavior of restorative materials under compressive forces decreases the risk of fracture and increases the success of the treatment.

Objective

To compare the compressive stresses of lithium disilicate crowns with three different finishing lines.

Methodology

In silico study of the simulation of a dental preparation on a lower right first molar with chamfer (0.6 mm), shoulder (0.5 mm) and deep chamfer (0.5 mm) finish lines. Using the SolidWorks®️ version 2017 software, the maximum stresses, minimum stresses, and location of the compressive force were collected on the Megapascal (Mpa) measurement scale.

Results

The chamfer type termination line (0.6mm) obtained a lower compressive stress compared to the other two shoulder type termination lines (0.5mm) and deep chamfer (0.5mm).

Conclusions

It was shown that the chamfer type finishing line (0.6mm) presented a better force distribution, determining greater reliability in the selection of this finishing line with the use of a lithium disilicate crown in a unitary manner.

Dynamic finite element modelling of the macaque mandible during a complete mastication gape cycle

Three-dimensional finite element models (FEMs) are powerful tools for studying the mechanical behaviour of the feeding system. Using validated, static FEMs we have previously shown that in rhesus macaques the largest food-related differences in strain magnitudes during unilateral postcanine chewing extend from the lingual symphysis to the endocondylar ridge of the balancing-side ramus. However, static FEMs only model a single time point during the gape cycle and probably do not fully capture the mechanical behaviour of the jaw during mastication. Bone strain patterns and moments applied to the mandible are known to vary during the gape cycle owing to variation in the activation peaks of the jaw-elevator muscles, suggesting that dynamic models are superior to static ones in studying feeding biomechanics. To test this hypothesis, we built dynamic FEMs of a complete gape cycle using muscle force data from in vivo experiments to elucidate the impact of relative timing of muscle force on mandible biomechanics. Results show that loading and strain regimes vary across the chewing cycle in subtly different ways for different foods, something which was not apparent in static FEMs. These results indicate that dynamic three-dimensional FEMs are more informative than static three-dimensional FEMs in capturing the mechanical behaviour of the jaw during feeding by reflecting the asymmetry in jaw-adductor muscle activations during a gape cycle. This article is part of the theme issue 'Food processing and nutritional assimilation in animals'.

Biomechanical influence of plate configurations on mandible subcondylar fracture fixation: a finite element study

Mandible subcondylar fractures have very high complication rate, yet there is no consensus on suitable plate design for optimal patient outcomes. Our study is aimed at comparing single mini, trapezoid, lambda, strut and double mini plates. A finite-element (FE) model of intact mandible was developed based on healthy CT-scan data, which was further virtually osteotomized and fixated with plates. The cortical and cancellous bones were assigned region-specific orthotropic and heterogenous isotropic material properties respectively. The models were subjected to six load cases representing the mastication cycle. Under opposite lateralities, the tensile and compressive mandibular strain distributions were found as the opposite, with tensile strains at the posterior border under ipsilateral molar clenching (RMOL) resulting in lesser mandibular strain in reconstructed mandible with single mini plate under RMOL but highest mandibular strain under the contralateral molar clenching (LMOL). Owing to the reduced mandibular strains under LMOL than RMOL, the contralateral chewing is preferred during the immediate post-surgery period for patients. Under LMOL, the peak von Mises stresses in the plate decreased with increase in the number of screws. Furthermore, the presence of two arms in double mini and trapezoid plates seems beneficial to neutralise the tensile and compressive strains across load cases.

Crack Propagation and Fatigue Performance of Partial Posterior Indirect Restorations: An Extended Finite Element Method Study

Dental ceramics are susceptible to slow, progressive crack growth after cyclic loading. The purpose of this study was to investigate the progressive patterns of cracks in two different types of CAD/CAM ceramic materials used with three different partial posterior indirect restoration (PPIR) designs and to determine the materials' failure risk using a fatigue test. Standard initial cracks were formed in Standard Tessellation Language (STL) files prepared for three different PPIRs. The materials chosen were monolithic lithium disilicate (LS) and polymer-infiltrated ceramic networks (PICNs). The extended finite element method (XFEM) was applied, and the fatigue performance was examined by applying a 600 N axial load. The cracks propagated the most in onlay restorations, where the highest displacement was observed. In contrast, the most successful results were observed in overlay restorations. Overlay restorations also showed better fatigue performance. LS materials exhibited more successful results than PICN materials. LS materials, which can be used in PPIRs, yield better results compared to PICN materials. While inlay restorations demonstrated relatively successful results, overlay and onlay restorations can be specified as the most and the least successful PPIR types, respectively.

Choosing the optimal mandible position for inferior alveolar nerve block (IANB) using finite element analysis

BACKGROUND

One of the most popular methods of local anesthesia in dentistry, inferior alveolar nerve block (IANB) involves the blockade of the inferior alveolar nerve (IAN) and lingual nerve (LN) in the pterygomandibular space. Despite the large number of works describing the contents of this space, the spatial displacements of the anatomical structures of this area at different positions of the mandible have not been sufficiently studied. The aim of our study was to study the spatial movements of the IAN and inferior alveolar artery (IAA) at various positions of the mandible using computer simulation and finite element analysis to find the safest way to conduct IANB.

MATERIALS AND METHODS

Reverse engineering was used to create a model of the cranial base and the mandible based on computed tomography (CT) data obtained from patient N (male, 24 years old), the arteries of the head and neck were designed from the data of multiphase angiography of patient M (female, 61 years old). Masticatory muscles, sphenomandibular ligament, temporomandibular joint and mandibular nerve were modeled in the SolidWorks software package based on an open database of anatomical structures. The finite element grid was generated in the Solidworks software. In the first series of experiments, the displacement of the mandible was modeled along the vertical axis down by 48 mm (maximum opening of the mouth), in the second series, the jaw was displaced vertically by 48 mm with a simultaneous transversal movement of 10 mm, in the third series, the jaw was displaced along the vertical axis down by 34 mm and transversally by 22 mm.

RESULTS

The largest distance between IAN and IAA was noted in the third series of experiments. The distance between the nerve and the vessel was minimal in the first series, with an open mouth without lateral displacements.

CONCLUSION

The generated computer model opens new possibilities for studying the dynamic anatomy of the pterygomandibular space. The results of this study can be used for further experimental and clinical trials to find the safest approach to the implementation of IANB, as well as applied in the practice of the educational process.

The effects of lingual buttons, precision cuts, and patient-specific attachments during maxillary molar distalization with clear aligners: Comparison of finite element analysis

INTRODUCTION

This study aimed to analyze the biomechanical effects of the combined use of clear aligners (CA) and auxiliaries (precision cuts, lingual buttons, and patient-specific attachments) on mesial tipping and extrusion of the premolars during maxillary molars distalization.

METHODS

Three-dimensional finite element method was employed to simulate clinical scenarios of CA with different auxiliaries for molar distalization. As such, 200 g of distal force was applied to the microimplants from the notches, lingual buttons, and hooks. Orthodontic tooth movement and the hydrostatic pressure in the periodontal ligament were compared.

RESULTS

Adding auxiliaries can provide the maxillary arch anchorage and promote the distal tipping of premolars and retroclination of maxillary incisors. In contrast, this effect was more pronounced in patient-specific attachment applications than in other types of auxiliaries. The independent application of the CA caused mesial tipping and extrusion of the premolar and also caused the incisor proclination.

CONCLUSIONS

The anchorage loss caused by the CA alone could be alleviated with the assistance of auxiliaries. Notably, patient-specific attachments further reinforce the anchorage of the anterior arch by incorporating anchor teeth as 1 anchorage unit.

Does the model reflect the system? When two-dimensional biomechanics is not 'good enough'

Models are mathematical representations of systems, processes or phenomena. In biomechanics, finite-element modelling (FEM) can be a powerful tool, allowing biologists to test form-function relationships in silico, replacing or extending results of in vivo experimentation. Although modelling simplifications and assumptions are necessary, as a minimum modelling requirement the results of the simplified model must reflect the biomechanics of the modelled system. In cases where the three-dimensional mechanics of a structure are important determinants of its performance, simplified two-dimensional modelling approaches are likely to produce inaccurate results. The vertebrate mandible is one among many three-dimensional anatomical structures routinely modelled using two-dimensional FE analysis. We thus compare the stress regimes of our published three-dimensional model of the chimpanzee mandible with a published two-dimensional model of the chimpanzee mandible and identify several fundamental differences. We then present a series of two-dimensional and three-dimensional FE modelling experiments that demonstrate how three key modelling parameters, (i) dimensionality, (ii) symmetric geometry, and (iii) constraints, affect deformation and strain regimes of the models. Our results confirm that, in the case of the primate mandible (at least), two-dimensional FEM fails to meet this minimum modelling requirement and should not be used to draw functional, ecological or evolutionary conclusions.

Macaca mulatta is a good model for human mandibular fixation research

Biomechanical and clinical studies have yet to converge on the optimal fixation technique for angle fractures, one of the most common and controversial fractures in terms of fixation approach. Prior pre-clinical studies have used a variety of animal models and shown abnormal strain environments exacerbated by less rigid (single-plate) Champy fixation and chewing on the side opposite the fracture (contralateral chewing). However, morphological differences between species warrant further investigation to ensure that these findings are translational. Here we present the first study to use realistically loaded finite-element models to compare the biomechanical behaviour of human and macaque mandibles pre- and post-fracture and fixation. Our results reveal only small differences in deformation and strain regimes between human and macaque mandibles. In the human model, more rigid biplanar fixation better approximated physiologically healthy global bone strains and moments around the mandible, and also resulted in less interfragmentary strain than less rigid Champy fixation. Contralateral chewing exacerbated deviations in strain, moments and interfragmentary strain, especially under Champy fixation. Our pre- and post-fracture fixation findings are congruent with those from macaques, confirming that rhesus macaques are excellent animal models for biomedical research into mandibular fixation. Furthermore, these findings strengthen the case for rigid biplanar fixation over less rigid one-plate fixation in the treatment of isolated mandibular angle fractures.

Assessment of labially impacted canines traction mode with clear aligners vs. fixed appliance: A comparative study based on 3D finite element analysis

Purpose: The objective of this study was to evaluate and compare the biomechanical differences between clear aligner and fixed appliance in the traction of labially impacted canines based on 3D finite element analysis.

Methods: A series of patient-oriented finite element models were constructed, including a maxillary dentition with a right labially canine, maxilla, periodontal ligaments, traction attachments, and clear aligners. The two most common clinical scenarios were investigated: Scenario A: impacted canine (distal) and Scenario B: impacted canine (mesial). For each clinical scenario, three traction models with clear aligners and one fixed appliance model were established.

Results: In all four models, the impacted canines exhibited similar initial displacement tendencies of mesially rotated in Scenario A and distally rotated in Scenario B, and with small differences in periodontal ligament stress magnitude. However, the sum of the periodontal ligament stresses of the anchorage teeth in the clear aligner mode was in the range of 56.28–76.21 kPa and in the fixed appliance mode was in the range of 6.61–7.22 kPa. The maximum value of initial displacement of the anchorage teeth in the clear aligner mode was in the range of 13.71–19.72 μm, while in the fixed appliance mode was 3.10–3.92 μm.

Conclusion: For impacted canines, clear aligner mode and fixed appliance mode have little difference in biomechanical effect. However, the anchorage teeth in the clear aligner mode endure higher stress and show a more pronounced displacement tendency. In addition, the biomechanical effects of different clear aligner traction models are various but not obvious.

Efficacy of various multi-layers of orthodontic clear aligners: a simulated study

BACKGROUND

In the invisible orthodontic treatment, composite thermoforming film materials have become the focus of orthodontic clear aligners. The orthodontic efficacy of clear aligners which consisted of multi-layers materials remains unclear. This study aims to evaluate the biomechanical effects of various multi-layers of clear aligners on en-mass retraction of maxillary anterior teeth.

METHODS

A patient-specific 3D non-linear finite element model numerical analysis was constructed to simulate the en-mass retraction of maxillary anterior teeth with clear aligner after extraction of the first premolars. Four kinds of multi-layers clear aligners with different proportion of film materials were simulated. The biomechanical responses of four different clear aligners on invisible orthodontics were calculated. The tooth displacement in all directions, the hydrostatic pressure of periodontal ligament, the orthodontic deformation of clear aligner, and the stress distribution of alveolar bone were compared and investigated.

RESULTS

In all experimental models, the maximum equivalent deformation of alveolar bone, the vector displacement of tooth and the compressive/tensile stress of periodontal ligament decreased with the increase of soft layer thickness. The elastic strain of clear aligners also decreased with the increase of the ratio of soft/hard layers.

CONCLUSIONS

The multi-layers clear aligner is better than the single-layer clear aligner in tooth movement, stress distribution of periodontal ligament and mechanical loading of alveolar bone, especially when the ratio of soft layer to hard layer is more than 50%. Moreover, the side effects of the multi-layers clear aligner are significantly less than those of the single-layer one.

Fracture Fixation Technique and Chewing Side Impact Jaw Mechanics in Mandible Fracture Repair

Lower jaw (mandible) fractures significantly impact patient health and well‐being due to pain and difficulty eating, but the best technique for repairing the most common subtype—angle fractures—and rehabilitating mastication is unknown. Our study is the first to use realistic in silico simulation of chewing to quantify the effects of Champy and biplanar techniques of angle fracture fixation. We show that more rigid, biplanar fixation results in lower strain magnitudes in the miniplates, the bone around the screws, and in the fracture zone, and that the mandibular strain regime approximates the unfractured condition. Importantly, the strain regime in the fracture zone is affected by chewing laterality, suggesting that both fixation type and the patient's post‐fixation masticatory pattern—ipsi‐ or contralateral to the fracture— impact the bone healing environment. Our study calls for further investigation of the impact of fixation technique on chewing behavior. Research that combines in vivo and in silico approaches can link jaw mechanics to bone healing and yield more definitive recommendations for fixation, hardware, and postoperative rehabilitation to improve outcomes. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Finite element simulations of the effects of an extraoral device, RAMPA, on anterosuperior protraction of the maxilla and comparison with gHu-1 intraoral device

OBJECTIVES

To investigate the effects of an extraoral device, right-angle maxillary protraction appliance (RAMPA), combined with a semi-rapid maxillary expansion intraoral device (gHu-1) on the anterosuperior protraction of maxillary bone.

MATERIALS AND METHODS

The finite element (FE) model included craniofacial bones and all sutures. The linear assumption was assumed for the FE simulations and the material properties of bones and sutures. The gHu-1 was simulated under screw activations equal to Δx = 0.25 and 0.5 mm in the lateral direction with and without RAMPA under a set of external forces {F1 = 2.94, F2 = 1.47, F3 = 4.44} N.

RESULTS

Displacement contours, nodal displacements of 12 landmarks, and von Mises stresses were compared. Combining RAMPA and gHu-1 (with Δx = 0.25 mm) resulted in changes in the displacement of the front part of the maxilla near the mid-palatal suture from (0.02, -0.1, -0.02) mm to (0.02, 0.3, 0.8) mm. For gHu-1 with Δx = 0.5 mm, the displacement of the same part changed from (0.04, -0.04, -0.2) mm to (0.04, 0.3, 0) mm. Similar trends were found in other locations.

CONCLUSIONS

The findings are in agreement with the previous cephalometric clinical data of an 8-year-old patient and prove the positive effects of RAMPA on the anterosuperior protraction of the maxilla when it is combined with the intraoral device gHu-1. In addition, RAMPA does not interfere with the lateral expansion generated by the intraoral device.

Orthodontic Tooth Movement Studied by Finite Element Analysis: an Update. What Can We Learn from These Simulations?

PURPOSE OF REVIEW

To produce an updated overview of the use of finite element (FE) analysis for analyzing orthodontic tooth movement (OTM). Different levels of simulation complexity, including material properties and level of morphological representation of the alveolar complex, will be presented and evaluated, and the limitations will be discussed.

RECENT FINDINGS

Complex formulations of the PDL have been proposed, which might be able to correctly predict the behavior of the PDL both when chewing forces and orthodontic forces are simulated in FE models. The recent findings do not corroborate the simplified view of the classical OTM theories. The use of complex and biologically coherent FE models can help understanding the mechanisms leading to OTM as well as predicting the risk of root resorption related to specific force systems and magnitudes.

Biomechanics of the mandible of Macaca mulatta during the power stroke of mastication: Loading, deformation, and strain regimes and the impact of food type

Mandible morphology has yet to yield definitive information on primate diet, probably because of poor understanding of mandibular loading and strain regimes, and overreliance on simple beam models of mandibular mechanics. We used a finite element model of a macaque mandible to test hypotheses about mandibular loading and strain regimes and relate variation in muscle activity during chewing on different foods to variation in strain regimes. The balancing-side corpus is loaded primarily by sagittal shear forces and sagittal bending moments. On the working side, sagittal bending moments, anteroposterior twisting moments, and lateral transverse bending moments all reach similar maxima below the bite point; sagittal shear is the dominant loading regime behind the bite point; and the corpus is twisted such that the mandibular base is inverted. In the symphyseal region, the predominant loading regimes are lateral transverse bending and negative twisting about a mediolateral axis. Compared with grape and dried fruit chewing, nut chewing is associated with larger sagittal and transverse bending moments acting on balancing- and working-side mandibles, larger sagittal shear on the working side, and larger twisting moments about vertical and transverse axes in the symphyseal region. Nut chewing is also associated with higher minimum principal strain magnitudes in the balancing-side posterior ramus; higher sagittal shear strain magnitudes in the working-side buccal alveolar process and the balancing-side oblique line, recessus mandibulae, and endocondylar ridge; and higher transverse shear strains in the symphyseal region, the balancing-side medial prominence, and the balancing-side endocondylar ridge. The largest food-related differences in maximum principal and transverse shear strain magnitudes are in the transverse tori and in the balancing-side medial prominence, extramolar sulcus, oblique line, and endocondylar ridge. Food effects on the strain regime are most salient in areas not traditionally investigated, suggesting that studies seeking dietary effects on mandible morphology might be looking in the wrong places.

β-catenin and Its Relation to Alveolar Bone Mechanical Deformation - A Study Conducted in Rats With Tooth Extraction

The aim of this study was to analyze the relationship between alveolar bone deformation and β-catenin expression levels in response to the mechanical load changed by dental extraction in adult rats. Twenty-four male rats (Rattus norvegicus albinus), Wistar linage, at 2 months of age, were used. The right upper incisor tooth was extracted, and euthanasia occurred in periods 5 (n = 6), 7 (n = 6), and 14 (n = 6) days after Day 0. In the control group (n = 6), the dentition was maintained. The euthanasia occurred within 14 days after day 0. After euthanasia, the rats of all groups had their left jaw with tooth removed and separated in the middle. The pieces were undergone routine histological processing and then the immunohistochemical marking were performed to label expression of the primary β-catenin antibody, which was evaluated by qualitative and quantitative analysis. One head by each group (control and experimental) was submitted to computerized microtomography. After the three-dimensional reconstruction of the skull of the rat in each group, the computational simulation for finite elements analysis were performed to simulate a bite in the incisors. In finite element analysis, the strain patterns were evaluated after the application of bite force. The results were analyzed considering the areas in which changes in the amount of deformations were detected. The action of the bite force in the experimental condition, resulted in a uniform distribution of the amount of deformations, in addition to lower amount of deformation areas, differentiating from the control group. Comparing with the control group, the levels of β-catenin signaled in the lingual bone of the middle third of the alveolar bone were raised in the periods of 5 and 14 days. The increased β-catenin positive staining intensity was concentrated on osteocytes and gaps of osteocytes. The findings of the present study were in accordance with our hypothesis that the condition of dental extraction can cause the expression of β-catenin and alter the regimes of bone deformation.

The Mechanical Effect of the Periodontal Ligament on Bone Strain Regimes in a Validated Finite Element Model of a Macaque Mandible

The primary anatomical function of the periodontal ligament (PDL) is to attach teeth to their sockets. However, theoretical and constitutive mechanical models have proposed that during mastication the PDL redistributes local occlusal loads and reduces the jaw's resistance to torsional deformations. These hypotheses imply that accurately modeling the PDL's material properties and geometry in finite element analysis (FEA) is a prerequisite to obtaining precise strain and deformation data. Yet, many finite element studies of the human and non-human primate masticatory apparatus exclude the PDL or model it with simplicity, in part due to limitations in μCT/CT scan resolution and material property assignment. Previous studies testing the sensitivity of finite element models (FEMs) to the PDL have yielded contradictory results, however a major limitation of these studies is that FEMs were not validated against in vivo bone strain data. Hence, this study uses a validated and subject specific FEM to assess the effect of the PDL on strain and deformation regimes in the lower jaw of a rhesus macaque (Macaca mulatta) during simulated unilateral post-canine chewing. Our findings demonstrate that the presence of the PDL does influence local and global surface strain magnitudes (principal and shear) in the jaw. However, the PDL's effect is limited (diff. ~200-300 με) in areas away from the alveoli. Our results also show that varying the PDL's Young's Modulus within the range of published values (0.07-1750 MPa) has very little effect on global surface strains. These findings suggest that the mechanical importance of the PDL in FEMs of the mandible during chewing is dependent on the scope of the hypotheses being tested. If researchers are comparing strain gradients across species/taxa, the PDL may be excluded with minimal effect on results, but, if researchers are concerned with absolute strain values, sensitivity analysis is required.