Timing selection for loosened tooth fixation based on degree of alveolar bone resorption: a finite element analysis

Abstract

Objective

This study aimed to evaluate timing of fixation to retard bone absorption using finite element analysis(FEA).

Methods

Volunteer CT images were used to construct four models of mandibles with varying degrees of alveolar bone resorption. By simulating occlusal force loading, biomechanical analysis was made on the periodontal membrane, tooth root and surrounding bone (both cancellous and cortical) of mandibular dentition.

Results

The von Mises stress value of the periodontal structures was positively related with the degree of alveolar bone resorption, and the von Mises stress at the interface between the periodontal membrane and tooth root was increased significantly in moderate to severe periodontitis models. The von Mises stress at the interface between the periodontal cortical bone and cancellous bone was increased significantly in the severe periodontitis model. And the von Mises stress value with oblique loading showed significantly higher than vertical loading.

Conclusion

Teeth with moderate to severe periodontitis, loosened tooth fixation can be used to retard bone absorption.

Introduction

Periodontitis has become an oral disease with very high incidence in China, and has generally, become the leading cause of tooth loss in adults [1]. Clinically, tooth loss caused by chronic periodontitis is often accompanied by severe soft and hard tissue defects, which affects appearance and function, and also brings certain difficulties to denture restoration and implantation [2]. Periodontal disease treatment is mainly synthetic serial treatment, with periodontal drug treatment and surgical treatment to control inflammation, with restoration treatment to adjust occlusion, eliminate tooth trauma and reduce the burden of periodontal supporting tissues. With loosened tooth fixation to disperse occlusal force, control pathological loosening and displacement, and improve chewing efficiency, the affected tooth can get physiological rest, which is beneficial to the control of periodontal inflammation and healing of damaged periodontal tissue, and the oral comfort of patient was improved [3]. After perfect periodontal treatment, tooth looseness of patients with moderate to severe periodontitis can be significantly improved [4].

Due to limitations of medical ethics, it is very difficult to directly measure the intraoral biomechanical characteristics of patients with periodontitis. FEA can be used to make stress analysis of the oral structure, shape, load and mechanical properties of materials for objective, accurate and true stress distribution [5]. Baghdadi D used FEA to explain the biomechanical behaviour of crowded lower front teeth in reduced Periodontium [6]. Many in vitro studies have proved splinting was effective in slowing bone resorption and have measured tooth mobility generated by different dental splint materials to determine their flexibility using FEA [7,8,9]. However, these studies did not indicate which timepoint of splinting was more benefical. Hence, the aim of this study was to provide a scientific biomechanical basis for the timing of fixation of loosened tooth.

Materials and methods

The volunteer

The experimental volunteer was selected from the 2nd generation of Chinese virtual person (living digital-human) "No.23" volunteer, with complete dentition.

Finite element models

CT image data of mandible of the volunteer was imported into the medical three-dimensional image reconstruction software Materialise Mimics 10.0 (Materialise Company, Belgium) in DICOM format to generate the bone tissue surface profile. The gap between the tooth root and mandible was filled with a 0.2 mm periodontal membrane structure, and in combination with various software tools and clinical anatomical features, the profile curve of each structure was manually plotted. The three-dimensional digital anatomy model of mandible was reconstructed (including the separation of cortical bone and cancellous bone, enamel, dentin-cementum complex, periodontal membrane, dental pulp, temporomandibular joint condyle and other structures), then imported into reverse engineering software Geomagic Studio12.0 (Geomagic Company, USA) in STL format, using polygon editing tools to optimize the mesh of each independent anatomical structure. Alveolar bone resorption in chronic periodontitis is generally described by the ratio of resorption area to the length of tooth root, and usually divided into three degrees. In order to unify bone resorption of each tooth position in the sample models, the critical value of bone resorption in varying degrees of periodontitis was used as the standard for division. The model was trimmed, and the total length of the measured tooth root was subtracted from the distance of enamel-cementum boundary and alveolar ridge to calculate the root length of the inner alveolar bone of each tooth. The height of alveolar bone was reduced to 0, 1/3, 1/2 and 2/3 of the inner length of the root bone. The edge was trimmed to make it smooth and continuous. Four models were established to simulate (a) non-periodontitis (b) mild periodontitis (c) moderate periodontitis and (d) severe periodontitis of the mandible models, and to simulate the stable state of bone resorption after periodontal treatment. The model was imported into the finite element pre-processing software HyperMesh 11.0 (Altair Company, USA) for materialization. After the Boolean operation, finite element models of a mandible with varying degrees of alveolar bone resorption were constructed (Fig. 1).

Fig. 1

Modeling Steps: A CBCT image of living digital-human No. 23. B the three-dimensional digital anatomy model of mandible. C The mesh refinement of mandible. D The mesh refinement of enamel. E The mesh refinement of dentin-cementum complex. F The mesh refinement of dental pulp. G The mesh refinement of periodontal membrane

Mesh independency analysis

Take the non-periodontitis model vertical loading as an example. The initial mesh size of teeth and surrounding tissues adopted 0.5 mm, and the maximum stress error percentage between 0.5 mm and 0.25 mm was 1.4% (Table 1). Given the calculation efficiency, the four modles (Fig. 2) were established using 0.5 mm mesh size. There were 283,886 tetrahedrons and 53,718 nodes in the non-periodontitis model, 274,913 tetrahedrons and 52,382 nodes in the mild periodontitis model, 273,965 tetrahedrons and 52,719 nodes in the moderate periodontitis model, and 261,095 tetrahedrons and 50,911 nodes in the severe periodontitis model.

Table 1 Mesh independency analysis

Fig. 2

Perspective view of finite element models of mandibles with varying degrees of periodontitis. Note In all the charts in this study, A normal mandible; B mandible with mild periodontitis; C mandible with moderate periodontitis; D mandible with severe periodontitis

Experimental preset conditions

All materials were considered linear-elastic, isotropic and homogeneous [10,11,12,13,14,15,16]. The applied material properties (elastic modulus and Poisson's ratio) were obtained from the literature (Table 2).

Table 2 Material mechanical parameters

Loading force

Different teeth play diverse roles in chewing function, same loading was applied to the teeth of same region. Therefore, the loading of 100 N, 50 N and 25 N were loaded on the occlusal surfaces of molars, premolars, and anterior teeth, respectively [10]. In this study, a certain distributed point on the occlusal surface of each tooth was selected for concentrate loading [17], and the chewing cycle was simplified to vertical loading and oblique loading. Under loading conditions, the edge of the mandible was fixed to prevent the jaw from displacement in the X, Y and Z directions [13, 15] (Fig. 3).

Fig. 3

A 1. The point of occlusal loading, 2. The fixed edge. B 3. Vertical loading. C 4. Oblique loading

Results

Biomechanical analysis of tooth root interface

During vertical loading, the von Mises stress distribution on the root surface was concentrated on the root labial side of the lower anterior teeth. The maximum Von Mises stress (MVMS) of mandibular dentition was concentrated on the incisor. The MVMS of the root interface of normal mandible was 13.43 MPa; the MVMS of the root interface of mild periodontitis was 15.60 MPa; the MVMS of the root interface of moderate periodontitis was 21.79 MPa; the MVMS of the root interface of severe periodontitis was 29.97 MPa (Fig. 4, Table 3). During 45° oblique loading, the stress was concentrated on the root surface of mesiolingual posterior teeth, and the MVMS was concentrated on the root surface of the mandibular first molar. The stress was 2–3 times than when vertical loading (Fig. 5, Table 4). In addition, during both vertical and oblique loading, the MVMS was increased with aggravation of periodontitis, and stress was increased more significantly with moderate to severe periodontitis.

Fig. 4

Cloud chart of a finite element analysis of mandibular root interface during vertical loading. A normal mandible; B mandible with mild periodontitis; C mandible with moderate periodontitis; D mandible with severe periodontitis

Table 3 Peak stress of root interface during vertical loading (MPa)

Fig. 5

Cloud chart of a finite element analysis of mandibular root interface during 45° oblique loading. A normal mandible; B mandible with mild periodontitis; C mandible with moderate periodontitis; D mandible with severe periodontitis

Table 4 Peak stress of root interface during 45° oblique loading (MPa)

Biomechanical analysis of periodontal membrane interface

During vertical loading, the von Mises stress distribution of the periodontal membrane was concentrated at the corresponding position of periodontal membrane on the labial side of the lower anterior tooth root close to the neck. The MVMS was concentrated in the incisor. The MVMS of the periodontal membrane of a normal mandibular was 1.35 MPa; the MVMS of the periodontal membrane with mild periodontitis was 1.65 MPa; the MVMS of the periodontal membrane with moderate periodontitis was 2.33 MPa; the MVMS of the periodontal membrane with severe periodontitis was 3.15 MPa (Fig. 6, Table 5). During 45° oblique loading, stress was concentrated at the corresponding position of the periodontal membrane on the lingual side of the lower anterior tooth root close to the neck. The MVMS was concentrated on the mandibular first molar, and the stress was 2–4 times than during vertical loading (Fig. 7, Table 6). It was found that during vertical and oblique loading, the MVMS increased with the aggravation of periodontitis, and stress of the periodontal membrane was increased more significantly in moderate to severe periodontitis.

Fig. 6

Cloud chart of a finite element analysis of a periodontal membrane of mandibular dentition during vertical loading. A normal mandible; B mandible with mild periodontitis; C mandible with moderate periodontitis; D mandible with severe periodontitis

Table 5 Peak stress of periodontal membrance during vertical loading (MPa)

Fig. 7

cloud chart of a finite element analysis of a periodontal membrane of mandibular dentition during 45° oblique loading. A normal mandible; B mandible with mild periodontitis; C mandible with moderate periodontitis; D mandible with severe periodontitis

Table 6 Peak stress of periodontal membrance during 45° oblique loading (MPa)

Biomechanical analysis of the periodontal cortical bone interface

During vertical loading, occlusal force was transmitted to the mandible via the periodontal membrane of the tooth. When the cortical bone of the mandible was stressed, the maximum Von Mises stress was located in the corresponding cortical bone area of the labial neck of the incisor root; the minimum stress was located in the corresponding cortical bone area of the neck of the canine root (Fig. 8, Table 7). With 45° oblique loading, stress was concentrated in the corresponding cortical bone area of the mesiolingual neck of the posterior tooth (Fig. 9, Table 8). It was found that during vertical and oblique loading, the MVMS was increased with aggravation of periodontitis. In mild to moderate periodontitis, stress was increased slowly; in severe periodontitis, cortical bone stress was increased more significantly.

Fig. 8

Cloud chart of a finite element analysis of a periodontal cortical bone of mandible during vertical loading. A normal mandible; B mandible with mild periodontitis; C mandible with moderate periodontitis; D mandible with severe periodontitis

Table 7 Peak stress of periodontal cortical bone during vertical loading (MPa)

Fig. 9

Cloud chart of a finite element analysis of a periodontal cortical bone of mandible during 45° oblique loading. A normal mandible; B mandible with mild periodontitis; C mandible with moderate periodontitis; D mandible with severe periodontitis

Table 8 Peak stress of periodontal cortical bone during 45° oblique loading (MPa)

Biomechanical analysis of the periodontal cancellous bone interface

When the cancellous bone of the mandible was stressed, the stress was mainly distributed in the corresponding cancellous area of the root apex. During vertical and oblique loading, the MVMS was increased with aggravation of periodontitis. The maximum Von Mises stress was located in the apical area of the second molar. In mild to moderate periodontitis, stress increased slowly. In severe periodontitis, stress of cancellous bone increased more significantly (Figs. 10, 11, Tables 9, 10).

Fig. 10

Cloud chart of a finite element analysis of a cancellous bone of mandible during vertical. A normal mandible; B mandible with mild periodontitis; C mandible with moderate periodontitis; D mandible with severe periodontitis

Fig. 11

Cloud chart of a finite element analysis of a cancellous bone of mandible during 45° oblique loading. A normal mandible; B mandible with mild periodontitis; C mandible with moderate periodontitis; D mandible with severe periodontitis

Table 9 Peak stress of periodontal cancellous bone during vertical loading (MPa)

Table 10 Peak stress of periodontal cancellous bone during 45° oblique loading (MPa)

Discussions

• 1. According to experimental results and Tables three to 10, the stress of each structure of periodontal tissue was increased with the degree of alveolar bone resorption. The change in stress of the periodontal membrane of the mandible was similar to the trend of change in the stress of the root surface. In moderate to severe periodontitis, stress was increased significantly, while the stress distribution of cortical bone and cancellous bone of mandible were similar. In severe periodontitis, stress increased significantly. Although the load was within the limit of the physiological tolerance of each tissue, stress fatigue would occur if a repeated load exceeded a certain limit. This result was consistent with alveolar bone resorption and tooth loosening in advanced periodontitis, and the influence of oblique force on stress distribution was greater than that of vertical force [18]. Reddy used the finite element method to analyze changes in the stress of the periodontal membrane of single maxillary incisor at different levels of alveolar bone and found that with decrease in the height of alveolar bone, stress concentration of tooth cusp increased significantly [19]. Some scholars found the maximum tensile stress of the periodontal membrane of a maxillary incisor in the lingual neck area increased with decrease in bone height, but stress decreased with the increase in width of the periodontal membrane, and changes in the height of alveolar bone had little effect on the cortical bone [20]. Based on previous studies, combined with our study, the stress distribution of the labial cortical bone and root in the anterior tooth area of the mandible were concentrated, and the increasing trend was significant with the absorption of alveolar bone. We found that more studies on splinting fixation of the lower anterior teeth were conducted in vivo and in vitro experiments than molars, which was consistent with clinical application. The finite element method is not only of great significance for study of periodontal tissue stress and deformation under loads, it also provides medical and biomechanical experts with a basis for their biomechanical behavior in oral treatment [21].

• 2. With alveolar bone resorption, the clinical crown is lengthened, which increases the axial force of the root surface. According to the results of finite element analysis, stress increase during oblique loading is more obvious than during vertical loading. Lateral force shall be minimized. Many scholars have used loosened tooth fixation to disperse the occlusal force of lower anterior teeth. Nikolaus used the finite element method to make stress analysis on the therapeutic effects of periodontal splints made of different materials [22], and many scholars have also made clinical evaluation on the treatment timing of loosened tooth fixation [23, 24]. Factors such as wear of occlusal surface and tooth movement in patients with periodontitis will affect the distribution of occlusal force and occlusal contact [25]. However, many studies have directly observed intraoral periodontal parameters for comparative analysis. In combination with study results, for the root surface and periodontal membrane interface, stress increase is particularly significant in moderate to severe periodontitis. In order to ensure uniform stress distribution and reduce stress concentration, and prevent lower anterior teeth from loosening and falling out, patients with moderate to severe periodontitis may undergo preventive internal fixation of lower anterior teeth to disperse occlusal force.

• 3. The mandible model of periodontitis in the study simulates bone level through alveolar ridge trimming. The mandibular alveolar ridge at the defect site has no obvious structural change except for loss of height. Under actual situations, alveolar bone resorption in periodontitis is not uniform and continuous horizontal resorption, and the depth of periodontal pockets on different sides of each tooth is not the same. According to literature reports, there are differences in the morphology and characteristics of a normal mandible and alveolar ridge trophic mandible, and bone characteristics and evolution, as well as changes in bone density shall also be taken into account for mandible changes [26, 27]. With rapid development of computer medical imaging technology and medical technology, application of the finite element method in the field of stomatology can achieve more accurate biomechanical analysis in vitro for various oral diseases, thereby simulating and guiding clinical practice [22].

  Many scholars [4, 28] assessed 162-227 splinted teeth, splint survival and number of plint-repairs retrospectively over 10 years to indicated that splinting can assist the retention of periodontally affected teeth by reducing their mobility. The significance of biomechanical analysis of mandibular models with varying degrees of alveolar bone resorption was to conduct mechanical analysis, and evaluation of the tooth and periodontal tissues of patients with periodontitis, to avoid stress concentration, retard bone absorption and then to prompt clinical optimization treatments, improve long-term survival rate of teeth.

• 4. Finite element analysis has some limitations that should be considered. Contact boundaries, isometric characters of materials, and static loading jeopardize the validity of analysis. Nevertheless, the simplifications adopted in this study are not expected to fundamentally alter the behavior of bone materials, and therefore, this analysis should help advance a better understanding periodontal changes of alveolar bone resorption.

• 5. As for future studies, to validate the results obtained in this simulation, further clinical and longitudinal follow- up studies are required. And it should be detail the biological behavior of all the materials to build vitro models closer to the clinic.

Conclusion

Teeth with moderate to severe periodontitis, loosened tooth fixation can be used to retard bone absorption based on the Finite Element Analysis, especially for lower anterior teeth.