Effects of Condylar Elastic Properties to Temporomandibular Joint Stress
Mandibular condyle plays an important role in the growth and reconstruction of the temporomandibular joint (TMJ). We aimed to obtain orthotropic elastic parameters of the condyle using a continuous-wave ultrasonic technique and to observe the effects of condylar elastic parameters on stress distribution of the TMJ using finite element analysis (FEA). Using the ultrasonic technique, all nine elastic parameters were obtained, which showed that the mandibular condyle was orthotropic. With the condyle defined as orthotropic, the occlusal stress was transferred fluently and uniformly from the mandible to the TMJ. The stress distribution in the isotropic model showed stepped variation among different anatomical structures with higher stress values in the cartilage and condyle than in the orthotropic model. We conclude that anisotropy has subtle yet significant effects on stress distribution of the TMJ and could improve the reality of simulations.
Stresses in the temporomandibular joint (TMJ) are important for maintaining normal structure and function . Finite element analysis (FEA) is often used in stress analysis of the TMJ. Two-dimensional FEA of the mandible was developed by Haskell et al.  and Maeda et al.  to investigate stress distribution in the TMJ. Three-dimensional models were subsequently developed to simulate the human mandible, including the TMJ, more accurately [4–8]. Some of these models were used to investigate the influences of craniofacial discrepancy  and the relationship between disc displacement and TMJ stress distribution [10–12]. Although FEA is a promising tool in the study of functional morphology of craniofacial structures, its reliability depends greatly on the accuracy of the input, that is, parameters such as load, geometrical form, and elastic properties.
A series of studies have shown that if bone properties are regarded as isotropic, it may result in misinterpretations of strain data, for bony structures of identical shape will resist loads differently if their elastic properties differ [13, 14]. However, in FEA of the TMJ, the condyle is usually assumed to be isotropic and defined with the same elastic properties as the mandible. Mechanical testing, however, has shown that the cancellous bone of the condyle is anisotropic [15, 16] and consists of parallel plate-like trabeculae primarily oriented in the vertical direction, perpendicular to the mediolateral condylar axis. In the horizontal direction, the plate-like trabeculae are interconnected with rods . To our knowledge, no studies have determined all isotropic elastic parameters of the mandibular condyle.
Many researchers use mechanical testing to obtain elastic parameters of biological tissues [18, 19]. However, mechanical testing can lead to nonuniform loading at the specimen ends due to surface irregularities, induced bending, and difficulties associated with measuring accurately. Comparatively, ultrasonic techniques offer advantages over mechanical techniques [20, 21]. Smaller specimens can be used, and anisotropic properties can be determined from a smaller region of the whole bone. Until now, ultrasonic techniques have been used to determine elastic properties of human and canine cortical femora , mandible , cancellous bone , cortical bone , dentate mandible , enamel and dentin , and craniofacial skeletons of primates . However, there are no reports on the application of ultrasonic techniques to the mandibular condyle.
We investigated the orthotropic elastic properties of the condyle with continuous-wave ultrasonography and then applied findings to the FEA model. Our objective was to determine the orthotropic parameters of the condyle and further to prove the effects of the condylar anisotropy to the level and distribution of the TMJ stress.
2. Materials and Methods
2.1. Ultrasonic Measurement of Condylar Elastic Properties
The freshness of samples is important for elastic properties measurement . Because of the limitations of obtaining fresh human condyles, this study used sheep joints. Compared with other animals, the size, anatomy, function, and joint accessibility of sheep are similar to humans  (Figure 1(a)). Six fresh condyles from sheep (aged months) were dissected and measured. Care was taken throughout the procedure to maintain the specimens in a moist state. The specimens were frozen in saline for preparation. Before measurement, the condyles were taken out and cut into -mm standard samples (Figure 1(b)). After the specimens were prepared, they were stored in a solution of 50% ethyl alcohol and 50% normal saline to maintain the elastic properties of the tissue . The density of each specimen was calculated based on Archimedes’ principle of buoyancy. On the samples, the sides were marked to denote the orientation with respect to the three-dimensional coordinate system, in which was the antero-posterior direction, was the vertical direction, and was the mesio-lateral direction.
We used an ultrasonic pulse transmission to measure the elastic properties [14, 23]. Ultrasonic waves were generated with a pulse generator, and two sets of mounted piezoelectric transducers were used (1.25 MHz longitudinal; 2.5 MHz transverse, Valpey-Fischer, Hopkinton, MA). Both longitudinal and transverse ultrasonic waves were passed through various axes (Figures 1(c) and 1(d)). Time delays were measured using a digital storage oscilloscope (54601-A, Hewlett-Packard, Austin, Texas) to make a phase comparison of the signal before and after its transmission through a specimen. Ultrasonic velocities were calculated by dividing the specimen thickness by the apparent time delay minus the system time delay. Organic glass with known ultrasonic velocity was used for system verification. The absolute error of the system was 0.11 m/s. Five measurements were performed for each specimen. Thus, 9 ultrasound velocities transmitted in the condyle were obtained, in which , , , , , and were longitudinal wave speed measured with the generation and sensing transducer in line, and ,, and were the transverse wave speed measured with the generation and sensing transducer at a right angle. Here, 1 indicated the antero-posterior direction, 2 indicated the vertical direction, and 3 indicated the mesio-lateral direction. Linear elastic wave theory yields relationships between various velocities through the specimen, and its elastic properties was as follows: , , , , , , , , , from which all the 9 independent elastic coefficients of the condyle were calculated. Further, according to the matrix relationship between elastic coefficients and technical constants expressed as follows , Young’s modulus , shear modulus , and Poisson’s ratio in each direction were obtained, which were further applied to FEA for defining the orthotropic material parameters of the condyle:
2.2. Construction of the FEA Model
A young male adult cadaver with normal interjaw relationship and natural dentition was used to construct the FEA model. The use of a cadaver conforms to a written protocol that was reviewed and approved by the Department of Anatomy of the Fourth Military Medical University. Multiplanar computerized tomography scan (Tomoscan SR 7000, Philip, Eindhoven, the Netherlands) images of the bone structure with 1.0 mm sections were obtained for hard tissue reconstruction. One-millimeter manually sliced and digitalized images of bilateral discs of the same cadaver were used for soft tissue reconstruction. A layer of 2 mm elements was established to simulate articular cartilage. Then, the tetrahedral elements with three degrees of freedom per node were adopted to build superficial and internal structure meshes for the model. During that process, the meshes in mandible were relatively sparse and those of condyle were densely plotted for it was one of the pivots in this research. Thus, a three-dimensional FEA model was established by ALGOR (Super Sap) software (ALGOR Inc., Plano, Texas. At the end of this process, the FEA model presented a mesh comprising 8292 elements and 10936 nodes (Figure 2(a)).
Superior disc surface and the mandibular angle area were restricted. Two different models with the same geometry, mesh, boundary, and loading conditions but different condylar material properties, one with isotropic and the other with orthotropic elastic parameters, were established. Isotropic material properties were homogeneous in all directions of the element. Orthotropic material had three material axes, which corresponded to those in ultrasonic measurements. The orientations of each axis in the orthotropic element are shown in Figure 2(b). Each of the three material axes was assigned a value for the elastic module, and each of the three planes defined by those axes received values for the shear modulus and Poisson’s ratio. The elastic parameters of the mandible, disc, and cartilage were determined by referring to previous studies [5, 10, 30]. For the load condition, stress vectors previously measured by us with a photoelastic method  were applied to each mandibular tooth position.
During observation of stress values, each of the three functional parts was divided into nine parts from the medial to the lateral part and from the anterior to the posterior part. Twenty-nine observation points were located (Figure 2(c)), in which point number 1 and number 2 were in the medial and lateral part of the condylar neck, number 3 to number 11 were in the condyle, number 12 to number 20 were in the cartilage, and number 21 to number 29 were in the disc. Maximal principal stress and minimal principal stress were considered to quantitatively compare the effect of anisotropy of the condyle.
2.3. Statistical Analysis
Elastic coefficients obtained were analyzed by analysis of variance. Dunnett's post hoc test was used to test for significant differences between directions for the elastic moduli, shear moduli, and Poisson’s ratios. TMJ stress values of a group of observation points in the condyle, cartilage, or disc between the isotropic and orthotropic models were compared by the between-subjects effects test of univariate analysis of variance with significance set at 0.05 using SPSS 11.5 software (SPSS Inc., Chicago, Illinois).
3.1. Elastic Parameters of Condyle
Ultrasound velocities transmitted in the condyle (), the elastic coefficients (), and the technical constants (, and ) that were applied to the FEA are shown in Table 1. All 9 orthotropic elastic coefficients of the condyle were successfully obtained from each small condyle sample. The results showed significant differences of elastic coefficients between directions, such as () and (), which indicated that the mandibular condyle was orthotropic material. The values of condylar elastic parameters (see Table 1) were different from that of cortical bone ( MPa, ) and cancellous bone ( MPa, ) . Young’s modulus of condyle in (vertical to transverse axis) direction, in which the condyle endures the primary force in functional movement, showed the highest value For the shear modulus, and had the highest value, which meant that the condyle has high endurance for shearing or twisting loads along the direction parallel to the transverse ridge.
3.2. Stress Distribution of TMJ
A three-dimensional FEA model, including the mandible, condyle, disc, and cartilage, was established in this study. The distributions of the maximal and minimal principal stress in the model with the orthotropic condyle are shown in Figures 3(a) and 3(b). After loading of occlusal force, the stress was transferred fluently and uniformly from the mandible to the TMJ. Stress of the anterior part of the TMJ was mainly compressive and that of posterior part was mainly tensile in the intercuspal position. The maximal and minimal principle stress in the condyle and cartilage of the TMJ model with the orthotropic condyle was lower than that of the model with the isotropic condyle (observation points 1–20, ), especially in the condylar neck and condyle area (observation points 1–11, ), whereas the stress level in the disc of both models was not significantly different (observation points 20–29, , Figures 3(c) and 3(d)). Thus, the stress distribution in the isotropic model showed stepped variation in different anatomic structures.
(c) Comparison of maximal principle stress in TMJ
(d) Comparison of minimal principle stress in TMJ
4.1. Elastic Property of the Condyle
The anatomy for sheep condyle in our research showed several characteristics that resembles to those of human, including diameter in mesial-lateral direction being longer than that in anterior-posterior direction, the long axis of bilateral condyle directing backwards, and bigger anterior condylar slope. Architecturally, the condyle also had a thin external cortex that envelops a central medullary component that is made up of trabecular bone and soft tissue marrow. The only difference was that the fossa and condyle of the sheep were flat, much like an adapted edentulous human TMJ. The condyles of 1-year-old sheep adopted in this research almost had the same size with that of adult human being. There were also thin layer of fibrocartilage covered on the condylar surface. Before the ultrasonic measurement, fibrocartilage layer had been carefully divested to ensure that elastic parameters obtained were indeed from the bone tissue of condyle. In a much detailed report on comparative anatomy research on TMJ of the goat and sheep , it was ponited out that the TMJ of goat and sheep showed many charicteristics just as those of human, such as a distinct absence of protection against anteroposterior forces, a medial and lateral ligament attaching the disc to the condyle, and the translatory capacity of the condyle. As a result, it was concluded that the goat and sheep should be strongly considered as the prime models for TMJ research and procedural training [32, 33]. Besides, it was also reported that the mean trabecular thickness in young sheep ( months) was 0.28 mm, which corresponded well with that of human (age 30–49, 0.27–0.34 mm) . Whereas, we also noticed that strains in the mandible are caused by muscle contraction and occlusal force. The conditions of isognathy versus anisognathy are expected to produce different patterns of mandibular strain during mastication. In sheep, the TMJ affords greater lateral movement of the mandible, because only one tooth row can occlude at a time, and the bite point acts as a fulcrum, which would inevitably lead to larger and differently oriented strains on the working side compared to the balancing side. The impacts of this different strain distribution character on the elastic parameters of the condyle still need further investigation.
The ultrasonic technique for measuring elastic parameters offers particular advantages over mechanical techniques, including application to small specimens and independence from friction and damage due to gripping specimens, making repetitive measurements on a single specimen possible. Thus, it has the potential to characterize three-dimensional elastic anisotropy from a single sample, which effectively reduces individual errors in measurements . Using this technique, the values of condylar elastic parameters were shown to be different from those of cortical bone and cancellous bone of the mandible . This finding could probably be explained by the different functional and remodeling characteristics of different bones . Our results showed that the highest Young’s modulus of condyle was in direction (vertical to its transverse axis). Young’s modulus in the axial direction was also shown to be higher than that of the transverse direction using mechanical tests of the human condyle .
4.2. Effects of Condylar Elastic Properties to TMJ Stress
Anisotropy of the mandible has been shown to affect interface stresses and peri-implant strain . In addition, FEA can be adversely affected when elastic properties are modeled imprecisely, and modelers should attempt to obtain elastic properties data about the species and skeletal elements that are the subject of their analyses . Our finite element calculations appeared to be the first to include regional anisotropy in TMJ stress analysis. We found that anisotropy of the condyle caused stress levels of the TMJ to decrease in the cartilage and condyle area. The stress in the TMJ, with the characteristic elastic parameters in the present study, was distributed uniformly in the condyle and transferred smoothly from the condyle to other parts of the TMJ. When the condyle was defined with the same parameters as the mandible, the stress values of the condyle were close to that of the condylar neck but much higher than those of the disc and cartilage. Thus, it seemed that the stress was not transferred continuously and gradually from the mandible to the condyle and further to the cartilage and disc, but rather descended in a step fashion.
It was approved by the present study that the mandibular condyle was orthotropic, with all nine elastic parameters having been obtained by the continuous-wave ultrasonic technique. The anisotropy character of the condyle has subtle yet significant effects on stress distribution of the TMJ, which could improve the simulation reality of finite element modeling.
This work was sponsored by the Nature Science Foundation of China (no. 39370915, 30670518) and institutionally supported by the Department of Oral Anatomy and Physiology of School of Stomatology of Fourth Military Medical University. The authors thank the acoustic research institute of Shaanxi Normal University and Professor Shang Zhiyuan for technical support and generous assistance in performing the ultrasonic tests.
V. C. Mow and A. F. Mak, “Lubricant of diarthrodial joints,” in Handbook of Bioengineering, R. Skalak and S. Chien, Eds., pp. 1–34, McGraw-Hill, New York, NY, USA, 1988.View at: Google Scholar
B. Haskell, M. Day, and J. Tetz, “Computer-aided modeling in the assessment of the biomechanical determinants of diverse skeletal patterns,” American Journal of Orthodontics, vol. 89, no. 5, pp. 363–382, 1986.View at: Google Scholar
Y. Maeda, M. Sogo, S. Tsutsumi, M. Okada, and T. Nokubi, “Biomechanical study of temporomandibular joint on its form and function—part I: condyle morphology in frontal section,” The Journal of Osaka University Dental School, vol. 33, no. 1, pp. 65–69, 1993.View at: Google Scholar
T. W. P. Korioth and A. G. Hannam, “Effect of bilateral asymmetric tooth clenching on load distribution at the mandibular condyles,” The Journal of Prosthetic Dentistry, vol. 64, no. 1, pp. 62–73, 1990.View at: Google Scholar
J. Chen and L. Xu, “A finite element analysis of the human temporomandibular joint,” Journal of Biomechanical Engineering, vol. 116, no. 4, pp. 401–407, 1994.View at: Google Scholar
K. Nagahara, S. Murata, S. Nakamura, and T. Tsuchiya, “Displacement and stress distribution in the temporomandibular joint during clenching,” Angle Orthodontist, vol. 69, no. 4, pp. 372–379, 1999.View at: Google Scholar
M. C. Castaño, U. Zapata, A. Pedroza, J. D. Jaramillo, and S. Roldán, “Creation of a three-dimensional model of the mandible and the TMJ in vivo by means of the finite element method,” International Journal of Computerized Dentistry, vol. 5, no. 2-3, pp. 87–99, 2002.View at: Google Scholar
K. Tanne, E. Tanaka, and M. Sakuda, “Stress distributions in the TMJ during clenching in patients with vertical discrepancies of the craniofacial complex,” Journal of Orofacial Pain, vol. 9, no. 2, pp. 153–160, 1995.View at: Google Scholar
B. Buranastidporn, M. Hisano, and K. Soma, “Articular disc displacement in mandibular asymmetry patients,” Journal of Medical and Dental Sciences, vol. 51, no. 1, pp. 75–81, 2004.View at: Google Scholar
E. B. W. Giesen and T. M. G. J. van Eijden, “The three-dimensional cancellous bone architecture of the human mandibular condyle,” Journal of Dental Research, vol. 79, no. 4, pp. 957–963, 2000.View at: Google Scholar
L. J. van Ruijven, E. B. W. Giesen, M. Farella, and T. M. G. J. van Eijden, “Prediction of mechanical properties of the cancellous bone of the mandibular condyle,” Journal of Dental Research, vol. 82, no. 10, pp. 819–823, 2003.View at: Google Scholar
S. S. Kohles, J. R. Bowers, A. C. Vailas, and R. Vanderby Jr., “Ultrasonic wave velocity measurement in small polymeric and cortical bone specimens,” Journal of Biomechanical Engineering, vol. 119, no. 3, pp. 232–236, 1997.View at: Google Scholar
R. B. Ashman, S. C. Cowin, W. C. van Buskirk, and J. C. Rice, “A continuous wave technique for the measurement of the elastic properties of cortical bone,” Journal of Biomechanics, vol. 17, no. 5, pp. 349–361, 1984.View at: Google Scholar
R. B. Ashman and W. C. Van Buskirk, “The elastic properties of a human mandible,” Advances in Dental Research, vol. 1, no. 1, pp. 64–67, 1987.View at: Google Scholar
F. Linde, “Elastic and viscoelastic properties of trabecular bone by a compression testing approach,” Danish Medical Bulletin, vol. 41, no. 2, pp. 119–138, 1994.View at: Google Scholar
W. C. van Buskirk, S. C. Cowin, and R. N. Ward, “Ultrasonic measurement of orthotropic elastic constants of bovine femoral bone,” Journal of Biomechanical Engineering, vol. 103, no. 2, pp. 67–72, 1981.View at: Google Scholar
M.-Q. Wang, M. Zhang, and J.-H. Zhang, “Photoelastic study of the effects of occlusal surface morphology on tooth apical stress from vertical bite forces,” Journal of Contemporary Dental Practice, vol. 5, no. 1, pp. 74–93, 2004.View at: Google Scholar
C. Bifano, G. Hubbard, W. Ehler, and A. N. Goss, “A comparison of the form and function of the human, monkey, and goat temporomandibular joint,” Journal of Oral and Maxillofacial Surgery, vol. 52, no. 3, pp. 272–277, 1994.View at: Google Scholar
J.-I. Ishimaru and A. N. Goss, “A model for osteoarthritis of the temporomandibular joint,” Journal of Oral and Maxillofacial Surgery, vol. 50, no. 11, pp. 1191–1195, 1992.View at: Google Scholar
A. M. O'Mahony, J. L. Williams, and P. Spencer, “Anisotropic elasticity of cortical and cancellous bone in the posterior mandible increases peri-implant stress and strain under oblique loading,” Clinical Oral Implants Research, vol. 12, no. 6, pp. 648–657, 2001.View at: Google Scholar