The generated model yielded a total of 1,854,710 tetrahedral elements with 350,289 nodes, and was imported into a FEA program ABAQUS 6.9 (ABAQUS Inc, Providence, RI). Dentine, pulp, enamel, periodontal ligament, dental follicle, cortical bone and cancellous bone were assigned relevant mechanical properties as indicated in
Table 1, according to values established in the literature [53] (
link), [63] (
link)–[66] (
link). Please note that although the mechanical properties of bone do vary according to the direction of force applied, these differences are nonetheless fairly small in all three directions [67] (
link), [68] (
link) and to a considerable extent lead to very close agreement between experimental and numerical results [64] (
link), [69] (
link). For this reason, use of linear isotropic elastic properties seems justifiable in anatomical modelling under physiological loading [70] (
link). As such, tissues were treated as homogeneous, isotropic and linearly elastic materials.
Sites for attachment of masticatory muscles were defined according to the established anatomical literature, in which the force from the following muscles was modelled: superficial and deep masseter; anterior, middle and posterior fascicles of temporalis; medial and lateral pterygoid; and the anterior belly of the digastric [71] , [72] (
link) (
Figure 3).
The temporomandibular joint is comparatively resilient, with the condyle impacting against a fibrous articular disk in the anterior region, and articular ligament material in the posterior region. To model this, two blocks of cortical bone were positioned on the articular surfaces of condyles, and the space in-between filled by a 2 mm thick layer of elastic material (
Figure 3), while the mechanical properties assigned to the anterior and posterior parts of this area were analogous to the anterior and posterior articular joint materials as indicated in
Table 1[64] (
link), [66] (
link).